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I.  INTRODUCTION  AND  OUTLINE 

1 . 1 HISTORICAL  BACKGROUND 

The  seventies  has  ushered  in  the  era  of  computer  communication  networks  that  are  essentially  networks 
of  computers  interconnected  for  the  purpose  of  resource  sharing.  The  concept  of  reso  ce  sharing,  wherein 
users  can  access  the  vast  software  and  computer  hardware  on  a global  basis  is  a fitting  extension  of  the  time- 
sharing systems  of  the  sixties.  The  unrelenting  pursuit  of  more  powerful,  faster  and  versa  :iie  computer  capa- 
bility has  resulted  in  the  union  of  time-shared  systems  sharing  each  other’s  resources  by  interconnecting  them 
over  a network  to  effect  a computer  communication  network. 

Time-shared  systems  are  not  in  the  process  of  being  phased  out,  but  rather  are  being  expanded  depending 
on  users  need.  Computer  systems  for  airline  reservations,  hotel  reservations,  military  fire  control,  credit  check- 
ing, on-line  banking  and  savings,  etc.,  are  in  effect  groups  of  remote  terminals  feeding  a central  processor,  where 
each  CPU  is  dedicated  to  a particular  problem.  The  quest  for  bigger,  fastei,  and  more  reliable  computers  is 
asymptotically  bounded  by  the  speed  jf  an  election  thru  a conduetoi ; hence,  CPU  access  times  or  cycle  times 
in  nano-seconds  is  the  upper  limit  and  physical  size  of  the  computer  becomes  somewhat  moot. 

Space  technology  paved  the  way  for  miniaturization  and  high  reliability  by  virtue  of  micro  logic.  The 
direct  consequence  of  micro  logic  in  the  computer  field  was  the  mini-computer;  this  increased  considerably 
the  volume  of  computers  along  with  their  associated  users  and  subsequently  increased  the  demand  for  software. 
The  demand  for  more  sophisticated  software,  taxed  the  capability  of  the  mini-computer;  hence  another  bound- 
ary was  being  approached,  that  being  a given  cotnputei  complex  o *d  only  be  efficient  for  certain  classes  of 
problems,  whether  the  complex  was  made  up  of  mini  or  maxi  computers  or  a combination  of  both.  The  prob- 
lem of  effecting  a computer  facility  w ith  maximum  throughput,  minimum  cost,  high  reliability,  remote  acces- 
sing, computational  versatility,  computational  efficiency,  etc  , was  rapidly  approaching  an  impasse.  Dr.  L.  G. 
Roberts IR067)  was  initiating  the  concept  that  existing  computer  systems  could  be  better  utilized  when  inter- 
connected in  a network  so  that  the  resources  (hardware,  software,  and  data)  could  be  shared. 

The  construction  of  a computet  netwoik  enabled  each  computer  system  oriented  to  some  specialized 
effort  to  be  interconnected  on  a network,  but  each  system  maintains  its  autonomy.  The  result  of  this  inter- 
connecting of  computer  systems  saves  duplication  of  effort  and  hopefully  minimizes  cost  to  the  user,  while 
maximizing  his  computational  efficiency. 

The  computer  communication  network  while  heralding  in  the  era  of  resource  sharing  brings  with  it  many 
new  problems  associated  with  computer  networks.  Computer  network  problems  reduce  essentially  to  the  prob- 
lem of  moving  messages  from  source .9  to  destination  D in  a minimum  time  frame  (error  free)  and  is  very  simi- 
lar t<'  a transportation  or  tiaveling salesman  problem. 


I'he  computet  network  should  look  transparent  to  the  h.dividuul  computer  system  and  provide  the 

necessary  message  ham capability  reqt tired  to  transmit  and  insure  rapid  delivery  of  error  free  messages 

Itotn  source  to  destination.  Since  computer  networks  are  not  fully  connected  for  economic  consideration, 
there  exists  a requirement  for  storing  and  forwarding  of  messages  between  relay  nodes  within  the  network 
and  as  a consequence  a routing  strategy  for  rapid  delivery  of  messages  that  can  adapt  to  network  traffic  loads 
. or  blockage  must  be  established.  Thus,  the  problem,  simply  staled  but  not  simply  solved,  is  the  storing  and 
lor  warding  of’  messages (puckotN)  along  some  given  network  configuration, 

A historical  icvicw  of  the  evolution  of  computer  networks  is  contained  in  the  work  by  l;ultz  |FU72| 

ami  MIS  bibliography  |N'HS7.1j.  I,  will  suffice  for  the  present  to  use  the  definition  by  Cole  |C07I  J for  a 

"compute,  network":  "A  computer  network  is  a set  of  interconnected  processors  which  can  be  utilized  joint- 

• |,r",l"tfllw  m:"",or'  b"'  w',idl  n»nn*illy  «rc  controlled  by  separate  operating  systems,  and  can  perform 
in  an  autonomous  mat  met 


llic  uiKilylic  heat  men  I o!  stochastic  How  of  message  traffic  in  connected  networks  of  communication 
w;is  si i it  1 1 c cl  by  Klein mck  |KLMJ.  Klcinmck  considered  i etworks  that  were  channel-capacity  limited, 
and  Ins  measure  of  performance  was  taken  to  be  the  average  delay  encountered  by  a message  in  passing  through 
die  net.  lie  addressed  <|ucst  jiois  pet  laming  m the  assigmnem  of  channel  capacities,  effect  of  priority  discipline, 
Jioice  of  routing  p.ocedmc.  ami  deMgn  n|  iopolog,Lal  M.ueture,  Allliougli  Klcinmck's  initial  efforts  preceded 
|>i esotit  day  computet  networks,  his  concepts  are  the  basis  lor  most  analytic  studies  in  this  area.  The  system  to 
which  he  directed  Ins  teelmujues  were  associated  vsiih  the  automatic  telegraph  switching  system  of  Western 
Union  |o-  Air  l;urcc  military  iialiic  | VI.K5.S |, 

Kleimock  used  the  discipline  ol  the  early  investigators  of  i|ueueing  processes,  Johannscn.1  l-rlang  [BR04K] 
and  numerous  other  contributors  to  the  early  development  of  tjueueing  theory  (Molina  |M022,  27]  and  O’Dell 
|0  1)20,27) ) and  incorporated  these  disciplines  along  with  modern  day  researchers  (sec  Klcinrock  [KL64]  for 
fmihe.  listing)  to  stmlN  communication  nets.  Iljs  efforts  are  notable  in  estahlisliing  some  of  the  mathematical 
tationale  used  in  evaluating  communication  networks. 


Tlk'  w,„k  of  Hu-  Kami  Corporation  in  I,  series  of  memotanda  tilled  “On  Distributed  Communications  to 
,i  I.iiw-I)iiI;i-K:i(c  Military  Cimiiiiiind  und  Control  System"  w;is  primarily  the  work  of  Itartm  IBA64A-64II] . 
Oil, ei  contributors  in  the  Kami  study  were  HiK-l.m  1 11064]  and  Smith  |SM04|  and  the  Initiation  of  store- 
and-loiwaid  messag.e  switched  digital  lielworks.  'Ibis  series  of  memoranda  brought  out  Karan's  definition 
of  a “packet"  aiong  with  his  “hot  potato"  routing  algorithm. 

I lie  Kami  series  was  primarily  interested  jp  the  survival  of  communications  after  an  attack  on  the  sys- 
tern.  Hie  series  crnicei  tied  itself  with  stochastic  adaptive  muling  techniques  that  used  informal  ion  from  mes- 
sages i Kissing  thmugh  the  network  to  adjust  muting  tables;  deleiminislic  techniques,  which  use  dynamic 


kck'K'iue  ti>  Jnli.miiH'ii’s  work  eiin  hi*  lomul  in  [Urmkmeyer  4 8 j . 
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programming  or  graph-theoretic  algorithms  to  recalculate  changes  in  the  tables  from  observed  changes  in  the 
network. 

The  groundwork  set  down  in  the  Rand  series  led  to  the  development  of  computer  communication  net- 
works, the  most  ambitious  effort  undertaken  iO  da^e  being  the  Advanced  Research  Projects  Agency  (ARPA) 
Computer  Network,  or  more  simply  the  ARPA  NET.  The  evolution  ol  this  network  is  rather  extensive  ar.u 
a very  good  bibliography  on  the  ARPA  NET  is  given  by  Wolf  [WO 73  J,  who  also  gives  an  excellent  overview  of 
the  ARPA  network. 

The  prime  mover  behind  the  creation,  growth,  and  development  of  the  ARPA  computer  system  is  cred- 
ited to  Roberts  [R067].  He  indicated  that  existing  computer  systems  could  be  more  effectively  utilized  by 
interconnecting  them  in  a network  so  that  their  resources  could  be  shared.  The  idea  of  resource  sharing  was 
not  entirely  a new  concept,  but  ins  .*ad  a natural  outgrowth  of  earlier  interconnected  systems,  such  as  the 
Semi-Automatic  Ground  Environment  ^SAGE)  air  defense  system  of  the  fifties  (Everett  57]. 

There  were  numerous  computer-connected  systems  that  existed  prior  to  ARPA  and  are  historically  re- 
viewed by  Fultz  [FU72J . One  of  the  early  systems  was  the:  American  Airlines  SABRE  Reservation  System 
[PL61 1,  that  led  to  the  present  systems  used  by  other  airlines,  hotels,  etc.  The  need  for  improved  military" 
data  communications  led  to  the  Automatic  Digital  Network  (AUTOD'N)  Communications  System  in  1%3 
[Millar  68].  This  syr,cm  depended  on  both  line  and  message  switching  facilities  with  added  constraints  of  sur- 
vivability and  vulnerability. 

In  addition  to  computer  systems  listed  by  Fultz,  there  is  a tutorial  wri  -‘-up  of  four  currently  operating 
computer-communication  networks  [Schwartz  72].  They  include  the  TYMNET  network,  the  G.E.  Information 
Services  Network,  the  NASDAQ  over-the-counter  stock-quotation  system,  and  the  Computer  Sciences  Infonct. 
These  systems  all  contain  some  mixture  of  computers  and  networks  and  possess  to  various  degrees  the  need  for 
store-y*  ' forward  capability.  This  store-and-forward  technique  of  sending  and  retrieving  data  as  opposed  to 
dedicated  connections  while  hanu.mg  data  is  the  uniquer.es  f Cv  mputer-communications.  It  should  be  noted 
that  in  the  TYMNET  network,  they  prefer  to  describe  the  function  of  their  network  as  virtual  line  switching 
as  opposed  to  message  switching. 

Computer  networks,  while  solving  possible  cost  effectiveness  of  computer  systems  by  bringing  massive 
hardware  and  software  resources  to  users  with  moderate  investments,  also  brought  massive  problems  of  network 
structure,  Handling  and  transmission  of  message  routing  procedures,  and  CPU  protocols.  This  then  is  the  chal- 
lenge of  the  seventies,  to  extend  to  more  users  the  vast  computer  resources  on  a global  basis,  with  minimum 
cost,  high  reliability,  remote  accessing,  maximum  throughput,  all  with  a minimum  of  effort  and  time  for  the 


user. 
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1 .2  PROBLEM  CONSIDERED 

The  intent  of  this  research  effort  is  to  establish  a tractable  adaptive  routing  procedure  applicable  to 
store-and-forward  computer-communication  networks.  The  routing  assignment  will  be  on  a node-to-node 
basis  and  will  essentially  be  an  extension  of  Baran’s  [BOEHM  and  BARAN  64]  “hot  potato”  routing  tech- 
nique with  priority  assignments.  The  concepts  of  adaptive  routing  are  not  new,  but  are  lacking  in  true  adapt- 
ability due  to  network  topology  changes  (new  links  or  blocked  links),  or  the  “deadly  embrace”  of  “looping” 
as  well  as  “ping-ponging  It  is  unlikely  that  there  can  be  an  all  encompassing  optimal  routing  technique  for 
all  network  designs  and  all  system  constraints.  The  definition  of  optimal  is  open  to  debate,  since  it  could  be 
conditioned  on  cost,  time-delay,  guaranteed  delivery,  minimum  path,  maximum  throughput,  etc. 

The  definition  of  optimality  that  will  be  adopted  here  will  be  in  the  context  given  under  dynamic  pro- 
gramming principles  vs  applied  to  optimum  paths  and  states: 

The  optimal  path  from  X to  A has  the  property  that,  for  any  node  Y along  the  path,  the  remainder  of 
the  path  is  the  optimal  path  from  Y to/».  If  we  define  N(X)  as  the  set  of  neighbors  of  X connected  to  A'  by 
links,  we  have 


D(X)  = Min  [d{X,  T)  + DfT)].  (1.1) 

YeN(X) 

The  procedure  will  be  to  simulate  some  routing  strategy  for  a given  network  topology  or  general  class  of 
nr* works  and  measure  the  effective  adaptability  of  the  routing,  when  it  is  compared  against  various  perfor- 
mance measures.  Performance  will  be  measure  1 on  the  basis  of  average  message  delay,  throughput,  and  per- 
cent of  undelivered  messages.  When  dealing  with  an  operational  network,  a measurement  of  average  message 
delay  includes  such  items  as:  queueing  delays,  line  and  node  blockages,  node  processing  delays,  and  message 
retransmission  due  to  channel  errors. 

The  measure  of  performance  will  generally  be  obtained  from  computer  simulation,  since  mathematical 
models  do  not  lend  themselves  to  analysis  of  network,  performance.  This  absence  of  a matiiematical  model 
to  effectively  evaluate  a network  in  closed  form,  means  one  has  to  represent  the  system  paiametricaily,  simu- 
late the  system,  analyze  the  results,  and  ascertain  what  changes  should  be  made. 

The  major  emphasis  in  this  research  will  be  to  demonstrate  how  some  relatively  simple  add-ons  to  already 
existing  adaptive  routing  techniques  can  decrease  the  average  message  delay  and  correspondingly  increase  mes- 
sage throughput  in  the  network.  A further  objective  of  this  study  will  be  to  assess  the  effects  of  priori  v assign- 
ment to  messages  that  have  reached  some  specified  aging  threshold  while  in  the  network  and  note  the  effects 
of  such  priority  assignment  on  network  performance.  The  technique  will  involve  working  on  a node-to-node 


“Looping"  is  when  a message  repeatedly  travels  the  same  set  of  nodes.  “Ping-punging"  is  where  a message  goes  back 
and  forth  between  the  same  two  nodes. 
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basis  within  some  specific  network  topology  and  optimize  the  routing  while  attempting  to  minimize  the  hard- 
ware and  software  requirements  at  the  nodes. 

1 i OUTLINE  OF  THE  DISSERTATION 

Chapter  II  is  an  overview  of  routing  strategies  and  routing  classifications.  Chapter  III  treats  piiority 
assignment  strategies  and  optimization  of  a defined  measurement  parameter,  throughput  factor,  0(t,k)  and 
how  it  relates  to  the  priority  threshold  setting  k.  A closed  foan  solution  for  0(t,k)  is  obtained  on  a 3-node 
network  assuming  infinite  buffers  and  fixed  routing,  The  problem  of  analytically  modeling  a store-and-for  ward 
system  with  adaptive  routing  is  also  addressed  in  Chapter  III. 

Chapter  IV  and  V are  simulation  studies  on  an  8*node  and  19-node  network  respectively.  The  simula- 
tions are  a means  of  verifying  the  performance  of  some  selected  algorithms  with  the  modifications  suggested 
in  this  study.  The  adaptability  of  the  routing  algorithms  were  tested  on  undamaged  and  damaged  networks 
of  the  class  mentioned  above.  Chapter  VI  is  a comparison  of  the  numerical  throughput  factor  tc  the  simulated 
one  for  the  19-node  AREA  network. 

1 his  work  is  summarized  in  Chapter  VII  along  with  final  conclusions  and  recommendations  for  future 

work. 
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II.  ROUTING  STRATEGIES 

An  excellent  overview  of  the  basic  routing  problem,  routing  requirements,  and  routing  classification  is 
given  by  Fultz  ( FU72]  and  a biief  summary  of  his  efforts  are  worth  reiterating  here. 

The  basics  of  the  deterministic  routing  problem  have  been  formulated  mathematically  and  some  general 
techniques  exist  for  solutions  to  some  specific  classes  of  problems  which  have  been  cited  by  Ford  and  Fulker- 
son [FG*2].  Their  solutions  are  associated  with  maximum  flow  problems,  the  Hitchcock  transportation,  pro- 
ject planning,  and  the  shot  test  route  problem. 

The  standard  switched  circuit  communication  network  has  a routing  problem  associated  with  finding 
free  communication  links  that  maintain  a fixed  path  for  the  duration  of  the  message  flow.  On  the  other  hand, 
in  a message  switching  computer  communion ‘ion  network,  the  messages  must  be  routed  under  a different 
strategy.  Message  routing  involves  queueing  delays  due  to  he  store  and-forward  aspects  of  louting  strategies. 
This  class  of  nets  utilizes  message  delay  estimates  as  a means  of  measuring  network  perfoimance;  hence,  rout- 
ing strategies  are  of  prime  concern. 

General  requirements  for  message  routing  can  be  listed  in  three  categories.  One,  the  routing  technique 
should  adapt  to  changes  in  network  topology  resulting  from  node  or  line  failures.  Second,  routing  strategy 
should  adapt  to  varying  source-destination  traffic  loads  to  insure  minimal  message  delays.  Third,  contingency 
plans  should  exist  for  dropping  of  undeliverablc  messages;  i.e.,  messages  should  be  dropped  after  some  pre-set 
time  lapse. 

Store-and-forward  routing  is  in  the  same  context  as  packet-switched  networks.  In  cither  vernacular,  the 
basic  unit  of  information  exchange  is  the  “packet”  and  varies  from  1008  to  1024  bits  depending  on  whether 
the  header  is  included.  The  ARPA  system  accepts  messages  up  to  8095  bits  in  length  and  divides  these  into 
1008  bit  packets. 

2.1  ROUTING  CLASSIFICATION 

Before  analyzing  some  of  the  popular  routing  algorithms,  particularly  those  pertaining  to  adaptive  tech- 
niques, it  would  be  instructive  to  list  the  routing  classification. 

It  was  stated  previously  that  there  were  three  major  categories  foi  classifying  routing  techniques  and 
they  were  given  by  Fultz,  Boehm  and  Mobley  [B066],  along  with  a quick  summary.  Fultz’s  list  augmented 
with  the  RAND  series  are  as  follows: 

A.  Deterministic  Routing 
1.  Flooding 

a.  All  links 

b.  Selective  links 
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2.  Fixed  routing 

a.  Table  look-up 

3.  Split  traffic 

a.  Optimum  bifurcated 

b.  Suboptimum  bifurcated 

4.  ideal  observer  ( scheduling  problem) 

a.  Present 

b.  Future 

3.  Ideal  routing  table  and  shortest  path 

6.  Dynamic  programming  (solution  of  the  shortest  path) 

7.  Butrimenko’s  method  (of  updating  shortest  paths) 

B.  Stochastic  Routing 

1.  Random  muting 

2.  Lolated 

a.  Local  delay  estimates 

b.  Shortest  queue  + bias 

3.  Distributed 

a.  Periodic  updating  (nearest  neighbor) 

b.  Asynchronous  updating  (Percolation) 

4.  Baran's  “hot  potato” 

5.  B.  wards  leai nmg 

6.  Negative  reinforcement 

7.  Bi-udaptive 

X.  Superposition 
(\  Mow  Control  Routing 

1.  Isarithmic  network 

2.  Butter  stoiage  allocation  schemes 

3.  Special  route  assignment  algorithms 

Any  attempt  to  go  into  any  significant  detail  of  all  the  possible  routing  schemes  in  the  various  categories 
would  be  too  much  of  an  undertaking  and  would  detract  from  the  issu  at  hand  and  can  be  found  in  the  cited 
references.  The  listing  is  to  he  p keep  techniques  in  their  proper  perspective  and  a given  routing  scheme  will  be 
elaborated  on  when  deemed  appropriate  for  a given  discussion. 

A brief  resume  of  some  ot  the  routing  techniques  will  be  given  next;  details  can  be  found  in  Fultz  [FU72J, 
Boehm  and  Mobley  |B()6f>.  (>‘>J.  oi  Fultz  and  Kleintock  |FU71  J. 
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2.2  DETERMINISTIC  ROUTING 

Deterministic  routing  is  where  the  routes  are  selected  from  the  network  structure,  based  upon  the  mini- 
mum time  delay  to  reach  a destination  under  ideal,  loop-free  conditions. 

2.2.1  Flooding 

Each  node  receiving  or  originating  a message  transmits  i«  over  every  outgoing  link,  after  checking  to  see 
that  it  has  not  previously  transmitted  this  same  message,  or  that  it  is  noi  the  destination  of  the  message.  In- 
efficiency is  the  price  paid  for  this  technique. 

2.2.2  Selective  Flooding 

Each  station  will  transmit  only  over  links  heading  in  the  general  direction  of  the  destination.  Storage  of 
tables  at  each  node  for  the  best  outgoing  link  is  the  price  paid  by  this  method.  In  addition,  if  nodes  and  links 
arc  disabled,  some  sources  and  destinations  would  be  connected  only  by  round-about  paths  excluded  from 
consideration  by  the  selective  Hooding  technique. 

2.2.3  Fixed  Rooting 

Fixed  routing  specifies  a unique  path  over  a selected  set  of  links  n(Nr  ...Ay)  for  a given  source  node  A', 
and  destination  node  V This  is  accomplished  by  a table  look-up  from  a routing  table  contained  at  each  node. 
Fixed  routing  reouircs  completely  reliable  nodes  and  links. 

2.2.4  Sp!  f Traflic 

This  techi.ique  has  been  reported  by  Fultz  IFU72J,  and  as  opposed  to  fixed  routing,  allows  traffic  to 
flow'  over  more  than  one  set  of  paths  between  a given  source-destination  node  pair.  This  splitting  of  traffic  is 
called  traffic  bifurcation.  The  suboptimal  bifurcation  is  used  since  it’s  more  simple  in  structure  and  requires 
less  computation  time  than  the  optimal. 

2.2.5  Ideal  Observer 

This  technique  is  essentially  a scheduling  problem.  Each  new  packet  entering  the  system  has  its  mini- 
mum route  computed  for  the  given  source-destination  pair.  This  newly  computed  route  is  based  upon  the 
complete  present  information  about  the  packets  already  in  the  network  and  their  known  routes.  If  vhc  ideal 
observer  has  information  about  future  events,  this  is  also  utilized  in  computing  the  route.  This  technique  is 
more  theoretical  than  practical,  but  serves  as  a basis  for  comparing  other  routing  techniques. 

2.2.6  Ideal  Routing  Table  and  Shortest  Path 

Let  X,  Y,  and  A denote  arbitrary  nodes  in  a network,  and  let  <J(X,Y)  denote  transit  time  along  the  link. 
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If  node  A is  the  destination,  the  ideal  routing  table  entry  gives  the  shortest  distance  from  node  A to 
node  A starting  along  link  AT,  and  may  be  represented  as  a function  R(X.Y)  of  A'  and  Y alone.  If  one  defines 

/)(>')  as  the  distance  along  the  shortest  path  from  Y to  A , it  follows  that 

R(X,  Y)  = J(X,  Y)  + D(  Y),  for  node  A.  (2.1) 

Hence,  if  /)(K)  is  determined  for  all  nodes  then  one  can  compute  the  table  entries  R(X,  Y)  of  the  ideal 

routing  table  iTorn  equation  ( 2.1 ).  The  report  by  Boehm  and  Mobley  [B066]  elaborates  on  many  of  these 
routing  procedures. 

2.2.7  Dynamic  Programming 

The  dynamic  programming  principle  of  optimality  states. 

The  optimal  path  from  A'  to.-l  has  the  property  that,  lor  any  node  Y along  the  path,  the  remain- 
der of  the  path  is  the  optimal  path  Iron.  )'  to  A.  If  we  define  AT  AY  as  the  set  of  neighbors  of  X connected  to 
X b\  links,  then 


/>(.V)  = Mm  |(/(A'.>  ) ♦ IM  >')]  U.2) 

Yt  A'l  X) 

Dynamic  programming  requires  storing  routing  tables  at  each  node  as  well  as  a map  of  the  network.  Any 
clianges  in  the  net  are  sensed  In  neighboring  nodes  that  can  communicate  these  changes  to  all  other  nodes.  Dy- 
namic programming  is  Heated  in  B Hainan  and  Diey  lus  | Bl.l  Ci2). 

2.2.S  Butrimcnko’s  Method 

Butnmenko’s  |MU(>4|  method  o!  updating  the  shortest  path  uses  perturbation  techniques  that  adjust  to 
simultaneous,  or  near  simultaneous,  failures  or  additions  of  nodes  and  links. 

The  Butnmenko  method  checks  lor  possible  link  or  node  failures  in  the  network  by  verifying  at  each 
node  A'f  along  a message’s  route  whether  or  hot  it  is  still  the  minimum  route  to  take. 

2.3  STOCHASTIC  ROUT  I Mi 

Stochastic  routing  operates  on  a probabilistic  set  of  decision  rules.  Routes  are  based  on  the  actual  state 
of  the  network  or  estimates  ot  the  present  state.  Stochastic  routing  as  well  as  deterministic  routing  were  ini- 
tially suggested  by  Boehm  |BOb(>| 
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2.3.1  Random  Routing 

Random  routing  uses  a decision  rule  for  next  node  to  visit  based  on  some  probability  distribution  over 
the  set  of  neighbor  nodes.  The  set  of  nodes  can  be  all  or  some  selective  set  that  are  in  the  general  direction  of 
the  message’s  destination. 

Prosser  (PR52)  and  Kleinrock  (KL64)  investigated  some  of  these  random  routing  techniques.  They  con- 
cluded that  although  such  routing  may  be  relatively  unaffected  by  small  changes  in  the  net  structure,  they  were 
highly  inefficient  in  terms  of  message  delay. 

2.3.2  Isolated 

The  source  of  information  available  for  developing  the  delay  table  estimates  determines  whether  the 
routing  technique  belongs  to  the  isolated  or  distributed  class.  If  the  only  information  available  to  a node  comes 
via  packets  flowing  through  it.  then  the  procedure  is  designated  as  isolated,  or  a local  delay  estimate. 

Shortest  queue  + bus  is  a technique  where  a packet’s  route  is  selected  by  placing  it  in  the  shortest  output 
channel  queue.  This  technique  is  used  in  the  ARP  A system  | Heart  70)  and  is  extensively  dealt  with  in  Fultz’s 
[FU72]  work.  This  pmceduie  will  be  further  developed  in  section  (2.6.2). 

2.3.3  Distributed 

Opeiates  in  the  same  basic  way  as  the  isolated  technique.  In  this  routing  ‘echnique,  neighboring  nodes 
exchange  status  information  as  regards  to  delay  and  routing  information.  This  information  exchange  is  accom- 
plished by  special  information  packets;  hence,  routing  data  percolates  throughout  the  network, 

2.3. ; Raran’s  “hot  potato” 

In  a marner  analogous  to  selective  flooding  theie  is  a technique  whereby  the  table  ai  each  node  has  not 
only  the  best  outgoing  link  for  each  destination,  but  also  the  second  best  through  «th  best,  enumeiating  all  n 
outgoing  links.  Thus  an  incoming  message  can  be  immediately  routed  out  on  the  best  available  link  at  any 
given  instant.  This  selection  of  outgoing  links  is  called  “hot  potato”  routing,  and  was  extensively  investigated 
by  Paul  Raran  in  the  RAND  Memoranda  series  [BA64A-641J. 

2.3. i Backwards  Learning 

In  the  report  by  Boehm  anu  Mobley  IB066),  the  value  of  the  hand-over  number  {UN)  contained  in  a 
message  received  along  the  incoming  link  {If)  estimated  the  shortest  time  currently  necessary  to  go  from  the 
source  S to  the  present  node  X.  The  backwards  learning  t echnique  (Used  by  Baran  |BA64]  in  the  “On  Distrib- 
uted Communications”  series)  utilizes  the  fact  that  HN  also  .stimates  the  shortest  time  T ( S.IJ)  currently  nec- 
essary to  go  from  X to  S,  using  U as  the  outgoing  link  (LO).  The  value  T {S.Lf) OLI)  of  the  routing  table  at  X 
is  then  updated  by  a simple  averaging  process. 


2.3.6  Negative  Reinforcement 

Negative  reinforcement  is  a class  of  techniques  that  penalizes  various  routing  table  entries  by  various 
amounts  whenever  the  delay  coming  into  a node  is  equal  '•>  the  delay  going  out  from  said  node  via  the  same 
link.  It  adds  some  constant  delay  value  to  its  old  routing  table  values  whenever  this  equainy  occurs.  The  in- 
tent of  tli is  operation  is  to  reduce  the  ‘'ping-pong”  effect 

2.3.7  Bi-adaptive 

Bi-adaptive  is  still  a further  synthesis  of  backwards  learning  and  negative  reinforcement,  where  table  en- 
tries can  only  be  increased  by  negative  leinforcement  and  decreased  by  backward  learning. 

2 j.8  Superposition 

Superposition  is  the  linear  combination  of  backwards  learning  and  negative  reinforcem  at. 

2.4  FLOW  CONTROL  ROUTING 

Flow  control  routing  involves  techniques  to  coordinate  the  network  message  routing  algorithm  and  the 
activities  of  die  software  at  source  and  destination  nodes  regarding  sequencing  of  messages,  network  connec- 
tivity, buffer  storage  allocution,  etc.  Flow  control  routing  was  studied  and  classified  by  Fultz  [FU72J. 

2.4.1  Isarithmic  Network 

Davies  (DA71  ] proposed  the  isarithmic  network  as  a solution  to  the  network  congestion  problem.  He 
observed  that  the  level  of  traffic  within  the  network  could  be  expressed  b>  the  number  of  packets  in  transit. 

His  solution  was  to  prevent  congestion  by  placing  a limit  on  die  total  number  of  packets  in  the  network.  Thus, 
he  defined  an  “isarithmic”  network  as  one  where  the  total  number  of  packets  is  held  const,  m.  The  idea  was 
basically  to  have  a fixed  number  of  packets  in  the  system  eitl  er  emptv  or  c irrving  data.  When  data  is  icmoved 
from  a packet  an  empty  one  is  available  io  acquire  new'  data. 

2A2  Buffer  Storage  Allocation  Schemes 

Buffer  storage  ullocat.on  is  the  technique  of  coordinating  network  routing  with  that  of  congestion  con- 
trol algorithms  to  eliminate  bottlenecks  at  the  buffer  queues.  1:1  foe t ive  buffer  storage  schemes  enable  high 
throughput  levels  for  the  net  as  well  as  acceptable  average  message  delays.  Buffer  techniques  are  discussed  by 
Zeigler  [ZL7I J. 

2.4.3  Special  Route  Assignment  Algorithms 

Route  assignment  algorithms  aie  the  attempts  to  avoid  establishing  instantaneous  alternate  routes  in 
response  to  rapid  changes  in  traffic  flow  within  the  network.  They  control  the  rate  at  which  traffic  is  allowed 
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to  increase  on  paths  in  a net  and  the  interval  before  additional  alternate  routes  are  established.  This  subject 
has  been  investigated  by  Fultz  [FU72J,  Kahn  and  Crowther  (KA71 J and  Frank,  et  al  [FR72J. 

2.5  ADAPTIVE  ROUTING 

The  present  day  concept  of  adaptive  routing  can  probably  be  attributed  to  the  RAND  memorandum  by 
Boehm  and  Mobley  (B067J.  Their  objective  was  to  reroute  messages  in  a military  communication  system  where 
tne  lines  of  communication  had  been  damaged  or  destroyed.  They  were  investigating  an  effective  adaptive 
routing  technique  that  would  adjust  the  routing  tables  of  the  message-switching  control  system  to  a changing 
network  topology.  Present  da>  computer  communication  networks  need  to  employ  the  same  adaptive  rout- 
ing strategies,  nut  for  safeguarding  against  attack,  but  for  keeping  data  from  getting  blocked  at  overloaded 
nodes. 

Summarizing  the  previous  sections  we  can  say  that  adaptive  routing  techniques  fail  into  two  basic  classes: 

1 ) Stochastic  techniques  that  estimate  the  time  required  to  reach  a destination  from  observations  of  messages 
passing  through  nodes,  and  2‘  deterministic  techniques  that  compute  the  time  required  to  reach  a destination 
under  ideal  conditions  based  on  network  structure.  Fu.t/  lis’s  a third  technique,  called  flow  control  routing 
that  incorporates  butler  storage  allocation  schemes  and  special  route  assignment  algorithms. 

Routing  requirements  tor  packet  switched  networks  can  be  listed  briefly  as  follows: 

1 . Message  routing  should  insure  rapid  and  error-free  delivery  of  messages. 

2.  The  routing  strategy  should  adapt  to  changes  in  the  network  topology  resulting  from  node 
and  communication  link  failures. 

3.  The  routing  technique  should  udap:  to  varying  source-destination  traffic  loads. 

4.  Packets  should  be  routed  around  nodes  that  are  congested  or  temporarily  blocked  due  to  full 
storage. 

5.  Packets  should  be  routed  to  their  destinations  via  the  fewest  intermediate  nodes. 

6.  A technique  for  detecting  “loops  and  ping-ponging“  should  be  incorporated  in  the  routing 
strategy. 

7.  The  routing  technique  should  be  as  simple  as  possible  to  minimize  hardware  and  software 
requirements  at  each  node. 

2.6  ROUTING  SCENARIOS  COMPARED  IN  THIS  RESEARCH 

The  following  algorithms  arc  briefly  outlined  along  with  their  associated  acronyms.  The  ou  dines  are 
intentionally  brief,  since  there  arc  ample  references  cited  in  this  and  previous  sections. 

Routing  tables  are  initially  calculated  by  computing  the  minimum  path  matrix  on  an  empty  network. 
Floyd’s  algorithm  (FL062)  is  generally  used  to  solve  the  shortest  paths  for  setting  up  the  node  delay  tables. 
(Frank  71  j has  developed  a computationally  more  efficient  routing  matrix  that  the  author  claims  would  yield 
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a near  optimal  r juting  strategy.*  These  initializing  algorithms  would  be  essentially  optimum  if  the  network 
either  did  not  become  too  congested  and  there  were  no  changes  in  net  topology  due  to  link  or  node  failures. 

In  all  cases  of  i outing  the  following  specifications  should  serve  as  a guide. 

1.  Route  selection  should  be  independently  performed  at  each  node,  based  on  nearest  neighbor 
information  and  traffic  patter'  in  the  immediate  area. 

2.  Individual  routing  should  be  globally  sensible. 

3.  Updating  should  be  conditioned  on  queue  size  at  the  next  best  nodes;  i.e.,  update  if  traffic 
(messages)  start  to  overflow  at  a given  queue. 

4.  Paths  with  fewest  nodes  to  destination  should  be  considered  first. 

2.6.1  ARPA  Routing  Strategy 

The  ARPA  network  utilizes  Inter  Message  Processor  (IMP’s)  that  are  autonomous  to  the  system.  Each 
IMP  maintains  its  own  routing  table  that  it  updates  every  half  second  using  nearest  neighbor  estimates  of  delay 
to  each  destination.  ARPA  routing  is  adaptive,  although  it  is  not  optimal.  There  is  also  an  instability  in  the 
ARPA  system  caused  by  such  things  as  node  A routes  messages  to  node  C via  node  B and  B beheves  the  best 
path  to  C is  via  node  A . This  instability  in  routing  is  generally  called  ping-pong  effect  and  in  th.  example  would 
be  ping-ponging  between  nodes  A and/i. 

The  routing  algorithm  directs  each  packet  to  its  destination  along  a path  for  which  the  total  estimated 
transit  time  is  minimum  [Heart  70 j.  This  path  is  noi  determined  in  advance;  instead,  each  IMP  individually 
decides  which  of  its  output  lines  should  transmit  a packet  addressed  to  another  destination.  The  selection  is 
made  by  a relatively  fast  and  simple  lookup  procedure,  wlieie.  for  each  possible  destination,  an  entry  in  tuc 
table  designates  the  appropriate  next  leg.  These  entries  reflect  line  oi  IMP  trouble,  traflic  congestion,  and  cur- 
rent subnet  connectivity.  This  routing  tabic  is  updated  evei\  half  second  as  follows: 

Each  IMP  estimates  the  delay  it  expects  a packet  to  encounter  in  rcachinii  every  possible  destination  ovei 
all  of  its  output  lines.  The  minimum  delay  estimate  is  selected  foi  each  destination  and  periodically  (about 
twice  a second)  passes  these  estimates  to  its  immediate  neighbors.  Each  IMP  then  constructs  its  own  routing 
table  by  combining  its  neighbors'  estimates  with  its  own  estimates  of  delay  to  that  neighbor.  The  estimated 
delay  to  each  neighbor  is  based  on  both  queue  lengths  and  recent  performance  of  the  connecting  communica- 
tion network.  For  each  destination,  the  table  is  then  made  to  speeily  that  selected  output  line  for  which  the 
sum  of  the  estimated  delay  to  the  neighbor  plus  the  neighbor’s  smallest  delay,  is  a minimum  for  a given 
destination. 

The  routing  table  is  consistently  and  dynamically  updated  to  adjust  for  changing  conditions  in  the  net- 
work. This  essentially  makes  the  system  adaptive  to  fluctuations  of  line  traffic.  IMPS,  and  congestion.  It  is 

* Optima!  routing  as  used  in  tins  study  implies  routing  that  mminii/cs  flu-  average  message  delay. 
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not  necessary  lor  an  IMP  to  know  the  topology  ol  the  network.  Indeed,  an  IMP  need  not  know  the  identity 
ot  its  immediate  neighbors.  This  enables  the  leased  circuits  to  be  reconfigured  to  2 new  topoiogy  without 
requiring  any  changes  to  the  IMP’s. 

In  the  ARPA  network  the  basic  unit  of  information  passed  between  any  pair  of  nodes  is  called  a “packet” 
with  maximum  size  ot  approximately  1000  bits.  Messages  that  originate  at  a HOST  computer  (any  number  of 
large  computers  serving  the  ARPA  net  are  called  HOST)  have  a maximum  length  of  approximately  8000  bits. 
The  IMP  breaks  these  large  message  packets  into  1000-bit  segments  or  packets.  These  packets  are  then  handled 
by  the  network  as  independent  entities  mtii  they  reach  their  destination  node.  Packets  of  a multi-packet  mes- 
sage have  headers  associated  wylt  them  lor  error  cnccking,  how  many  packets  in  a message  and  their  ordering, 
and  what  their  destination  node  <s.  When  packets  reach  their  destination  node,  they  are  reassembled  before 
transt erring  to  the  destination  HOST.  Single  packet  messages  are  given  higher  priority  than  multipacket  mes- 
sages so  that  the  network  can  support  interactive  users. 

2.6.2  Shortest  Queue  + Bias  *■  Periodic  Updating  (SQ  + B + PUD) 

the  ARPa  syjtem  uses  an  algorithm  for  routing  called  shortest  queue  + bias  + periodic  updating 
(S(J  + B + PUD).  This  scheme  uses  minimum  delay  vectors  to  propagate  from  node  to  node  any  changes  in 
the  system.  The  routing  strategy  uses  a combination  of  the  present  table  entries  and  queue  lengths  at  the  node 
to  choose  the  “best”  outgoing  link,  and  its  scenario  is  as  follows: 

For  each  point  in  time,  calculate  minimum  time  delay  (vector)  to  go  from  say  node  / to  node  J:  i.e., 

Min  (/./)( Fig.  2.1). 

Node  / will  ask  ail  of  its  neighboring  nodes Kj  what  is  their  min.  delay  vector;  i.e.,  (A.,/). 


All  Kj ’s  are  connected  directly  to 
/,  but  not  necessarily  to  J. 


Take  minimum  of  all  neighboring  nodes  Ki  and  this  establishes  the  following: 


(2.3) 


Min (/,/)  = Min] (A\,/)  + M/.A.)  + £(/*.) | . 
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where 

( KJ)  = path  length  from  node  A}.  to  node  J, 

9.(1,  K. ) = link  time  from  node  / to  node  A^. , 

£?1/,A'.)  = queue  length  at  node  / on  output  link  to  node  A'. , 
and 

A.  = a neighboring  node  of/. 

New  routing  table  now  uses  the  new  delay  values 

<H=|  (A-,../) + «(/.*:,.)}  . (2.4) 

Hence,  the  new  routing  table  becomes: 

(2.5) 

The  next  node  A is  chosen  such  that 

QU, A\,)  + = Min  j Q U,K. ) + A/.  ( , (2.6) 

i 

Hence 

Min  (/,/)  = .V(/  + B + /m  (2.7) 

2.6.3  Preas,:gned  Links  (PAL) 

Preassigned  links  (PAL)  is  a technique  where  outgoing  links  for  incoming  messages  (either  external  or 
internal)  are  assigned  as  soon  as  messages  get  on  the  queue;  i.c.,  buffer  storage.  Outgoing  links  are  chosen 
from  the  routing  tables;  hence,  for  no  updating,  link  assignment  if  fixed  routing.  When  updating  is  em- 
ployed, link  assignments  change  according  to  what  constitutes  the  new  minimal  path. 
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It  would  take  considerably  more  programming  to  check  messages  in  queue  and  to  recalculate 
optimum  path,  than  it  takes  to  check  available  outgoing  links  against  those  messages  that  have  that  link 
assignment. 

2.6.4  Backwards  Learning  (BL) 

In  the  report  by  Boehm  and  Mobley  [B066] , the  value  of  the  hand-over  number  (HN)  contained  in 
a message  received  along  the  incoming  link  (LI)  estimated  the  shortest  time  currently  necessary  to  go  from 
the  source  5 to  the  present  node  X.  The  bacbi'ards  learning  technique  fused  by  Baran  in  the  “On  Distrib- 
uted Communications”  series)  utilizes  the  fact  that  HN  also  estimates  the  shortest  time  71(5,  LI)  currently 
necessary  to  go  from  X to  5,  using  LI  as  th?  outgoing  link  (LO).  The  value  7(5,  Z./)^  D of  the  routing 
table  at  X is  then  updated  by  a simple  a\eraging  process: 

ns.  = ns.U)oiri  + (K(//A'-r0)]  • |ttV-r0],  (2.8) 

where 

!k  j , if  I IN  Tn  < 0 =*■  learning  constant 
k2,  if  HN-Tq  > 0 =*  forgetting  constant, 


and 


0 < kvk2  < 1. 

*he  new  estimate  becomes  a certain  fraction  A'(0  < K < 1)  of  the  value  between  the  old  and 
the  current  HN.  This  procedure  includes  a “learning  constant”  k j , used  to  decrease  T(S,  Lf)()y  and  2 
“forgetting  constant  k y used  to  increase  it.  If  more  sophisticated  statistical  estimations  are  employed, 
these  would  require  extra  information  on  each  table  entry,  thus  increasing  the  computer  cycle  requirement 
at  each  node  and  defeating  the  simplicity  strived  for  in  this  study. 

Backwards  learning  (BL)  is  then  defined  as: 

Tu  - Tn  ♦ A [HN-  Tq ] , (2  9) 


where 


A ■ »/i 
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The  major  disadvantages  of  the  backwards  learning  technique  are  that  it  utilizes  no  direct  informa- 
tion on  the  progress  of  messages  through  the  system,  and  the  indirect  information  supplied  by  the  above- 
mentioned  method  is  often  inadequate  or  misleading.  BL  updates  towaid  S and  not  toward  Z>,  and  in 
cases  where  there  is  node  failure  either  due  to  hardware  or  software  failure,  BL  will  not  adjust  to  these 
new  constraints  and  looping  will  occur  or  “ping-ponging”  between  nodes  until  the  message  has  reached 
its  maximum  time  to  be  in  the  system  and  is  subsequently  discarded.  “Ping-ponging”  occurs  whenever 
LI  = LO. 


NOTE:  HN  is  the  actual  time  of  message  transit  from  point  S to  X and  T0  is  value  in  table  at 
to  go  from  X to  S. 

2.6.5  Negative  Reinforcement  (NR) 

In  the  following  topology,  if  in  going  from  node  S to  node  D via  /,  and  message  is  then  sent  back 
from  J to  / (Fig.  2 2),  the  new  table  vabe  at  J becomes* 

TN  (J,  I , D)  = T0  (J,  I,  D)  + 2 8(4  /)  (2.10) 


and  at  node  / 


rN  (A  J,D)  = TQ  (/,  J,  D)  ♦ 2 8(4/7.  (2.11) 

Thus  NR  penalizes  various  routing  table  elements  whenever  LI  = LO. 

The  disadvantage  of  negative  reinforcement  is  that  it  has  no  way  of  adjusting  to  improvements,  such 
as  replacement  of  a damaged  link  or  the  creation  of  a new  link.  Routing  table  entries  could  have  little  or 
no  relation  to  the  actual  transmission  times  between  nodes. 


J 


Fig.  2.2  - 6 Node  Network 
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2.6.6  Superposition  (SP) 

Superposition  =*  Backward  Learning  and  Negative  Reinforcement.  This  would  be  equation  (2.9) 
used  in  combination  with  equations  (?.  1 0)  and  (2.11),  where  the  latter  two  equations  are  used  to  offset 
“ping-ponging.” 


2.6.7  Bi  Adaptive  (BA) 

Bi  Adaptive  =*  Backward  Learning  and  Negative  Reinforcement,  such  that  in  equation  (2.8) 

K t - 1/2,  1/3,  1/4,  etc.  (some  assigned  value)  and  =0. 

2.6.8  No  Updating  (NU) 

On  a lightly  loaded  net  the  ideal  routing  tables  would  be  essentially  optimum  routing.  The  no 
updating  routine  (NU)  was  used  to  evaluate  the  performance  of  other  routing  algorithms. 

2.6.9  Dynamic  Programming  (DY) 

This  programming  technique  is  analogous  to  th.  “ideal”  observer  and  assumes  that  its  knowledge  of 
what  links  and/or  nodes  are  inoperable  does  not  change.  With  this  quasi  all  encompassing  overview  of  the 
network,  this  *vpe  of  programming  should  yield  the  optimum  path  from  any  source  S to  any  destination 
D via  any  intermediate  nodes N.  with  a minimum  time  delay. 

2.6.10  Last  M Nodes  Visited  (LMNV) 

LMNV  is  t technique  in  the  program  routing  strategy  to  prevent  the  message  from  returning  to  the 
last  node  it  has  visited.  This  -cheme  can  be  extended  to  include  the  last  M nodes  visited  (A/=  1,2,.. .,  A") 
and  will  be  designated  as  (LMNV).  Preventing  messages  from  taking  routes  previously  visited,  is  intended 
to  reduce  message  transit  time  by  cutti  ig  out  the  “looping”  or  “ping-pong”  effects. 

2.6. 1 1 Best  Three  Links  (BTL) 

There  is  a tag  associated  with  each*  message  such  that  when  it  is  prevented  from  taking  the  last  node 
visited,  it  will  be  given  an  option  to  take  one  of  three  best  links  from  a given  node  towards  its 
destination  node  D.  The  use  of  last  M nodes  visited  (LMNV)  and  best  three  links  (BTL)  are  generally 
used  together  iri  the  progrr  nming.  If  the  message  becomes  trapped  in  the  network  where  it  can’t  proceed 
without  backtracking  over  previous  routes,  it  may  be  “killed”  or  discarded  from  ‘he  network  arid  tallied 
as  an  undelivered  message.  The  point  at  which  it  is  discarded  from  the  net  is  deiermined  by  a clocking 
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mechanism*  that  keeps  a tally  of  how  long  a given  message  lias  been  in  the  system  and  when  this  tally 
reaches  some  prescribed  A/DFm  ax  and  the  message  has  not  been  delivered,  then  it  is  discarded  from  the 
system. 

2.6.12  Priority  Assignment  Based  on  Threshold  (km) 

An  aging  parameter  has  been  incorporated  in  the  program  whereby  after  a message  has  reached  some 
percentage  of  A/DjTmax  it  is  given  a priority  while  in  the  queue  so  that  it  can  be  dispatched  through  the 
network  on  a priority  basis.  Once  a message  has  been  upgraded  to  a priority  level  it  continues  to  have 
this  priority  at  all  succeeding  nodes.  The  priority  assignment  is  non-preemptive,  in  that  a message  being 
serviced  at  a given  node  will  continue  to  be  serviced  until  completion,  regardless  of  its  priority  status. 

Note  also  that  messages  are  not  generated  with  a priority,  they  are  only  upgraded  to  a priority  after  having 
spent  some  time  X%  of  A/DT^  AX  in  the  system  (Z)T  > km). 

The  idea  behind  the  upgrading  or  priority  assignment  is  to  enhance  the  delivery  of  messages  already 
in  the  system.  Messages  already  in  the  system  represent  some  investment  as  regards  to  time,  servicing,  and 
use  of  system  capacity.  Once  this  investment  has  been  made  it  would  be  desirable  to  expedite  the  handling 
of  such  messages  when  they  become  stalled  in  the  network,  due  to  looping,  congestion  at  various  nodes, 
etc.  The  level  of  setting  the  priority  threshold  km  will  be  accomplished  by  simulation  and  ascertaining 
where  message  throughput  peaks  when  network  loading  p increases,  i.e.,  what  level  of  thresholding  is 
necessary  to  peak  the  throughput  and  minimize  the  time  delay. 

2.6.13  MODl 

The  routines  associated  wiih  last  M nodes  visited  (LNMV)  and  best  three  links  (BTL)  will  be  desig- 
nated as  modification  sub  one  (MOD^). 

2.6.  i 4 MOD\  + 507? 

When  the  aging  parameter  that  sets  up  priorities  is  included  with  any  programming  routine,  it  will  be 
added  on  as  some  percentage;  thus  MOD^  + 50%  indicated  that  an  aging  parameter  of  k = 509?  was  used. 

It  should  be  noted  here  that  the  techniques  of  preassigned  links  (PAL),  MOD{ , and  MOD^  + 5095- 
are  strategies  being  offered  in  this  study  to  improve  existing  routing  algorithms.  Further,  these  relatively 
simple  add-ons  are  essentially  asynchronous  and  serve  as  congestion  control  and  as  devices  for  “unblocking” 
heavily  loaded  networks.  The  effectiveness  of  the  above  mentioned  seher?s,  as  regards  to  minimizing  the 
average  message  delay  while  increasing  total  message  throughput,  will  be  demonstrated  in  the  simulation 
runs  of  Chapters  IV  and  V. 

* 


Clocking  is  not  synchronous,  but  is  a numerical  incrementation  given  to  each  packet  in  the  system  during  every 
operational  cycle  of  the  simulation  joutinc. 
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III.  PRIORITY  ASSIGNMENT  STRATEGIES 

The  analysis  of  queueing  theory  is  essentially  derived  from  probability  theory  and  as  such,  queueing 
theory  concerns  itself  with  congestion  due  to  random  flows.  The  parameters  of  concern  in  a queueing 
system  differ  according  to  how  the  queues  in  a system  are  structured  and  the  mode  of  operation  of  such 
systems.  The  general  lack  of  knowledge  of  the  system  specifications  renders  extensive  analysis  to  be 
almost  intractable  except  for  a few  idealized  systems.  The  limitations  on  the  analysis  is  directly  due  to 
limitations  of  the  mathematics  itself.  The  system  can’t  be  completely  specified  due  to  no  systematic 
analytic  modeling  and  the  analytical  modeling  has  to  be  too  restrictive  if  it  is  to  be  solved.  The  conse- 
quence of  this  dilemma  results  with  the  mathematical  modeling  of  physically  unrealizable  systems. 

This  chapter  attempts  to  bring  into  focus  some  of  the  standard  analysis  applicable  to  a two-class 
priority  network.  The  analysis  addresses  the  procedures  for  establishing  the  priority  threshold  and  its 
subsequent  effect  on  the  network  as  regard  to  message  delay  and  message  throughput. 

A review  of  queueing  theory  along  with  some  of  the  standard  queueing  formulations  are  given  in 
Appendix  A. 

3 1 PRIORITY  QUEUES 

The  subject  of  priorities  needs  to  be  furil.er  investigated  due  to  the  special  effects  of  assigning  such 
priorities  and  the  unusual  impact  such  assignments  have  on  the  analysis  problem.  Priority  queueing  assumes 
that  the  class  of  messages  now  arriving  at  a given  service  facility,  need  special  handling  either  within  the 
queue  or  with  their  service  requirements.  The  priority  classification  can  be  established  with  some  parameter 
p (1, 2, . . . , k),  where  1 denotes  the  class  with  highest  priority  and  k the  lowest.  The  discipline  by  which 
a server  selects  the  next  message  and  serves  it  is  termed  the  “priority  discipline.”  The  decision  rules  where- 
by any  priority  scheme  must  act  on  are  generally  covered  by  the  following  rules: 

A.  Whether  the  server  continues  or  discontinues  the  servicing  of  a message  already  being  serviced. 

B.  What  message  should  be  next  selected  once  the  server  is  free  to  take  one. 

The  service  selection  may  be  made  exogenously,  i.e.,  dependent  only  upon  knowledge  of  priority 
class  to  which  a message  belongs.  Alternatively,  the  service  selection  may  be  endogenously,  i.e.,  the  deci- 
sion could  be  based  on  the  existing  state  of  the  system,  such  as  type  of  message  last  serviced  or  the 
waiting  time  of  the  messages.  Endogenous  priority  disciplines  are  the  least  known  disciplines  and  are 
seldom  treated  in  the  literature;  it  is  this  class  that  we  attempt  to  analyze  later  in  this  Chapter. 
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With  exogenous  priority  disciplines  the  decision  rule  for  servicing  the  next  message  depends  only  upon  its 
priority  class.  A message  of  the  /th  class  if  present  is  always  taken  for  service  before  a unit  of  the  /th  class 
(/  < /).  In  the  event  a message  of  the  /th  class  is  in  service  and  a message  of  ith  class  {i  < j)  arrives,  further 
alternatives  may  be  invoked. 

A.  Preemptive:  The  service  of  the  / th  class  message  is  immediately  interrupted  and  the  service  is 
given  the  /' th  class  message. 

B.  Non-preemptive:  The  service  of  the  /'th  class  message  is  continued  to  completion. 

In  this  study,  there  will  be  no  preemption;  hence,  no  further  discussions  upon  the  manner  by  which 
the  / th  class  message  is  serviced  on  its  re-entry,  need  be  further  examined.  Service  discipline  for  any  class 
will  be  on  a first-come-first  served  (FCFS)  basis  and  is  the  same  as  first-in-first-out  (FIFO),  where  both 
systems  are  servicing  the  oldest  message  (in  the  queue)  first.  Jiunn  Hsu  [HS7I]  gives  an  overview  of  some 
of  the  standard  service  disciplines.  James  Martin  IMA72J  gives  a very  comprehensive  treatment  of  queueing 
effects  related  to  the  design  and  analysis  of  data  transmission  systems. 

There  will  he  only  two  types  of  systems  addressed  in  this  report.  The  Poisson  input,  exponential 
service,  and  single-server  [M/Mf  1 ) or  the  Poisson  input,  general  service,  and  single-server  (M/G/\ ).  The 
notation  of  M is  for  the  memoryless  Markovian  process  ami  this  is  indicative  to  Poisson  or  exponential 
distributions. 

3.2  WAITING  TIM  I:  DISTRIBUTION  IN  A TWO  CLASS  PRIORITY  SYSTF.M 

In  Appendix  A the  waiting  time  distribution  for  the  steady  state  single  class  queue  was  given  by 
Eq.  (A.  16).  The  problem  here  is  to  obtain  the  waiting-time  distribution  in  a two  class  priority  system. 

We  should  like  to  avoid  attempting  to  solve  what  could  be  a very  difficult  queueing  problem.  The 
intent  here  is  to  use  some  of  the  already  existing  theory  to  get  an  approximation  of  what  to  expect  when 
the  simulation  is  performed.  Miller  |M 1 59]  made  some  very  significant  contributions  in  his  study  of 
"Priority  Queues.”  He  solved  the  waiting-time  distribution  for  a two  class  priority  in  a head-of-thc-line 
priority  queue.  His  analy  sis  dealt  with  the  M/M/ 1 and  M/G/ 1 single  server  system. 

Kesten  and  Runnenburg  j KHS5  / } have  deiived  a generalized  characterization  for  the  Laplace  trans- 
form of  the  steady  state  waiting-time  distribution  of  a K class  priority  system.  Their  efforts  have  been 
duly  referenced  by  most  writers  investigating  priority  queues.  Although  Kesten  and  Runnenburg  devel- 
oped what  appears  to  be  a very  informative  piece  of  theoretical  analysis,  they  apparently  made  an  error 
when  taking  the  inverse  for  the  two  priority  class  case. 

In  1965  a group  of  researchers  at  Stanford  Research  Institute  [NEEoSj  noted,  and  corrected  the 
error  on  the  inversion  formula  of  Kesten  and  Runnenburg.  The  waiting-time  distribution  for  the  priority 
case  as  given  by  Kesten  and  Runnenburg  becomes  (/  > 0): 

( 1-^aQ/ 

P 


Hj  (/)  = 1 ->41  exp 


.P  < 1 


(3,1) 
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where  the  service  distributions  for  both  priority  classes  are  equal  and  given  by 


(*)  = F2(x)  = 1-exp  {-x/p)  . 


(3.2) 


The  corrected  waiting-time  distribution  for  the  non-priority  class  for  the  same  service  distribution 
is  given  by  [NEE65] , 


H2{>t)  = P 'Ple~0U  + f(?)dr,  for  p2  > p{ 

L P-P,  2rr  J * 


P~P 


1~P  mi  -) 

= — f Hr)  Jr,  for  p2  < , 

2rr  < 1 


2rr 


(3.3a) 


(3.3b) 


where 


P = X/p  . 

P,  = Xj/pj  , 

a = |(p-p,)  (1"P)]/PP  , 

«1*  (I-/pT  >V 

q2=  ci+ysr  )2/i  1 

/(r ) - e'^yOft^rHr-ajVrfr-a), 


and 


//,  (/)  = H2  (/),  for  the  single-class  case. 

Since  both  distributions  are  for  exponential  head-of-the*line  priority  queueing  system,  the  analytic 
model  will  be  for  a M/M/ 1 system.  The  inversion  formula  for  the  M/G/ 1 system  is  not  available  in  the 
literature  and  its  solution  would  be  too  esoteric  for  this  discussion  and  would  not  be  any  more  informa- 
tive than  the  A//A//1  analysis.  R.  H.  Davis  [DA65]  has  developed  an  explicit  formula  for  the  waiting- 
time distribution  of  any  number  of  parallel,  negative-exponential  servers  subject  to  Poisson  demands  from 
any  number  of  priority  classes. 
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3.3  PRIORITY  ASSIGNMENTS 

The  priority  assignment  will  be  based  on  how  long  a message  has  been  in  the  system,  such  that  when 
Pr  (wait  in  Q > T)  the  message  is  given  a priority.  There  will  be  only  two  classes  of  messages;  those  of 
priority  (p  ~ I)  and  those  without  (/>  = 2).  The  probability  that  a message  reaches  some  age  T and  is  given 
a priority  will  be  established  from  the  “waiting  time”  distribution.* 

The  priority  process  to  be  considered  in  this  study  will  be  of  the  MjGj  1 type,  where  the  service  dis- 
tribution is  a constant.  The  constant  service  distribution  is  a result  of  all  messages  being  of  equal  length 
and  having  identical  service  requirements.  The  constant  server  requires  analysis  by  the  Imbedded  Markov 
chain.  Also,  the  analysis  for  a priority  process  has  some  additional  difficulties,  since  the  state  description 
has  to  incorporate  information  about  each  priority  class,  and  thus  becomes  multidimensional. 

The  waiting-time  distributions  for  a M/M/\  system  will  be  used  to  get  an  analytic  approximation  to 
the  M/G/ I simulation. 

The  selection  of  priority  assignments  in  this  investigation  is  based  upon  the  waiting-time  distribution 
as  opposed  to  service-time  distribution.  Numerous  authors  have  addressed  the  problem  of  priority  selection 
based  o.i  minimum  service  time  (e.g.,  [Phipps  56) . [Jackson  60  and  62j . and  [Morse  5*1 ).  Cole  [007 i] 
based  his  priority  selection  upon  message  length,  yielding  the  same  results  as  a selection  based  on  service 
time. 

Conway,  et  al  [CON67J  asserts  that  the  optima!  priority  assignment  for  different  service  time  require- 
ments would  be: 


"l  "2  ' "n 


where 


- random  variable  for  service-time  of  a k th  class  message, 
r.j.  = some  linear  cost  or  other  weighting  factor  associated 
with  the  waiting-time  for  messages  of  k 1,1  class. 

Cox  and  Smith  [Cox6l]  assert  there  is  some  distinction  to  be  made  regarding  a message's  queueing- 
time and  its  waiting-time.  The  former  being  the  time  from  his  entering  the  system  to  the  instant  at  which 
his  service  begins,  and  the  latter  is  the  time  f;om  entering  the  system  to  the  instant  at  which  his  service  is 
completed.  Waiting-time  equals  queueing  time  plus  service  time.  In  many  systems  queueing-time  and 
service-time  are  statistically  independent,  allowing  one  to  obtain  properties  of  waiting  time  from  those  of 
queueing-time,  and  vice-versa. 

* | (ok*  71  I investigated  the  tjse  where  priority  assignment  was  based  or.  message  length  distribution. 
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Message's  queueing-time  or  waiting-time  are  of  direct  interest  when  there  is  an  economic  loss 
(message  expires)  if  it  is  kept  queueing.  The  effect  of  too  long  a delay  either  quereing  or  waiting  results 
in  giving  a priority  assignment  to  a message  once  it  has  been  in  the  system  more  than  km  time  units.  The 
effect  of  this  priority  assignment  will  hopefully  expedite  delivery  of  messages  that  have  already  spent 
ample  time  in  the  system. 

The  introduction  of  a priority  system  of  any  form  has  no  effect  on  the  proportion  of  time  for 
which  the  server  is  busy.  This  is  solely  determined  by  the  total  service-time  of  the  messages  arriving  over  a 
long  period. 

The  priority  assignment  will  be  designated  as  p ( k ) and  defined  as 

p(k)  = Pr{Dr  > kMDTMAX),  0 < k < 1 . (3.5) 


where 


Dj  = total  average  message  delay  • 
k = priority  threshold  setting. 

3.4  MINIMIZING  THE  AVERAGE  MESSAGE  DELAY 

The  average  message  delay,  as  has  been  suggested  earlier,  is  one  of  the  essential  measurements 
deemed  of  vital  importance  in  a store-and-forward  communication  system.  The  work  of  Cobham  ICOB54] 
formulates  the  equation  for  finding  the  expected  waitirg  time  for  messages  of  each  priority  level  of  a single- 
channel (unsaturated)  wailing  line  system.  His  general  equation  as  would  apply  to  an  exogenous  priority 
discipline  has  the  form 


fo  '2  •ffl(0/(l-Op.,Xl-Op)  • 0.6) 


N 

p.  = A./p.  and  A = £ A.  , 
1 
k 

°p  »t  < '-"o  = °- 

1 

N 

B(t)  = (I/A)  £ \ 

1 


where 
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Bf(t)  is  the  cumulative  service  time  distribution  function  for  the  ith  priority.  Thus  for  a two  class 
priority  system  (p  = 1,2),  the  average  message  delay  for  the  priority  (p  = 1)  message  becomes 

(3.7) 


where 


K'o  = -i-X  J t2  dB  (()  = average  time  to  complete 
° existing  service. 

For  the  non-priority  (p  = 2)  messages,  the  delay  is 


W2  = / 0-P,)  (J-Pj-Pi)  - (3-8) 

When  priority  assignments  are  established  within  the  system,  then  the  priority  message  rate  becomes: 


\ = P(k)\ , (3.9) 

where  p (A:)  is  given  by  Eq.  (3.5),  and  the  non-priority  message  rate  is 


\2  = \{\-p(k)).  (3.10) 

The  unconditional  delay  at  an  arbitrary  node  becomes: 


where 


W = p (Ac)  + (1  -P(k))W2 


1 -P 


1~P(*-)P 

!-p, 


,P  < 1 


(3.11) 


1-p  = 1-pj -p2  . 

In  this  investigation,  all  messages  while  assumed  to  be  independent  and  are  of  equal  length,  ncnce 
have  equal  servicing.  When  pj  = /i2,  Kq.  (3.1 1 ) is  reduced  to 


W = W(FCFS). 


(3.12) 
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It  appears  at  first  glance  that  Eq.  (3.1 1)  cannot  be  minimized  by  some  optimal  selection  of  k in 
establishing  the  p ( k ) term.  Subsequent  developments  will  bear  out  that  the  average  message  delay  can  be 
reduced  over  the  entire  network,  but  cannot  be  directly  assessed  from  the  single  node  formulation  given 
by  (3.11). 

The  average  message  delay  over  the  network  is  a function  of  the  routing  procedures,  network  utili- 
zation factor  p,  and  network  topology.  The  allowance  of  an  adaptive  routing  procedure  vs  fixed  routing, 
makes  the  indirect  effects  of  priority  assignment  mathematically  intractable;  however,  there  is  a technique 
to  analytically  assess  the  effect  of  priority  assignment  designated  as  the  throughput  factor. 

3.5  THROUGHPUT  FACTOR  0 (r,  k) 

The  thrc  jghput  factor  0(/,  k)  will  be  defined  as  the  ratio  of  the  delivered  messages  to  the  total 
number  of  messages  handled  by  the  network.  Messages  finding  '.he  queue  full  (balking)  will  not  be 
included.  It  is  desired  to  obtain  the  optimum  value  of  k that  maximizes  0(r,  k)  for  all  p < 1. 

A closed  form  expression  containing  the  parameters  of  interest  can  be  obtained  from  a simplified 
network.  Fig.  (3.1)  represents  the  3-node  tandem  network  under  consideration.  Non-priority  messages  with 
lengths  drawn  from  an  exponential  distribution  with  mean  p - 1,  will  arrive  at  nodes  M and  S at  rates 
X,'  and  X,'  respectively.  It  is  assumed  that  the  inter-arrival  times  of  messages  at  either  node  are  indepen- 
dent. 


S M D 

Mg.  3.1  - Three  Node  Tandem  Net 


We  will  assume  that  1/2  of  X2'  messages  have  M as  a destination,  and  the  remaining  half  have  D as 
a destination.  We  desire  to  see  what  messages  coming  into  At  have  D as  a destination. 

Neglecting  the  actual  servicing,  the  probability  that  when  a message  leaves  node  5 it  is  a priority 
message  is  given  by 

p(k ) = Pr(kMDTMAX  <MDTmax)  ,0  < k < 1 (3.13) 

and  using  the  distribution  given  by  Eq.  (3.1),  we  get 


p(k)  = Xj'cxp  |-(1-X2')fan| -Xj'exp 


(3.14) 
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We  define  another  variable  called  the  “geriatrics”  factor  g (m),  as  the  probability  of  a message  with 
source  S and  destination  D,  expires  before  it  reaches  M. 

Then 


g(m)  * Pr(t  > M>rMAX)  , (3.15) 

and  since  all  messages  into  node  S have  no  priority,  we  can  use  H j ( t ) as  the  distribution  for  a single  class 
system,  i.e.,  let  p?  = 0.  Hence  Eq.  <3.1 5)  becomes: 

g(m)  = l-//j  (m)  - X2'  c (3.16) 

Fq.  (3.16)  shows  that  g(m)  is  a function  of  the  maximum  delay  time  m and  message  rale  at 
node  S.  Further  g (m)  is  not  a function  of  the  threshold  parameter  k , but  effects  the  initial  rate  of  un- 
delivered messages  as  a function  of  m and  X 2' . We  note  that 

p(k)  + g(m)  < 1 . (3.17) 


this  completes  the  analysis  at  node  S. 

At  the  second  node  M.  the  following  variables  are  defined: 

X j ' = arrival  rate  of  new  external  messages  at  node  M 
p ( k ) X2*  / 2 = arrival  rate  of  internal  priority  messages  at  node  3/  (3.18) 

(1  - p(k)  -£(w))X2 72  = arrival  rate  of  non-priority  messages  at  node  M. 

The  non-priority  message  rate  is  designated  by  X2  anti  given  by 

X2  - (!-/>(*)  ~*(m))X2 72  + X/  (3.19) 

and  the  priority  message  rate  is  designated  by  Xj  and  expressed  as 


Xj  = p { k ) X272  . 


(3.20) 
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The  total  arrival  rate  X of  messages  into  M becomes 

X = (1  -P{k)-g{m))\2'll  + X./  i p (A)X2'/2 

= X2'/2  - ^(m)X2,/2  + X/  . (3.21) 

We  further  note  that  when  the  selection  of  m and  the  corresponding  value  of  X^  is  such  that 

gim)\'l2  * 0 , 

then  (3.22) 

X * X^/2  + X/ . 

The  network  flow  pattern  is  as  shown  in  Fig.  (3.2).  To  calculate  the  probability  that  a message 
entering  at  node  M is  undelivered  at  node  D,  three  cases  must  oe  considered. 
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Fig.  3.2  - Flow  Pattern  in  a 3 Node  Net 


Case  1 

Messages  that  arrived  from  an  external  source  or  from  node  5 that  had  no  wait  at  node  S. 
We  define  the  following  random  variables  that  are  assumed  to  be  independent 


T j = time  spent  at  node  S 
r2=  time  spent  at  node  M 
T - Tj  + T2  = total  time  spent  in  the  system. 


(3.23) 
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From  probability  theory  [PARZEN  62]  we  can  express  the  following  conditional  probability  of  a 
random  event  A,  given  a random  variable  X by  the  Stieltjes  integral  as 

P(A)  = f_  Pr(A\X  = x)  clFx  (x) , (3.24) 


then 


P(A)  = Pr(T  > m)  = £_p,  (0  JFT|  (()  , 


where 


(3.25) 


P,(0  = Pr(T  > m | r,  =0), 


and 


Ft\  (0“{i  til' for  any  rcal  t' 


Thus,  the  probability  that  a message  expires  before  reaching  D given  that  it  had  no  wait  ai  node  S 
becomes: 


since 


P(A)  = Pr{T2  > m) 
= 1 - //2  (m) 


(3.26) 


T - T2  for  Tx  = 0 , 


//2  (/)  is  the  waiting-tim;  distribution  of  the  non-priority  messages  in  an  exponential  FCFS  priority 
queueing  system  and  is  given  by  Eqs.  (3.3a)  and  (3.3b).  Hence  Eo.  (3.26)  becomes: 


P(A)  = 


P P c -0™  + / , for  p2  > p, 


P “P 


(3.27) 


/ , for  p2  < Pj 
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where 


/ = ^ f1  e~"  y/(a2 -r)(r  - a, ) /r  (r -a)dr  , 

2n  "'a, 

a * (P-P,)(l  -P)/P.  (3.28) 

«,= 

“2=  < 1 + v/^l  >2  ' 


Case  2 

Priority  messages  coming  into  node  M. 

Using  the  conditional  probability  given  in  Eq.  (3.24)  we  define  P{B)  as 
P(B)  = Pr(T  > m)  = p2  (1)  dFT^  (?) 


(3.29) 


where 


p2(/)  = Pr(T  > m | T,  =f) 

= Pr(T2  > m - 7',  | F,  = r)  (3.30) 

= f>-(r2  > m - ( | T,  = /)  . 

and  due  to  the  independence  of  T j and  7^,  Eq.  (3.30)  becomes: 

P2  (t)  = X(  ■Xl)(m_')  . (3.31) 


The  random  variable  associated  with  the  distribution  function  Fj  ( t ) in  Eq.  (3.29)  is  that  of  a 


priority  message  and  exists  only  in  the  interval  from  km  to  m.  The  probability  density  function  (pdf)  of 
a priority  message  is  given  by 

0 , t 

(I  -Xj'Je-d  'X2')' 


/t  (0  = 


I 


|e"0  -X2>m  -<l  -X2')m| 

0,  / < km. 


othc«wise 


(3.32) 


31 


Hence, 


(1  ^lMm-Oe-O-Xj  )t 

P(B)=  f : ; dt 

km  1 'X2  )km  )™  ] 


Xj(l-  X2,)e*(1’xi)w[e'(xi  " X2>A:m-e"(XrX2')ml 
(Xj  - ) [e“  • *2  >*m-  e**1 ’ X2  > *] 


, 0 < A:  < 1 


(3.33) 


We  note  that  Eq.  (3.33)  holds  for  all  values  of  0 < k < 1 and  P(B)  - 0 for  k = 1,  as  can  be  noted  by  the 
limits  on  the  integral. 


Case  3 

A non-priority  message  that  had  to  wait  at  node  S. 

We  define  P(C)  as: 

p(0=/nr>m)=  J p3(f)dFT3(f), 

where 


Pj(0  = Pr(T>m  | Tx  » f), 
= jRr^  > m - r) , 


(3.34) 


(3.35) 


and 


10,  t>  km 

(1  - Xj')*'  ^ ” *2  [X2(l- X2  )^m)]  , otherwise  . 

0,  f < 0. 


The  random  variable  associated  with  the  pdf  given  in  Eq.  (3.34)  exists  only  in  .‘he  interval  specified  by 
Eq.  (3.35). 

Now 


thus 


p3(t)  * Pr(r2>m-t ), 


/>2-p, 

f~a(m  ■ 0+  /t  for  p^  > p 

p3(f> . ^ p-p,  1 

/,  ior  , ■C  p.  , 


(3.36) 
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where  / is  defined  by  Eq.  (3.28). 

The  solution  for  P(C)  as  given  in  (A.27)  and  (A.28)  allows  us  to  write  it  in  the  following  form: 


PiO  -* 


/j  + /2,  for  p2  > pj 
/2,  for  p2  < Pj 


0 < * < 1 


(3.37) 


where 


(1  - X2’)(P2-  p,)<f  “m 

(p  - Pj)(l  - Xj1-  a)(l  - ” ^2 


[1  -<f<>  ' ^-a)km]  , 


and 


I 


2 


(3.38) 


(1  - Aj'/vl-  P) 
2<I(1  - P*0  - *2 


-a2  e"mr  /(aj-rXr-a,) 

(1  - Xj' - r)[r(r-  a)) 


[1  _ e-(l  - X2'~  r)*m] 


The  throughput  factor  <t>(l,k)  is  a normalized  function  and  related  to  the  percent  oi  undelivered 
messages  v(t,k ) by  the  following: 


k)  = 1 - v{t,k).  (3  .39) 

v(t,  k)  is  the  percentage  messages  handled  by  the  system,  but  not  delivered.  From  Fig.  (3.2)  we  define  the 
probability  r(t,  k)  as  the  probability  that  a message  with  source  S and  destination  D is  undelivered.  Let  S(t,  k) 
define  the  probability  that  a message  with  source  M and  destination  D is  undelivered. 

Further,  let  u(k)  = - g{m)-p(k)).  By  Eqs.  (3.13),  (3.16),  and  from  the  definition  of  it  follows 

that  u(k)  is  the  probability  that  a non-priority  message  has  experienced  som»  ielay.  Then  r(t,k)  can  be 
expressed  by  the  following  mutually  exclusive  independent  events: 

r{t,k)  = g(r?i)  + (\-\2)P(A)+p(k)P(B)  + u(k)P{C).0<k<\,  (3.40) 


where 


(1  -X^)  = (1  - p)  = the  probability  that  a message  entering  the  system 
at  node  S finds  the  queue  empty  at  5. 


S(t,k)  = P(A)  ,0<*<  1 


(3.41) 
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The  percent  of  all  undelivered  messages  v(t,  k)  with  destination  D is 

» (t,  k)  = [(/\2'  /2)  r(t,  k)  + X/  S{t,  k ))  /(X2'/2  + X/  ) , (3.42) 

hence,  the  throughput  factor  0(/,  fc)  becomes 

<P«,k)  = 1 - [(X2' /2)Kr.A:>  ^ XirS;^Ar»/(X2'/2+Xi').  (3.43) 

It  should  be  noted  that  Fqs.  (3.40)  and  (3.41)  are  continuous  functions  of  k.  However,  it  was  shown  for 
Case  2 and  Case  3 that  P{B)  = 0 for  k - 1 and  P(0  = 0 for  k = 0. 

We  would  like  to  find  the  optimal  value  of  the  threshold  parameter  setting  k that  maximizes  4>{t,k)  for 
all  p < 1.  The  optimization  will  be  investigated  in  the  following  section. 

3.5.1  Numerical  Solution  for  {k  ) 

It  is  apparent  that  even  if  v(t,  k)  is  a convex  function  of  k , it  defines  k as  an  implicit  function  of  f, 
and  k cannot  be  readily  obtained  (if  at  all)  as  an  explicit  function  of  t.  Thus  it  becomes  necessary  to  employ 
a numerical  solution  to  find  k . 

F.q.  (3.42)  was  solved  using  Simpson’s  rule  [Krylon  621  to  evaluate  the  integrals  given  by  (3.27),  (3.33), 
and  (3.37).  Computer  runs  were  made  on  a fixed  value  of  X^  and  X,'  and  solving  (3.42)  for  values  of  k - 0.1, 
0.2,  . . .,  0.9,  and  some  appropriate  value  of  m = (3/Z)7'max).  The  selection  of  m was  chosen  to  make  %{m) 
reasonably  small,  i .c.,g{m)  < 5.0  percent.  All  runs  v,  ire  made  using  128  points  in  Simpson’s  rule  and  a value 
of  m = 20. 

Fig.  (3.3)  shows  the  fraction  of  undelivered  messages  v(t,k)  versus  priority  setting  k.  The  utilization 
factor  p was  calculated  from 


p = X2  /2  - g(m)  X2  /2  + Xj* 


(3.44) 


where 


£(m)  = X2  e (*  ~ ^2  ‘ 

When  m is  fixed,  g(m)  is  directly  effected  by  the  selection  of  X2  , hence  the  curves  given  in  Fig.  (3.3) 
depicts  the  influence  of  gim)  on  the  selection  of  k that  minimi  es  v(r,£).  The  range  on  p was  from  .37  to 
.49,  but  the  range  on  g(m)  was  from  35  percent  down  to  .02  percent.  The  results  indicate  that  when  the 
value  of  the  “geriatrics”  factor  g(m)  is  very  high,  there  is  no  well  defined  k value  to  minimize  v(t,  A),  as  can 
be  seen  by  the  flatness  of  the  bottom  curve. 
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Fig.  3.3  - Fraction  of  Undelivered  Messages  Vs.  Priority  Setting 


Table  (3.1 ) is  a tabulation  of  the  curves  in  Fig.  (3.3)  and  gives  the  associated  values  of  X,'  and  X2'  selected 
for  each  run.  It  should  be  noted  how  g ( m ) varies  with  the  selection  of  X2'  and  the  subsequent  deepening  of 
the  caternary  in  the  curves  as  g (m)  decreases.  Thus,  when  g ( m ) is  very  small  the  selection  of  kop  is  well 
defined  and  represents  very  fine  tuning  on  the  network  parameter.  The  reason  for  this  fine  sensitivity  is  due 
to  the  very  small  value  of  v ( t,k ) < .021  percent,  hence  the  selection  of  knp  represents  what  appears  to  be  a 
vernier  setting  on  v ( t.k ).  The  average  value  of  kop  in  this  set  of  runs  was  0.5.  The  curves  have  been  scaled 
against  the  ordinate  values  as  indicated.  Curves  are  assumed  numbered  from  bottom  to  top  as  1,  2 n. 


TABLE  3.1 

FRACTION  OF  UNDELIVERED  MESSAGES  v(t,k) 
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Curve  #1 

Curve  #2 

Curve  #3 

Curve  #4 

Curve  #5 

k 

>• 

ii 

o 

V = .i 

>* 

n 

b 

<-ri 

V’-i 

X,'  = .1 

x2'  = .95 

X2-‘  = .8 

X2  = .85 

K - 7 

X2  = .6 

0.1 

.3216 

. 1 248E-01 

.4037E-01 

.1527E-02 

.1768E-03 

0.2 

.3205 

. 1 244E-01 

.3957E-01 

.1444E-02 

.1643E-03 

0.3 

.31% 

.1221E-01 

.3907E-01 

.1415E-02 

.16 19E-03 

0.4 

.3189 

. i 2 ! IE-01 

.3875E-01 

1410E-02 

.1633E-03 

0.5 

.3183 

.1210F.-01 

.3859E-01 

.1424E-02 

.1676E-03 

0.6 

.3178 

.1219E-01 

.3858E-01 

.1456E-02 

.1751E-03 

0.7 

.3177 

. 1 243F.-01 

.:.880E-0l 

.1515E-02 

.1868E-03 

0.8 

.3179 

.1289E-01 

.3944E-01 

.1617E-02 

.2045E-03 

0.9 

1 

.3195 

.1380E-01 

.4097E-01 

.1789E-02 

.23 10E-03 

Fig.  (3.4)  are  the  curves  associated  with  Table  (3.2).  This  selection  of  runs  were  made  over  a wider  range 
of  network  intensity  p,  but  for  the  same  m = 20  as  before.  The  curves  display  the  same  trend  as  the  previous 
results,  indicating  a single  low  point  over  the  range  of  k.  Flat  curves  are  still  indicative  of  no  real  sensitivity  to 
the  selection  of  kQp  and  deep  caternaries  indicate  precise  location  of  kQp  for  minimum  v(t,k).  The  average 
value  of  kQp  for  this  run  was  0.375. 
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Fig.  3.4  - Fraction  of  Undelivered  Messages  Vs.  Priority  Setting 


TABLE  3.2 

FRACTION  OF  UNDELIVERED  MESSAGES  v(t,  k) 


k 

Curve  #1 

V ■ 4 

x2'  = a 

Curve  #2 
X/  = .1 
Xj’  = .8 

Curve  #3 
X/  = .7 

X2  " -4 

Curve  #4 

X.'  = -5 

X2'  = .4 

0.1 

.09692 

.01284 

.1253 

.2076E-02 

0.2 

.09332 

.01244 

.1241 

.1921E-02  : 

0.3 

.09134 

.01221 

.1238 

.1901E-02 

0.4 

.09071 

.01211 

.1239 

.1908E*02 

0.5 

09129 

.01210 

.1239 

.1917E-02 

0.6 

.0v302 

.01220 

.1239 

.1924E-02 

0.7 

.09599 

.01243 

.1240 

.1928E-02 

0.8 

.1003 

.01289 

1241 

. 193 1 E-02 

0.9 

.1063 

.01380 

.1241 

.1932E-02 
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Table  (3.3)  are  the  values  for  Fig.  (3.5)  and  indicate  that  family  of  curves  are  generated  for  the  same 
value  of  p,  and  in  this  case  for  p - 0.6.  It  can  be  observed  from  Eq.  (3.44)  that  there  are  numerous  selections 
of  X j*  and  X2  possible  for  a given  p,  the  only  restriction  being  that 

X/  + X2'/2<  1.  (3.45) 

This  selection  of  X^  and  X2  keeps  p < 1 on  each  link  of  the  3-node  net  given  in  Figs.  (3.1)  and  (3.2). 
Since  g(m)  > 0 in  all  cases,  this  allows  X}  + X2'/2  = 1,  without  violating  the  requirement  of  keeping  p < 1. 
The  average  value  of  kQp  for  this  set  was  0.4. 


TABLE  3.3 

FRACTION  OF  UNDELIVERED  MESSAGES  v{t,k) 


k 

Curve  #1 

V = 5 

V ■ 

Curve  #2 
X/  = .15 

X2  ^ 

Curve  #3 
xi'  = 2 

x2’  = .8 

0.1 

.2123E-03 

,%88E-01 

.1 183E-01 

0.2 

.2073E-03 

.9504E-01 

.11 1 IE-01 

0.3 

.2072F.-0.i 

.9418E-01 

.1081E-01 

0.4 

.2075E-03 

.9364E-01 

.1074E-01 

0.5 

.2077E-03 

.9338E-01 

.1084E-01 

0.6 

.2079E-03 

.9345E-01 

.1 1 10E-01 

0.7 

.2079E-03 

.9401E-01 

. » 157E-01 

0.8 

.2079E-03 

.9544E-01 

.1234E-01 

0.9 

.2079E-03 

.9859F.-01 

.1354F.-01 
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Fig.  3.5  - Fraction  of  Undelivered  Messages  Vs.  Priority  Setting 

Recalling  that  P(A)  defined  the  probability  of  undelivered  messages  at  D that  arrived  from  an  external 
source  or  from  node  5 had  no  wait  at  S (Case  1).  P(B } defined  the  probability  of  undelivered  messages  that 
were  in  a priority  class  coming  into  node  M (Case  2).  P(Q  was  for  non-priority  messages  that  had  to  wait  at 
node  S.  in  all  cases  tested,  P(A)  and  P(B)  decreased  with  increasing  k (not  tabulated),  while  P(C)  increased 
with  k.  From  the  definitions  given  it  is  clear  that  P{A)  < P{C)  and  P(B)  < P(C)  for  all  0 ^ A:  < I . 
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IV.  SIMULATION  ON  AN  8-NODE  HIGHLY  CONNECTED  NETWORK 

The  first  series  of  simulation  runs  are  for  the  8-node  network  of  Fig.  (4.1).  Appendix  B contains  the  flow 
diagrams,  Fortran  listing,  explanation  of  variables,  and  listing  of  the  simulation  runs  used  on  the  selected  network. 

RAND  program  originally  had  no  queueing  or  preassigned  links  (PAL)  and  generated  messages  such  that 
the  network  had  constant  loading.  However,  messages  also  had  a dependency  on  whether  a given  node  was  busy, 
hence  message  volume  (MV')  would  not  always  generate  the  same  total  number  of  messages. 

A buffer  was  included  in  the  RAND  urogram  to  facilitate  a common  basis  for  comparing  the  various  rout- 
ing schemes.  RAND  had  to  use  “Hot-Potato”  since  there  were  no  buffers  (queues).  Buffers  could  hold  up  to 
30  messages  at  each  node. 
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I ig.  4.1  Highly  Connected  8-Node  Network 

4.1  COMPUTATION  OF  MESSAGE  VOLUME  (MV) 

Message  volume  (MV)  is  expressed  as  a percentage  of  links  carrying  messages  per  unit  time.  The  number 
of  messages  generated  in  each  time  interval  is  equal  to  [AfF(No.  of  links  in  network)/ 100].  Count  the  number 
of  links  from  a given  net  topology  and  double  this  number  to  account  for  duplexing. 
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Thus  if  number  of  links  in  topology  - 12,  then  24  becomes  the  total  number  of  links  in  that  net. 

If  MV  - 10,  then  = |2.4|  =2  messages  per  unit  time  generated.  If  running  time  of  simulation 

has  T-  500  time  units,  then  total  number  of  messages  generated  during  that  run  is  1,000.  For  any  given  MV 
and  T , the  total  number  of  messages  generated  in  the  system  will  always  be  the  same.  However,  the  total  num- 
ber of  messages  handled  by  the  system  can  be  less  than  or  equal  to  total  number  of  messages  generated. 

For  the  8-node  net,  there  are  20  links  and  total  number  of  duplexed  links  = 40.  For  MV  - 100,  then 
= 40  messages  per  unit  time.  If  simulation  time  = 500  time  units,  there  would  be  40(500)  = 20,000 
messages  generated. 

Messages  are  generated  randomly  with  a random  source  5 and  random  destination  D such  that  D ^ 5. 
Messages  are  generated  on  each  computer  cycle  and  the  message  volume  is  directly  proportional  to  the  value  of 
MV  selected  when  initiating  the  run. 

Messages  are  allowed  some  maximum  circulation  time  MD1'max  and  if  they  are  not  delivered  before  that 
time  is  reached,  they  are  dropped  from  the  sy  stem  and  tallied  as  undelivered  messages.  Dropping  messages  from 
the  system  after  some  maximum  cycling  time  would  hopefully  relieve  network  congestion  and  allow  the  remain- 
ing messages  to  get  delivered.  The  longer  messages  remain  in  the  system  without  being  delivered,  the  greater  the 
average  message  delay  becomes. 

If  the  queue  at  S is  full  when  a new  message  i;  generated,  then  a busy  is  detected  and  the  message  is 
dropped,  and  counter  NZ  incremented.  It  should  be  noted  that  messages  finding  the  queue  full  are  different 
from  messages  undelivered.  When  a significant  number  of  messages  find  the  queue  full  (say  10%  of  total  mes- 
sages generated),  then  this  becomes  an  indication  of  the  network  approaching  saturation. 

The  total  number  of  messages  generated  for  any  given  network  loading  was  the  same  for  all  routing  algo- 
rithms on  any  specified  net  configuration.  Thus  one  can  add  the  number  of  messages  delivered,  as  tabulated  in 
the  message  throughput  Table  (4.3),  to  the  number  undelivered  (value  given  in  parenthesis)  and  obtain  the  total 
messages  handled  by  the  system.  The  difference  in  the  totals  is  attributed  to  the  number  of  messages  finding 
the  queue  full  (tabulated  in  Chap.  V)  and  the  messages  in  transit.  Hence,  for  light  net  loading,  messages  delivered 
plus  undelivered  will  be  about  the  same  for  all  program  routines.1 

The  program  has  an  initialization  (warm  up)  of  T-  MDT MAX  and  is  also  the  maximum  time  a message  can 
be  in  the  network.  Thus  for  a message  time  > MDTmax  = 50,  those  messages  are  considered  undeliverable. 

4.2  CBJFCTIVliS  OF  THK  SIMULATION 

The  objectives  of  this  simulation  study  was  to  find  an  adaptive  routing  algorithm  that  could  maximize  mes- 
sage throughput  while  minimizing  average  message  delay.  Since  the  routing  would  he  stochastic  in  nature  and 
essentially  utilize  information  about  the  state  of  the  network  via  node  delay  tables,  a systematic  analytical  method 
of  solution  was  not  possible.  Ideally  we  would  want  a routing  algorithm  that  could  adapt  to  network  topology 
changes,  either  due  to  congestion  at  some  given  nodes,  or  adjust  to  physical  changes  in  node  or  link  structure. 

I.  Also  (average  limej  x (messages  dclivercdl  + (A/DT^^ ) \ |number  of  undelivered  messages)  a constant  for  a giver. 

value  of  MV. 
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100 
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The  stochastic  techniques  rely  on  information  exchanges  between  nodes,  and  this  additional  traffic  reduces 
the  network  capacity  to  send  data.  Since  each  node  must  make  inquiries  of  its  neighboring  nodes  to  calculate 
new  delay  table  updates,  this  requires  some  type  of  processor  to  be  available  at  each  node.  Another  objective 
of  this  study  was  to  minimize  hardware  and  software  requirements  and  still  maintain  some  high  degree  of  oper- 
ating  efficiency. 

The  simulation  affords  a method  to  test  the  effectiveness  of  various  routing  schemes  and  possibly  how  to 
optimize  some  aspects  of  these  routing  strategies  to  increase  throughput.  The  initial  computer  runs  were  on  a 
small  scale  computer  (DDP-24)  and  a highly  connected  8-node  network  on  which  many  combinations  of  routing 
strategies  were  employed  (see  Appendix  B4.1  and  B4.2).  It  was  envisioned  that  after  these  initial  test  runs  were 
made,  selected  algorithms  giving  the  best  performance  would  be  tried  on  a 19-node  ARPA  network  and  run  on 
a C DC -3 800  computer: 

4.3  NO  LINK  OR  NODE  FAILURES 

Table  A 1 gives  a listing  of  some  specially  selected  routing  algorithms  as  applied  to  the  8-node  network. 

This  sampling  was  taken  to  show  the  strong  points  and  weaknesses  of  various  routing  techniques. 

With  no  destruction  of  the  network,  it  was  assumed  that  the  ideal  routing  tables  would  give  optimal  results 
up  to  moderate  network  loading.  Thus  the  no  upucting  (NU)  routine  was  used  as  a basis  to  compare  all  others. 

Routing  tables  were  initialized  b>  computing  the  minimum  path  matrix  on  an  empty  network.  These  ini- 
tializing algorithms  would  be  optimum  if  either  the  network  did  not  become  too  congested  or  there  were  no 
net  failures. 


TABLE  4.1 

PROGRAM  ALGORITHMS  FOR  DDP-24  SIMULATION 

1. 

Dynamic  Programming  (DY) 

2. 

No  Updating  (NU) 

3. 

Backward  Learning  (BL) 

4. 

Bi-Adaptive  (BA) 

5. 

Superposition  (SP) 

6. 

ARPA  Shortest  Que  + Bias  + Periodic  Updating  (ARPA) 

7. 

Last  M Nodes  Visited  (LMNV) 

8. 

Best  Three  Outgoing  Links  (BTL) 

= MOD] 

9. 

Aging  Parameter  (X%  of  MDT MAX ) 
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Fig.  (4.2)  depicts  the  data  plots  of  Table  (4.2).  These  were  selected  from  all  runs  made  in  Table  B4.1  as  being 
the  most  informative.  Routines  such  as  negative  reinforcement  (NR),  bi-adaptive  (BA)  and  no  updating  (NU), 
gave  identical  results  for  all  value ^ of  message  volume  (MV),  even  up  to  100%  network  loading.  The  quantity 
[2/5  MV]  is  a scaling  factor  to  i.void  plotting  distortions.  For  a 20*link  network  and  an  MV  given  as  a percent- 
age, the  fraction  2/5  appears  in  the  brackets  and  \x)  = greatest  whole  integer  less  than  X. 


Fig.  4.2  - Message  Volume  - ( 2/5  MVJ  Per  l/nit  Time 
Average  Message  Delay  for  8 Node  Net 

Backward  learning  appeared  to  be  not  too  effective.  This  was  largely  due  to  the  “ping-pong”  effect  between 
successive  nodes.  Backward  learning  utilizes  no  direct  infomiation  on  the  progress  of  messages  passing  through 
the  system,  and  its  indirect  method  of  computing  delay  values  is  often  misleading  or  inadequate  or  both.  The 
drawback  on  the  BL  technique  will  become  more  apparent  as  it  tries  (unsuccessfully)  to  adapt  to  network  top- 
ology changes. 

Superposition  (SP)  of  NR  and  BL,  where  BL  can  increase  or  decrease  table  entries,  shows  only  a modest 
improvement  over  BL  as  far  as  average  message  delay  is  conceined.  Table  (4.3)  shows  the  message  throughput  for 
the  various  routing  algorithms  and  the  quantities  in  parentheses  indicate  the  number  of  undelivered  messages. 

It  can  be  noted  from  Table  (4.3)  that  the  SP  loutine  eliminates  the  number  of  messages  noi  being  delivered, 
and  as  a consequence  the  total  number  of  messages  delivered  increases.  The  average  message  delays  for  the 
ARPA  and  ARfA  *■  MOD j + 50%  2 aging  indicate  substantial  gains  in  both  minimizing  average  delay  and  maxi- 
mizing throughput.  Fig.  (4.3)  in  conjunction  with  Tables  (4.2)  and  (4.3)  illustrate  the  effectiveness  of  these  routines. 

The  curves  in  Fig.  (4.3)  were  plotted  against  the  no  updating  (NU)  routine  to  indicate  how  an  adaptive 
routine  enhances  message  delays  and  throughputs  when  the  message  volume  increases  substantially.  The  ARPA 
routine  was  comprised  of  a shortest  queue  + bias  + periodic  updating  technique;  although  this  was  a sizable  in- 
crease in  programming  and  bookkeeping,  the  gains  were  significant  enough  to  warrant  such  an  investment.  The 


2.  Aging  parameter  same  as  threshold  setting  k in  Chap.  III. 


AVERAGE  TIME  FOR  MESSAGES  IN  AN  8-NODE  HIGHLY  CONNECTED  NETWORK 


MESSAGE  THROUGHPUT  IN  AN  8-NODE  HIGHLY  CONNECTED  NETWORK 


modification  to  the  ARPA  routine  (LNV,  BTL  and  50%  aging)  indicated  some  further  improvement  as  regards 
to  message  delay  and  throughput. 
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Fig.  4.3  - 8-Node  Network  With  50%  Disabled  Links 

After  all  ;.te  routines  listed  in  Appendix  B4.1  were  run,  it  became  apparent  as  to  what  rouhnes  were  the 
most  effective.  None  of  the  modifications  when  applied  to  the  BTL  doctrine  would  significantly  increase  its 
effectiveness.  The  various  combinations  of  last  M nodes  visited  and  aging  parameters  were  applied  to  the  ideal 
routing  tables  (NU)  anti!  the  values  that  “peaked”  our  performance  measures  were  selected  as  being  optimal.  The 
combination  of  last  node  visited  (LNV)  and  best  3 output  links  (BTL)  gave  the  minimum  average  message  delays 
and  maximum  message  throughput.  This  combination  of  LNV  and  BTL,  henceforth  called  MUD j , was  applied 
to  all  the  routing  strategies  listed  and  resulted  in  joine  significant  improvements. 

The  aging  parameter,  whereby  at  some  percentage  of  A 1T)T  MAX  a message  is  given  priority  on  the  queues, 
peaked  both  message  delay  and  throughput  at  a value  of  50%/  Since  there  was  no  direct  method  ( f calculating 
the  optimal  values  of  LMNV  and  aging  parameter,  it  had  to  be  determined  via  simulation.  Once  these  optimal 
values  were  determined  they  were  applied  to  the  better  routing  strategies  such  as  NU  and  ARPA  and  these  be- 
came the  adjusted  routines  given  in  Table  (4.2)  zsMOD^  + 50%  w<\ARPA  + MOD^  + 50%. 

The  preliminary  conclusion  would  be  that  the  combination  of  LNV  + BTL  = MODx  when  applied  with 
the  50%  aging  parameter  should  yield  the  minimum  message  delay  and  maxim;  . throughput  for  the  network 
under  consideration.  It  was  further  felt  that  <i  .his  technique  were  applied  to  any  of  the  listed  routing  algorithms 
of  Appendix  B,  it  should  enhance  their  perfonnan''-*  This  enhancement  has  been  effectively  verified  from  the 
tabulated  results  of  Tables  (4.2)  and  (4.3).  The  next  Sicp  in  confirming  the  optimal  adaptive  routing  technique 
would  be  to  see  how  it  fares  when  applied  to  a damaged  work. 


3. 


50%  aging  was  dose  to  the  predicted  value  for  optimal  k given  in  Chap.  111. 
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4.4  SIMIJL  :a\  WITH  50%  LINK  DESTRUCTION 

The  1 follow  through  on  the  8-node  highly  connected  network  would  be  to  see  how  the  adaptive 
routines  adjust  ,o  high  network  destiuction.  The  destruction  applied  only  to  the  links  (50%  outage)  and  was 
used  to  simulate  a highly  disasterous  situation  where  half  of  the  network  links  (or  lines)  were  inoperative. 

Figure  (4.3)  shows  the  connectivity  with  50%  link  destruction  on  the  8-node  network.  This  level  of  destruc- 
tive was  just  short  of  disconnecting  the  net.  Table  (4.4)  lists  the  same  set  of  selected  runs  as  were  given  for  the 
no  link  or  node  failure  case.  Although  there  was  an  extensive  set  of  runs  made,  as  indicated  by  Table  (B4.2)in 
Appendix  B,  the  algorithms  selected  represented  a good  cross  section  in  terms  of  our  performance  measurements. 
The  dynamic  (DY)  algorithm  would  be  the  presumed  theoretical  optimum  in  this  case,  in  that  it  would  imme- 
diately sense  link  or  node  failures  and  recompute  the  routing  tables.  The  intensity  of  network  traffic  was  in- 
creased to  the  point  where  the  number  of  messages  finding  the  queue  full  (not  tabulated)  exceeded  10%  of  the 
total  messages  generated  and  this  was  found  to  be  for  a MV  = 40;  hence,  [1/5(40)]  = 8 became  the  maximum 
intensity  that  was  used  for  this  set  of  runs.  The  identical  set  of  links  were  disabled  for  all  iterations  of  any  given 
routing  algorithm,  thus  all  algorithms  were  compared  for  identical  network  topology  and  message  intensities. 

K is  apparent  from  Table  (4.4)  and  the  curves  of  Fi^.  (4.4)  that  backward  learning  gets  aimos*  hopelessly 
bogged  down  and  this  is  reflected  in  its  average  message  delays  and  throughputs.  Table  (4.4)  also  iljustrates  how 
the  no  updating  (NU)  algorithm  fails  to  adjust  *o  the  net  loading,  since  after  initializing  its  tables  they  are  no 
longer  updated.  The  ARPA  routine  appears  to  be  fairly  adaptive  to  the  network  topology  and  loading  effects. 
The  modifications  to  the  ARPA  routine  (LNV  + BTL  + 50%  Aging)  improved  the  performance  of  the  routing 
strategies.  We  should  taka  cognizance  of  (ARFA  ♦ MOD j + 50%)  routine  in  the  area  of  MV  = 6 and  7,  which 
represented  moderate  loading.  The  MOD  j + 50%  addition  slightly  surpassed  the  DY  routine  in  minimum  average 
message  delays;  hence,  DY  is  not  a theoretical  lower  bound  (see  Table  (4.5)). 


I ig.  4.4  - Message  Volume  = (1/5  MVJ  Per  Unit  Time  Average  Message 
Dehy  lor  Disabled  8-Node  Net 


AVERAGE  TIME  FOR  MESSAGES  ON  AN  8-NODE  NETWORK  WITH  50%  LINK  DESTRUCTION 


48 


TABLE  4.5 

MESSAGE  THROUGHPUT  FOR  AN  8-NODE  NETWORK  WITH  50%  LINK  DESTRUCTION 


(1/5  MV] 

DY 

1 

BL 

ARPA 

ARPA  + 
MODj  + 50% 

1 

500  (0) 

471  (227) 

300  (0) 

500  (0) 

2 

1000  (0) 

790  (388) 

1001  (0) 

1002  (0) 

3 

1500  (0) 

939  (500) 

1505  (0) 

1501  (0) 

4 

2001  (0) 

1180  (414) 

2002  (0) 

1999  (0) 

c 

2497  (0) 

1282  (484) 

2497  (0) 

2496  (0) 

6' 

2997  (0) 

1016  (588) 

3003  (0) 

3008  (0) 

7 

3 4*4  (0) 

1207  (501) 

3460  (0) 

3480  (0) 

8 

’*>8  (0) 

*-  

1053  (570) 

3499  (4) 

3512  (0) 

The  DY  and  ARPA  routines  did  not  have  the  congestion  control  afforded  by  the  50%  aging  parameter  that 
gave  priority  to  messages  already  in  the  system  and  was  set  at  0.5  MDTmax.  This  congestion  control,  although 
not  directly  a routing  algorithm,  indirectly  enhances  some  of  the  delay  table  updates  by  removing  some  of  the 
node  congestion. 

4.5  ROUTING  TABLES 

A look  at  some  of  the  selected  routing  tables  may  give  some  insight  into  the  effectiveness  of  the  various 
routing  algorithms.  Initially  all  routing  tables  will  have  the  same  values.  When  message  intensity  starts  increas- 
ing, the  various  adaptive  routines  will  update  the  tables  and  depending  on  the  effectiveness  of  these  routines 
the  new  routing  tables  will  or  will  not  yield  optimal  path  assignments.  This  example  also  illustrates  why  the  BL 
routine  has  such  a high  percentage  of  undelivered  messages.  Tables  for  higher  message  intensities  continue  to 
show  the  futility  of  BL  routing  The  tables  for  the  ARPA  routine  along  with  MOD^  + 50%  show  almost  iden- 
tical routings  for  the  same  network  loading.  We  used  the  extreme  cases  to  make  the  illustration  more  informa- 
tive. 

It  should  be  reiterated  here  that  BL  was  duely  noted  by  its  creators  (B066J  as  being  something  less  than 
optimum.  They  introduced  NR  and  SP  as  a means  to  counter  the  undesirable  effects  of  BL.  The  use  of  BL 
in  this  study  was  to  dramatize  the  fact  that  a routing  algorithm  although  adaptive  may  not  be  very  effective. 
Hence,  our  concern  is  not  only  for  adaptive  routines  but  effective  ones,  when  measured  against  the  performance 
criterion. 
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We  will  define  our  routing  tables  with  the  following  notations: 

A'ij,  k)  = routing  tables  at  node  /,  with  destination  /, 
via  node  k 


and 


Al(j,  k)  = A 1 
A2(j.  k ) = A2 


thus 


An(i,  k)  = An 


A - (/l1,  A2,  ...  , An)  = routing  matrix 


An  example  of  bow  the  BL  routine  gets  bogged  down  with  very  light  network  loading  can  be  quickly 
seen  by  comparing  its  routing  tables  to  those  of  the  DY  routine.  The  tables  are  used  as  follows. 

Assume  a message  is  at  some  node  / and  has  a destination  node  /,  where  / #=  i.  Then  go  to  the  appropriate 
table  for  A’ij,  k)  and  select  the  minimum  (j,  *)  element  value,  this  gives  the  next  intermediate  node  k.  If 
k * j,  then  continue  until  k ~ j and  this  will  be  the  final  destination  node.  If  there  are  more  than  one  minimum 
(i,  k)  value,  take  any  one  and  proceed  until  destination  is  reached.  The  final  tally  of  intermediate  nodes  will 
give  the  delay  time  for  the  message  in  going  from  / to  /. 

An  example  will  quickly  demonstrate  the  above  procedures.  Take  node  1 as  a source  node  and  node  4 
as  destination.  Fig  (4.3)  shows  these  two  nodes  to  be  neighboring  nodes,  and  hence  a message  should  take 
only  one  time  unit  io  go  from  1 to  4.  From  Table  4.6  the  DY  routine  gives  the  following  routing: 

for  A ‘{j,  k ),  where  i = 1 and  / = 4 
A'(4,k)=  1, 
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and  this  corresponds  to  a k - 4,  and  k - j - 4,  hence  message  is  at  destination  in  I time  unit. 

Table  4.7  gives  the  BL  routing  for  the  identical  network  loading,  and  its  routing  tables  direct  a message 
from  node  1 to  4 as  follows: 


Al{4,  k)  - 2=*k  = 5 #/  (continue  until  k = /') 

/1S(4,  k)  = 3=»*  = 6 or  7 
A6(4,  k)  = 2^>k  = 2 or  3 
/42(4,  /c)  = 3=>fc  = 6 
A6{ 4,  A)  = * 2 or  3 

/43(4, 

Message  is  now  looping  between  nodes  2,  3,  and  6 and  will  not  get  delivered  since  there  is  no  alternate 
route  to  take.  If  at  node  5,  the  message  had  taken  7,  then  it  would  have  traveled  the  following  route. 


Al(4,  *)-2-*-5#/ 
i45(4,  k)  ■ 3+k  - 6 or  7 
A1  {A,  k)  ■ 6 +k  -5  or  6 
A6( 4,  *)  ■ 2~k  * 2 or  3 

and  the  next  set  of  nodes  are  the  same  as  the  previous  set,  except  this  time  there  Is  an  additional  entry  at 
k ■ 7,  and  the  message  will  still  not  get  delivered. 

Thus  we  can  see  how  some  routines  although  adaptive  can  result  in  very  poor  performance  as  regards  to 
average  message  delay  and  throughput. 
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TABLE  4.6 

ROUTING  TABLES  FOR  DY  ROUTINE  WITH  50%  DISABLED  LINKS  AND  [1/5  MV]  * 5 


A'a.k) 

A2(j,  k) 

A: '(/.*) 

a\  k) 

N 

k4  5 

N 

C 3 6 

\ 

k 2 6 8 

\ 

1 1 8 

1 

0 0 

i 

4 3 

1 

4 : 

J 3 

1 

1 3 

2 

4 3 

2 

0 0 

2 

l : 

2 3 

2 

4 3 

2 

3 3 

3 

1 2 

3 

0 0 0 

3 

4 2 

4 

1 3 

4 

3 4 

4 

4 < 

* 2 

4 

0 0 

5 

3 1 

5 

3 2 

5 

3 : 

2 4 

5 

2 4 

6 

4 2 

6 

2 1 

6 

2 

1 3 

6 

3 3 

7 

4 2 

7 

3 2 

7 

3 : 

2 4 

7 

3 4 

8 

2 4 

8 

2 3 

8 

3 ; 

3 1 

8 

3 1 

A5(j,  k ) 

460.  k) 

/ 

I7 (A  *) 

,18C.  *) 

/\ 

1 6 

2 3 

5 

7 

\ 

i 

X 

5 6 

/\ 

3 4 

1 

1 3 

1 

4 4 

2 

3 

1 

2 3 

i 

4 2 

2 

4 2 

1 2 

3 

3 

2 

3 2 

2 

2 4 

3 

4 2 

3 

2 1 

3 

3 

3 

3 2 

3 

1 3 

4 

2 4 

4 

4 3 

3 

4 

4 

3 4 

4 

3 1 

5 

0 0 

5 

3 3 

1 

2 

5 

1 2 

5 

3 3 

6 

3 1 

6 

0 0 

0 

0 

6 

2 1 

6 

2 4 

7 

3 2 

7 

3 3 

2 

1 

7 

0 0 

7 

3 4 

8 

3 3 

8 

3 2 

4 

4 

8 

4 3 

8 

0 0 
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TABLE  4.7 

ROUTING  TABLES  FOR  BL  ROUTINE  WITH  50%  DISABLED  LINKS  AND  (1/5  MV]  = 5 
Al(j,  k)  A2(j,k)  A V/,  k)  A*(j,k) 

8 

37* 

6 
19 
44 

41 
2 

48 

42 


A^fi.k)  Ab(j,k)  A7fi.k)  A*(j.  k) 

/\  1 6 7 2 3 5 7 _ 5 6 j 

1 1 34  34  1 44  38  5 40  1 40  80 

2 25  15  40  2 6 16  32  46  2 27  9 

3 27  20  39  3 14  8 28  15  3 16  3 

44  3 3 4 2253  4 66 

5 20  23  12  5 20  28  1 19  5 8 2 

6 47  20  46  6 44  41  7 32  6 22  16 

7 44  49  21  7 44  44  3 44  7 43  18 

8 37  14  2 8 9 4 2 20  8 15  10 
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4.6  SUMMARY  OF  8-NODE  SIMULA!  ION 

The  8-node  simulation  on  the  DDP-24  was  essentially  employed  to  check  out  some  of  the  selected 
routing  algorithms.  It  should  be  further  noted  that  not  all  available  adaptive  routines  were  tested,  but  instead 
some  representative  sampling  of  some  known  algorithms  were  used  to  demonstrate  strong  points  and  weaknesses. 

The  most  desirable  features  sought  for  in  an  adaptive  routing  algorithm  would  be  as  follows. 

1.  Minimum  hardware  requirement  at  each  node. 

2.  Minimum  software  and  inter-node  traffic. 

3.  Adaptability  tc  changes  in  network  topology  and  traffic  intensity. 

4.  Maximum  message  throughput. 

5.  Minimum  average  message  delay. 

6.  Minimum  pumber  of  undelivered  messages. 

It  is  very  conceivable  with  present  day  technology  that  with  unrestricted  use  of  items  1 and  2,  a very 
effective  adaptive  routine  could  and  may  already  be  developed.  R. taxing  hardware  and  software  requirements 
would  certainly  add  to  the  initial  overhead  and  general  operating  cost.  Item  3 is  really  the  crux  of  problem,  for  were 
it  not  for  the  changing  of  net  topology  (either  link  or  node  outages)  and  if  traffic  intensity  were  always  moderate 
p«  1,  then  any  effective  minimum  path  algorithm  would  certainly  suffice.  The  performance  measures  of 
items  4 and  5 arc  indicative  to  this  study  and  may  differ  for  different  networks.  While  item  6 may  be  toler- 
ated if  very'  small  in  some  networks,  a network  dedicated  to  emergencies  or  military  use,  etc.,  may  not 
tolerate  even  I lost  message 

Our  initial  runs  on  the  small  scale  computer  and  8-node  network  gives  us  some  insight  into  the  selection 
of  algorithms  to  be  tried  out  on  the  19-node  ARPA  network  using  a CDC-3800.  It  appears  that  the  technique 
of  last  node  visited  (LNV)  and  best  three  output  links  (BLT),  when  combined  (MODj)  yields  some  effective 
results.  The  additional  feature  of  the  aging  parameter  appears  to  “peak”  the  network  performance  when  a 
value  of  50%  is  employed. 

There  is  a distinction  to  be  made  between  throughput  (total  messages  handled  in  the  system)  and 
throughput  factor  0(/f  /cXfraction  of  messages  handled  and  delivered).  Table  (4.3)  shows  for  a (2/5  MV]  = 32, 
the  throughput  was  quite  different  for  various  routing  algorithms.  The  SP  had  a throughput  of  6038  messages 
and  a throughput  factor  0 (/.  k)  of  1.0.  However,  the  ARPA  and  ARPA  + MOD j + 50%  aging  had  over  twice 

the  total  throughput  as  the  SP  and  also  a 0(/f  k)  - 1.0.  Thus,  one  would  like  to  maximize  0(/,  k)  for  all 
p < 1,  plus  maintain  the  highest  total  throughput  rate.  (Generally,  when  the  average  message  delay  is 
minimized  and  0(r.  k)  = 1.0,  for  all  p<  1.  then  the  total  throughput  will  be  at  its  peak. 

The  dynamic  (I)Y)  routine,  although  not  a theoretical  optimum,  but  still  not  possible  to  employ  on  a 
real  network,  w ill  still  be  used  as  a basis  for  comparing  other  algorithms.  The  set  of  runs  for  the  19-node 
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ARPA  net  that  were  run  on  the  CDC-3800  are  tabulated  in  Appendix  C.  Note  that  runs  were  made  on  the 
ARPANET  with  a uniform  input  message  distribution  and  a Poisson  input  distribution.  Many  iterations  of 
the  selected  algoriui/ns  were  made  to  see  if  performance  still‘“peaked”  with  the  same  algorithms. 
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V.  CDC3800  SIMULATION  ON  A 19  NODE  ARPA  NETWORK 

The  same  rules  and  constraints  were  employed  in  these  simulation  runs  as  were  used  in  the  previous 
chapter.  The  network  was  a 19  node  ARPA  network  and  the  computer  a CDC-3800.  The  algorithms  were  run 
with  no  destruction,  then  with  10%  of  the  links  disabled,  and  finally  with  1 node  disabled.  Again  only  a few 
selected  runs  will  be  tallied  and  discussed,  since  inclusion  of  all  runs  would  make  for  a prohibitive  amount  of 
listing.  Appendix  C lists  the  runs  made  on  this  net  configuration. 

5.1  19  NODE  ARPA  SIMULATION 

Table  (5.1)  lists  the  algorithms  chosen  for  the  19  node  ARPA  network  of  Fig.  (5.1).  This  series  of  runs 
was  made  with  no  link  or  node  failures,  with  first  a uniform  input  message  distribution,  and  then  with  a 
Poisson  distribution. 
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TABLE  5.1 

PROGRAM  ALGORITHMS  FOR  CDC-3800  SIMULATION 

1 . Dynamic  Programming  or  No  Updating  (NU) 

2.  Backward  Learning  (BL) 

3.  No  Updating,  50  percent  Aging  (NU  + 50%) 

4.  ARPA  Shortest  Queue  + Bias  + Periodic  Updating  {ARP A) 

5.  MOD{  to  ARPA  (MOD{  + ARP) 

Fig.  (5.2)  shows  the  curves  for  the  tabulated  v.lues  of  average  message  delay  given  in  Table  (5.2).  The  NU 
routine  shows  a slightly  better  performance  than  the  MOO{ , but  the  throughput  with  the  MOD^ , given  in 
Table  (5.3)  was  slightly  better  at  high  traffic  intensity.  The  BL  routine  still  has  its  difficulties  as  shown  by  its 
average  message  delays  and  throughput  and  also  the  number  of  undelivered  messages  (undelivered  messages  are 
in  parenthesis). 


MESSAGES  PER  UNIT  TiME 

Fig.  5.2  - Average  Message  Dclaysfor  19  Node  Net  With  Uniform  Input  Distribution 
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TABLE  5.2 

AVERAGE  TIME  FOR  MESSAGES  WITH  A UNIFORM  INPUT  DISTRIBUTION 


MV 

NU 

Bl 

MOD{ 

NU  + 
ARPA 

ARPA 

ARPA  + 
MODl 

3 

2.50 

2.52 

2.50 

2.52 

2.56 

2.55 

6 

2.96 

5.19 

2.96 

3.77 

3.44 

3.31 

9 

4.57 

9.67 

4.57 

7.79 

6.30 

6.35 

12 

8.44 

9.65 

8.41 

9.93 

10.25 

10.89 

15 

11.73 

11.22 

11.07 

10.94 

11.33 

12.22 

TABLE  5.3 

MESSAGE  THROUGHPUT  FOR  MESSAGES  WITH  A UNIFORM  INPUT  DISTRIBUTION 


MV 

NU 

BL 

NU  + 

MOD]  + 50% 

NU  + 
ARPA 

ARPA 

ARPA  + 
MOD  ^ + 50% 

3 

750(0) 

750(0) 

750(0) 

750(0) 

750(0) 

750(0) 

6 

1502(0) 

1036(0) 

1502(0) 

1455(41) 

1499(1) 

1499(2) 

9 

2203(0) 

860(679) 

2203(0) 

1660(274) 

2023(94) 

1987(102) 

12 

2367(227) 

10KX67J) 

2365(278) 

1815(476) 

1608(428) 

1657(420) 

15 

2255(521) 

1124(619) 

2328(495) 

1826(669) 

1692(575) 

1610(597) 

The  ARPA  data  simulation  plotted  in  Fig.  (5.2)  was  slightly  aided  by  ihe  MOD{  + 50%  aging  addition  to 
it.  The  periodic  updated  ARPA  (ARPA  + PUD)  is  the  algorithm  of  interest,  since  it  is  used  on  an  existing 
network. 

Table  (5.3)  shows  that  a small  advantage  at  moderate  network  loading  (MV  = 12)  is  obtained  when  the 
MOD  j + 50%  is  used  with  the  ARPA  routine.  Table  (5.4)  shows  the  extent  of  network  loading.  At  zMV~  12, 
over  1/3  of  the  total  messages  generated  for  any  algorithm  were  not  able  to  get  into  the  system.  At  a MV  - 15, 
almost  50%  of  the  total  messages  generated  found  the  queue  full,  this  is  in  addition  to  the  number  of  undelivered 
messages.  Hence,  it  can  be  seen  that  with  a MV  - 12,  the  network  has  a high  message  intensity. 
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TABLE  5.4 

NUMBER  OF  MESSAGES  FINDING  THE  QUEUES  FULL 


Fig.  (5.3)  shows  the  plots  associated  with  Tables  (5.5)  and  (5.6)  when  the  message  input  process  was  Poisson. 

In  the  Poisson  case  messages  were  identically  distributed  with  an  average  message  rate  at  each  node  / of  \ and 

19 

the  total  message  rate  into  the  system  per  unit  time  was  I A.. 

/'  = 1 1 


0 4 0 !2  16  20  24 

MES  AGES  PER  UNIT  TIME 


Fig.  5.3  - Average  Message  Delay  for  19  Node  Net  With  Poisson  Input  Distribution 
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TABLE  5.5 

AVERAGE  TIME  FOR  MESSAGES  WITH  A POISSON  INP  DISTRIBUTION 


\ 

MV 

NU 

BL 

BL  + 50% 

BL 

MOD  + 50% 

_i_. - 

“ l58~  n 

5 

2.57 

2.58 

2.57 

2.57 

0.316 

10 

3.39 

6.99 

3.39 

3.39 

0.474 

14 

6.73 

8.48 

7.20 

7.20 

0.632 

18 

10.04 

10.14 

9.93 

9.93 

0.790 

23 

11.71 

10.36 

10.91 

10.91 

TABLE  5.6 

AVERAGE  TIME  FOR  MESSAGES  WITH  A POISSON  INPUT  DISTRIBUTION 


When  the  input  distribution  was  Poisson  the  average  number  of  messages  at  an  arbitrary  node 
i = X-=L{x.).  When  all  nodes  were  identically  distributed  /:(*)  - i for  the  19  node  ARPA  net  and  since 
there  were  68  duplexed  links  in  the  ARPA  net,  we  have 


68  MV 
100 


19X,. 


Thus,  we  can  gc.  valent  value  of  MV  to  go  with  thv  selected  X.  . 


The  plots  of  Table  (5.7)  was  an  attempt  to  improve  on  the  BL  routine  with  50%  aging.  Noteworthy,  was 
the  increase  in  throughput  and  resulting  decrease  in  undelivered  messages  as  tabulated  in  Table  (5.7).  The  NU 
r utine  was  still  used  as  a basis  for  comparison. 


TABLE  5.7 

MESSAGE  THROUGHPUT  FOR  MESSAGES  WITH  POISSON  INPUT  DISTRIBUTION 
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~1 

\ 

MV 

NU 

BL 

flL+50% 

NU  ARPA 

ARPA 

C . 1 

0.158 

5 

738(0) 

738(0) 

738(0) 

738(0) 

738(0) 

0.316 

10 

1528(0) 

965(320) 

1528(0) 

1413(65) 

1518(7) 

0.474 

14 

1974(77) 

713(779) 

1919(106) 

1634(317) 

1999(116) 

0.632 

18 

1913(337) 

657(998) 

2132(184) 

1876(462) 

1745(451) 

0.790 

23 

1767(484) 

565(1  152) 

1844(470) 

1854(626) 

1651(630) 

The  updated  ARPA  routine  was  further  adjusted  by  a DP  = 2 factor.  The  DP  factor  is  a constant  bias 
that  is  applied  to  reduce  looping.  A more  e\tensive  treatment  of  the  effectiveness  of  the  DP  bias  is  given  in 
Fultz  [FU72].  Fultz  further  suggested  the  optimum  routine  as  far  as  his  study  was  concerned  was  the  short- 
est queue  + bias  + periodic  updating  (SQ  + BIAS  + PUD ) and  this  algorithm  was  the  one  extensively  used  on 
the  19  node  network  as  a basis  for  comparing  the  various  MOI)^  + X%  aging  modifications.  Fig.  (5.3),  along 
with  Tables  (5.6)  and  (5.8),  indicate  there  is  some  slight  improvement  in  average  message  delay  and  throughput 
respectively  when  MOD j + 50%  iging  was  applied  to  the  ARPA  algorithms. 


TABLE  5.8 

MESSAGE  THROUGHPUT  FOR  MESSAGES  WITH  A POISSON  INPUT  DISTRIBUTION 


\ 

MV 

NU  ARPA 
+ MOD , 

NU  ARPA 
+ MODx  + 50% 

ARPA  + 
MOD]  +50% 

0.158 

5 

738(0) 

738(0) 

738(0) 

0.316 

10 

1493(25) 

1527(3) 

1528(0) 

0.474 

14 

1917(155) 

1940(137) 

2199(1) 

0.632 

18 

1833(374) 

1729(395) 

1639(462) 

0.790 

23 

1771(564) 

1693(582) 

1796(545) 

5.2  ARPA  NETWORK  WITH  3 LINKS  DISABLED 

Table  (5.9)  lists  the  average  message  uelays  (sec  Fig.  5.4)  associated  with  3 links  (10%)  disabled  and 
these  delays  are  plotted  in  Fig.  (5.5).  The  dynamic  (D  Y)  routine  was  used  as  a standard  in  this  case  since 
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with  link  outages  there  would  have  to  be  some  adjustments  made  on  the  routing  tables.  The  ideal  roi,!»«g 
tables  used  with  no  updating  (NU)  would  not  be  an  unbiased  comparison,  for  it  would  not  adjust  its  »es  for 
the  disabled  links.  Tables  (5.10)  and  (5.1 1)  indicate  that  the  total  throughput  was  lowered  considerably  for  all 
the  routing  algorithms. 


TABLE  5.9 

AVERAGE  TIME  FOR  MESSAGES  ON  19  NODE  ARPA  NETWORK  WITH  THREE  OF  LINKS  DISABLED 
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Fig.  5.5  - Average  Message  Delay  on  19  Node  Net  With  Three  Links  Disabled 


TABLE  5.10 

MESSAGE  THROUGHPUT  FOR  19  NODE  ARPA  NETWORK  WITH  THREE  OF  LINKS  DISABLED 


\ 

NMSG 

ARPA 

MOD x + ARPA 

MODx  + 50% 
+ ARPA 

0.053 

1.0 

266(0) 

266(0) 

266(0) 

0.105 

2.0 

482(5) 

487(1) 

488(2) 

0.125 

2.4 

592(4) 

594(2) 

592(1) 

0.145 

2.8 

667(6) 

665(8) 

671(3) 

0.165 

3.1 

619(128) 

625(115) 

625(114) 

0.200 

3.9 

595(257) 

564(291) 

676(203) 

TABLE  5.11 

MESSAGE  THROUGHPUT  FOR  19  NOTE  ARPA  NETWORK  WITH  THREE  OF  LINKS  DISABLED 


NMSG 

DYN 

BA 

DYN  + MODx 

TO3T 

1.0 

265(0) 

261(10) 

265(0) 

0.105 

2.0 

479(0) 

464(28) 

479(0) 

0.125 

2.4 

587(0) 

f6l(35) 

587(0) 

0.145 

2.8 

670(4) 

562(35) 

667(2) 

0.165 

3.1 

673(71) 

673(90) 

686(56) 

0.200 

■ J 

3.9 

626(225) 

227(237) 

633(228) 
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Fig.  (5.4)  shows  which  links  were  disabled  and  indicates  that  the  network  had  its  cross  link  capacity 
reduced  by  2/3.  Since  all  messages  have  a random  source  and  destination,  then  the  probability  of  message  being 
in  the  8 node  group  with  destination  for  the  1 1 node  group  using  link  (16,  8)  becomes  (8/19)  (11/18)  = .257. 
The  same  probability  holds  when  the  source  group  is  in  the  1 1 node  set  and  sends  messages  via  link  (8,  16)  to  the 
8 node  group. 

When  link  (16,  6)  is  used,  then  the  probability  of  an  arbitrary  message  coming  into  the  net  and  uses  link 
(16,  6)  in  either  direction  becor..es  (7/19)  (12/18)  = .246.  Since  link  (16,  8)  has  the  greatest  probability  of 
being  used,  it  establishes  the  limiting  *'alue  of  message  flow  that  can  be  tolerated.  Hence,  link  capacity  = 0.257 
x (number  of  messages)  and  since  Us  hnk  capacity  is  unity,  the  number  of  messages  per  unit  time  becomes 
1/0.257  = 3.89  messages/unit  time  as  the  upper  limit.  For  an  identically  distributed  19  node  network, 

= 3.89/19  = .205=*  p = 1,  or  saturation  ievel.  The  value  of  >.(  = .200  was  the  maximum  value  selected  for 
this  set  of  runs  to  avoid  overloading  the  network  and  Table  (5.9)  reflects  the  selection  of  X;.  used  on  this  network 
with  three  links  disabled. 

The  DY  + MOD  L + 50%  aging  \.wd  some  improvement  as  compared  to  ARPA  in  average  message  delay 
and  throughput,  but  as  previously  stated  the  DY  routine  requires  instant  knowledge  of  node  or  link  failures, 
hence  could  not  be  considered  as  a practical  routine  in  a real  network.  The  ARPA  + MOD { + 50%  aging 
showed  some  slight  improvements  in  throughput  and  also  a fewer  number  of  undelivered  messages. 

5.3  ARPA  NETWORK  WITH  ONE  NODE  DISABLED 

Fig.  (5.6)  ’mws  the  node  disabled  for  the  runs  tabulated  in  Tables  (5.12)  and  (5.13)  and  plotted  in 
Figs.  (5.7)  and  (5.8).  The  ARPA  routine  still  faired  quite  well  and  the  ARPA  + MOD j + 50%  aging  helped 
increase  throughput  (Table  5.15)  slightly  at  high  message  intensities,  but  the  ARPA  + MOD j performed  better 
in  this  case.  The  disabling  of  one  node  had  a noticeable  effect  on  message  throughput,  since  the  disabled  node 
would  service  half  of  the  cross  network  traffic.  Also,  by  disabling  node  17,  it  took  3 links  with  it,  resulting  in 
a 1/3  reduction  of  available  links  for  cross  net  traffic.  All  algorithms  were  adjusted  such  that  node  17  would 
not  send  or  receive  messages,  but  routings  that  previously  routed  thru  node  17  would  have  to  be  adaptive. 
Applying  50%. aging  to  BA,  NR,  and  SP,  aided  them  in  them  throughput.  Table  (5.14),  but  only  slightly  in  mes- 
sage delay  (Table  (5.12)  and  Fig.  (5.8)). 


TIME  UNITS 


llg.  5.7  - Average  Message  Delay  on  19  Node  w 

Net  With  One  Node  Disabled  ? 4 


0---0-0  ARPA  ♦ MOD, 

A-  -A-  -A  ARPA  ♦ MOD,  ♦ 50% 


4 6 8 

MESSAGES  PER  UNIT  TIME 


I ig.  5.8  Average  Message  Delay  on  19  Node 
Net  With  One  Nod?  Disabled 
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TABLE  5.12 

AVERAGE  TIME  FOR  MESSAGES  ON  A 19  NODE  ARPA  NETWORK  WITH  ONE  NODE  DISABLED 
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TABLE  5.14 

MESSAGE  THROUGHPUT  FOR  A 19  NODE  APPA  NETWORK  WITH  ONE  NODE  DISABLED 


NMSG 

BA 

FA  + 50% 

NR 

NR  + 50% 

SP 

SP  + 50% 

0.053 

1 

251(0) 

251(0) 

251(0) 

251(0) 

250(1) 

251(0) 

0.105 

2 

443(1) 

443(1) 

443(0) 

443(1) 

436(2) 

443(2) 

0.165 

3 

729(0) 

729(0) 

729(0) 

729(0) 

686(19) 

659(47) 

0.211 

4 

942(1) 

942(1) 

942(1) 

942(1) 

710(152) 

705(151) 

0.263 

5 

1216(0) 

1216(0) 

1216(0) 

1216(0) 

716(313) 

631(433) 

0.316 

6 

1418(4) 

1425(4) 

1418(4) 

1425(4) 

749(494) 

778(460) 

0.368 

7 

1575(50) 

1591(36) 

1595(35) 

1591(36) 

629(793) 

594(827) 

0.421 

8 

1530(203) 

1507(191) 

1530(203) 

1507091) 

626(905) 

740(820) 

0.526 

10 

1364(570) 

1411(553) 

1364(570) 

1411(553) 

708(1205) 

669(1150) 

TABLE  5.15 

MESSAGE  THROUGHPUT  FOR  A 19  NODE  ARPA  NETWORK  WITH  ONE  NODE  DISABLED 


N 

NMSG 

ARPA 

ARPA 

MODl 

ARPA 

MOD{  + 50% 

i 

0.053 

1 

250(0) 

250(0) 

250(0) 

0.105 

2 

443(0) 

443(0) 

443(0) 

0.165 

3 

730(0) 

730(0) 

730(0) 

0.211 

4 

942(0) 

942(0) 

944(0) 

0.263 

5 

! 2 1 5(0) 

1215(0) 

1212(0) 

0.316 

6 

1432(0) 

1430(1) 

1425(5) 

0.368 

7 

1490(123) 

1468(110) 

1427(137) 

0.421 

8 

1475(250) 

1472(247) 

1420(284) 

0.526 

10 

1277(724) 

1320(670) 

1345(679) 

SUMMARY  OF  19  NODE  ARPA  NETWORK 

The  results  of  the  19  node  ARPA  network  do  not  demonstrate  the  dramatic  effects  of  program  modifi- 
cation as  in  the  8 node  network.  The  8-node  network  was  a highly  connected  one,  in  addition  to  being  very 
symmetrical  in  its  topology.  As  a consequence  of  network  design,  it  may  be  stated  with  some  reservation  that 
the  high  connectivity  would  lend  itself  to  higher  throughput,  smaller  message  delays,  and  fewer  undelivered 
messages,  even  with  some  small  percentage  of  its  links  disabled. 
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rhe  19  node  ARPA  network  is  not  highly  connected  and  appears  vulnerable  to  any  disabling  of  nodes 
and  links.  Higher  connectivity  results  in  greater  overhead  cost,  hence  the  ideal  topology  of  N nodes  with 
h\N  • 1)  links  fully  duplexed  would  be  very  costly.  With  the  19  node  net  structuring,  the  ARPA  routine  works 
fairly  well  and  any  modifications  to  it  improve  it  only  slightly.  However,  the  ARPA  routine  takes  a lot  of 
updating  and  bookkeeping,  hence,  from  the  standpoint  of  minimum  computer  hardware  and  software  it  could 
not  be  viewed  as  optimum. 

A preliminary  conclusion  at  this  juncture  would  be  to  state  that  MODj  with  or  without  the  aging  aspect 
tends  to  help  3ny  routine  thus  far  studied.  It  is  further  envisioned  that  a more  simpler  algorithm  than  the 
ARPA  used  in  conjunction  with  MODj  and  the  aging  parameter  could  yield  effective  results  for  most  net 
topologies,  as  demonstrated  by  BL  + 50%.  Backward  learning  is  a rather  simplified  concept  as  well  as  easy  to 
implement  in  software.  If  we  allow  that  simplicity  in  software  implies  minimum  hardware,  then  we  could  assert 
that  some  routine  similar  in  implementation  to  BL  might  have  a significant  effect  on  the  ARPA  network. 

It  is  noteworthy  that  a value  of  50%  for  the  aging  factor,  proved  to  be  optimum  on  both  net  configurations, 
with  and  without  destruction.  Further,  this  value  of  50%  is  in  keeping  with  our  preliminary  predictions  in 
Chapter  III,  that  throughput  peaked  in  the  vicinity  of  a priority  threshold  setting  of  k - 50%. 

Chapter  VI  will  compare  some  numerical  computations  against  some  selected  simulation  data  of  this 
Chapter.  The  comparison  will  show  how  the  predicted  thioughput  factor  agrees  with  the  simulated  data  when 
using  the  optimal  value  of  priority  threshold  k. 
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VI.  THEORETICAL  MODEL 

Communication  networks  do  not  lend  themselves  to  tractable  mathematical  models,  particularly  when 
the  networks  are  inter-connected.  The  store-and-forward  procedures  of  computer  communication  further  com- 
plicates the  analysis.  When  adaptive  routing  vs.  fixed  routing  is  incorporated,  the  mathematical  modeling  is 
virtually  impossible.  The  problem  becomes  obfuscated  even  further  when  finite  buffers  (queues)  are  incorporated, 
along  with  finite  limits  on  channel  capacities.  There  is  in  fact  no  known  mathematical  analysis  to  describe  a 
physical  computer  communication  network  under  all  loading  conditions. 

The  only  approach  to  this  seemingly  insoluble  dilemma  is  by  some  appropriate  mathematical  approximations. 
The  technique  is  generally  one  of  analyzing  a single  node  for  the  parameters  of  interest  and  extending  the 
analysis  over  the  entire  set  of  nodes  comprising  the  network.  The  theoretical  model  then  becomes  a close 
approximation  to  actual  operating  conditions  within  the  network. 

6.1  AVERAGE  MESSAGE  DELAY 

it  has  been  stated  in  prior  chapters  that  one  of  the  measuremeni  parameters  of  interest  is  the  average 
message  delay.  The  average  message  delay  gives  an  excellent  over-view  on  network  performance  when  viewed 
by  the  user.  High  message  delays  would  discourage  users  who  are  attempting  to  use  the  system  and  would 
tend  to  reduce  subscribers  in  such  a system.  For  a fixed  investment  in  initial  overhead  for  setting  up  a com- 
puter network,  fewer  subscribers  would  tend  to  increase  the  cost  per  subscriber  to  the  point  of  further  reducing 
network  usage.  Hence,  in  a real  system,  the  average  message  delay  puys  a very  important  role  and  could 
directly  effect  the  network  utility,  or  lack  of  it. 

The  average  message  delay  has  to  be  calculated  on  the  basis  of  whether  or  not  there  is  some  priority 
classification  and  on  the  average  number  of  nodes  visited. 

6.1.1  Average  message  delay  with  no  priorities  in  the  system 

The  delay  experienced  by  a message  has  been  developed  from  queueing  theory  and  can  be  developed  by 
two  different  theoretical  approaches:  one  through  differential  difference  equations  due  to  Erlang,  and  the 
other  through  integral  equations  studied  by  (Lindley  52].  There  is  an  excellent  resume  of  queueing  theory  by 
Saaty  (SA57]  giving  most  of  the  formulations  generally  used  in  queueing  analysis. 

In  a single  server  queue  with  Poisson  arrivals  and  exponential  servicing  the  average  delay  in  the  steady 
state  is  given  by: 


W = 


f tP(< t) 


— , P < I 


(6.1) 


\ 
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where 


n<T)  = £ p„ H„(r)=  • (6.2) 

n=o 

and 


n -1  _ - 

*n(T)  = l-£  W///  ^T. 

i«o  L J 


(6.3) 


When  the  service  time  is  arbitrary  and  the  input  process  is  still  Poisson,  the  steady  state  delay  becomes: 


IV  = — £ — (l+(i/i)2]  (6.4) 

-/i(l  P) 

Where  s is  the  standard  deviation  of  the  service  time,  hence  for  constant  servicing,  s-o  arid  the  average 
delay  becomes: 


W = —?■— 


(6.5) 


It  should  be  noted  here  that  the  assumption  on  service  discipline  was  FCFS,  however,  when  the  queue 
discipline  is  changed,  the  distribution  of  waiting  times  change  but  not  the  average  waiting  time.  One  can 
further  note  from  equations  (6.1)  and  (h.5)  that  the  average  waiting  time  is  directly  effected  by  the  service 
distribution.  The  average  delay  for  constant  servicing  is  1/2  that  for  exponential  servicing  for  the  same  utiliza- 
tion factor  p. 

The  average  delay  given  by  the  above  formulations  are  associated  with  the  delay  at  a single-server  or  node. 
If  we  knew  the  average  number  Ar  of  nodes  visited  by  a given  message,  then  the  total  average  delay  in  the  system 
of  an  arbitrary  message  would  be: 


DT  = N W.  (6.6) 

liquation  (6.6)  has  been  further  simplified  where  it  is  assumed  that  the  delay  is  the  same  at  all  nodes 
visited  and  are  independent  random  variables  identically  distributed. 

6.1.2  Average  delay  for  a two  class  priority  system 

Priority  queueing  assumes  that  the  class  of  messages  arriving  at  a given  service  facility,  need  special 
handling  either  within  the  queue  or  with  their  service  requirements.  The  priority  classification  can  be  established 

with  some  parameter  p{\,2 k ),  where  1 denotes  the  class  with  highest  priority  and  k the  lowest.  The 

service  discipline  will  be  non-preemptive  and  FCFS  within  the  priority  groups. 
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Cobham  [COB54]  formulated  a generalized  equation  for  finding  the  expected  waiting  time  for  messages 
of  each  priority  level  in  a single-server  (unsaturated)  waiting  line  system.  His  equation  takes  the  form: 

©o 

*P  = TX  jf  jXI-o,,)  , (6.7) 


where 


N 

Pi  = \IPi  and  X =£x,  , (6.8) 

1 

k 

°p  =YlPi  < 1 °o-°  ' 

1 

N 

Bin  = (1A  • (6.9) 

1 


Bj(t)  is  the  cumulative  service  time  distribution  function  for  the  i^1  priority.  Thus  for  a two  clas  priority  sys- 
tem (/)=!, 2),  the  average  message  delay  for  the  priority  message  becomes: 


W 


1 


(6.10) 


where 

oo 

lt'0  - -^-X  f rdB(t)  2 average  time  to  complete  existing  service  . 

~ •'o 

For  the  non-priority  (/?=2)  messages,  we  obtain 


»;2  * ltV<1*M1'Vp2>  . (6.11) 

When  priority  assignments  are  established  within  the  system,  such  that 

P(k)  = 1\(Dt  - km)  . (6.12) 

where  p(k),  k.  and  m arc  as  defined  in  Chapter  111.  then 

X/>(A)  = X,  . 

and  (6.13) 

X(M*))  = x2 
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The  unconditional  de’ay  becomes: 


W = p{k)Wl  + (l-zKV))^  (6.14) 


where 


1 *Pj*P2  = 1-P- 


Cole’s  formulation  {(0  71))  has  essentially  the  same  form  and  is  given  by: 


W = MFCFS) 


where 


(6.15) 


x = a .Xj  + (i-a)  *2 


and 


W(FCFS)  lVo/(l-p). 

For  the  case  where  servicing  is  the  same  for  both  classes,  i.e.,  ^en  e<luat*ons  (6.14)  and  (6.15) 

reduce  to 


W = W(f'CFS)  f (6.16) 

since  x = x ^ , when  Pj=P2  *n  ecluahon  (6.15). 

The  average  message  delay  Df  cannot  be  expressed  like  that  for  the  single  class  system  given  by  equation 
(6.6),  but  instead  must  be  given  by 

Dt  = n Wx  + m W2  . (6.17) 

where  ft  = average  number  of  nodes  visited  as  a priority  message 
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and 


m - average  number  of  nodes  visited  as  a non-priority  message. 

Equation  (6.17)  assumes  messages  are  independent  random  variables,  identically  distributed  for  a given 
message  class.  As  with  equation  (6.6)  the  average  number  of  nodes  visited  is  also  a random  variable,  wlnse 
distribution  has  to  be  determined.  The  determination  of  the  distribution  function  for  the  number  of  nodes 
visited  by  a message  within  a communication  network  has  not  yet  been  solved. 

6.2  EFFECT  OF  A HIGH  UTILIZATION  FACTOR  p 

The  discussion  and  analysis  up  to  this  point  has  dealt  with  steady  state  formulations  (p  < 1).  When  the 
utilization  factor  p-*l  or  p > 1,  then  the  system  has  no  steady  state  solution  and  the  average  waiting  time 
approaches  infinity.  P.M.  Morse  [M058]  has  studied  the  transient  case  where  p > 1 and  develops  some  formu- 
lations for  the  time-dependent  probabilities.  Although  there  exists  a mathematical  model  for  over  saturated 
systems  (p  > 1),  these  formulations  are  more  esoteric  than  of  practical  concern. 

In  a real  physical  communication  network  the  conditions  when  p > 1 are  intolerable  from  the  standpoint 
of  prohibitive  queueing  delays.  Although  a real  system  can  become  over  saturated  this  is  generally  a short  term 
phenomena  since  users  will  tend  to  subside  when  networks  arc  over-loaded  resulting  in  the  utilization  becoming 
unsaturated  again,  i.c.,  p < 1.  Thus  in  the  final  analysis  the  transient  case  is  just  that,  a condition  that  will  not 
persist  very  long  and  steady  state  will  be  the  rule,  rather  than  the  exception. 

6.3  EFFECT  OF  PRIORITY  ASSIGNMENT 

The  justification  for  priority  assignment  generally  stems  Irom  desiring  to  minimize  the  average  message 
delay  and  maximize  the  system  throughout.  The  majority  of  priority  disciplines  give  high  priority  to  those 
classes  with  low  mean  servicing  time.  The  proof  that  such  a rule  leads  to  optimality  requires  the  a priori 
arrival  and  service  distributions  of  each  class  are  known  and  have  stationary  statistics. 

The  technique  employed  in  this  research  suggests  that  the  overall  performance  of  the  system  may  be 
improved  if  a message,  after  reaching  some  queueing  delay,  is  given  in  a Ron-preemptive  priority.  The  idea  was 
to  minimize  the  queueing  delay  by  an  optimal  selection  of  priority  threshold  p{k).  Where  p(k)  has  been  defined 
as 


p{k)-Pf{km  < / < m),  0 < k < 1 


(6.18) 


where 


m - maximum  allowable  message  delay. 
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However,  it  appears  at  first  glance  from  equation  (6.14)  that  the  average  message  delay 


W 


,P<  1 


(6.19) 


is  invariant  to  the  selection  of  p(k),  equation  (6.19)  is  equation  (6.14)  repeated  for  convenience.  Since  we 
defined  the  priority  rate  as 


\ = pm . 

and  when  the  servicing  ‘or  both  priority  classes  are  equal 


(6.20) 


Px  = P2  = P ■ 


Then  equati-.i  (6.19)  reduces  to 


IV 

IV  = — = I \\FCFS) 

1-P 

since 

X1 

p.  = — = pm/p  = pm 

p 

There  appears  to  be  some  contradiction  with  (6.22)  and  the  results  of  the  simulation  given  in  Chapters  IV 
and  V tha.  show  average  message  delay  is  reduced  whrtn  priority  assignment  are  initia  ed.  The  simulation  gives 
the  average  message  delay  over  the  network  as  opposed  to  an  arbitrary  node  given  by  (6.17)  and  (6,19).  The 
missing  parameter  that  makes  up  for  this  apparent  discrepancy,  is  the  average  number  of  nodes  visited.  Note 
in  equation  *6.17)  the  random  variables  n and  m repiescnting  the  number  of  nodes  visited  by  a priority  and 
non  priority  message  respectively.  These  two  random  viriablcs  are  the  quant  ies  that  can  not  be  readily  assessed 
by  analytic  rrrans  end  are  directly  effected  by  the  updating  routines  that  in  aim  are  effected  by  priority 
assignment. 

The  priority  assignmu..:  helps  to  relieve  conge  . ion  at  the  nodes,  hence  reducing  over  all  queueing  delry. 
The  reduction  in  queueing  delay  at  previously  congested  nodes  change  the  routing  of  messages  and  consequently 
the  number  of  nodes  visited  The  final  effect  results  in  a smaller  vaiue  of  average  message  delay  Df  when  viewed 
over  the  entire  network.  The  adaptive  routing  process  negocs  mathematical  modeling,  therefore  one  can  only 
obtain  a heuristic  analysis. 
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6.3.1  Balking  and  reneging 

Balking  is  the  condition  where  an  arriving  message  may  not  join  the  queue  because  of  the  length  of  the 
existing  queue.  In  this  study  the  queues  are  of  finite  size,  hence  balking  is  attributed  to  those  messages  finding 
the  queue  full.  The  utilization  factor  directly  effects  balking  by  fiiling  the  queues  faster  than  servicing  can  dis- 
patch the  messages.  For  a fin  te  queue  of  size  N,  with  Poisson  input  and  exponential  servicing,  P.M.  Morse  gives 
the  following  probability  of  finding  n messages  in  the  queue 


Pn  - 


1 _ 


1 "P 


‘A- 

,A'+1 


p\p<  1 


(6.23) 


where 


A = queue  size  - 

The  mean  number  in  the  system  (in  queue  and  in  service)  is 

') 

(p«D  (6.24) 

(P  -**  1 ) 

(P»D  . 

The  piohability  of  balking  becomes: 

PS  (6.25) 

The  number  of  ! :essap,es  finding  the  queue  full  is  heavily  dependent  on  p.  If  for  a fixed  service  rate  p, 
the  message  arrival  rate  X increases,  then  the  fraction  of  messages  bulking  given  by  PN  increases  rapidly, 
approaching  M 1/p)  when  p is  much  larger  than  unity.  Thus  in  an  over  saturated  system  (p  » I)  only  (1/p) 
of  the  messages  get  into  the  system.  Conversely  when  p <X  1,  the  system  is  being  under-utilized  and  the 
probability  of  balking  approaches  zero.  Balking  can  also  be  viewed  from  the  probability  that  the  queue  is 
empty  given  by 


1-P 


1-p 


S H 1 


A' 


= z < 


n= 0 


KA’-H  )p'V+  A) 

(1-pK1-pV4 


p 4-p 


1 


A'(A'+2Kp-l ) 

AW  1/p) 


/ 

o 


1-p 


I-p 


:V+1 


(6.26) 
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From  (6.26)  one  can  see  that  for  p « 1 PQ  -*  1,  hence  no  balking  and  for  p » 1,  PQ  = l/p^  and 
balking  is  highly  probable.  The  number  of  units  balking  cannot  be  readily  calculated. 

A renege  is  when  a message  due  to  its  aging  nny  expire  from  the  system.  The  probability  that  a message 
will  renege  is  given  by 


Pf(DT>  MDTmax)  = reneging, 


(6.27) 


where 

DT  - n h'j  + ,v  W2  (sec  6.17) , 


and 


MDTmax  is  the  maximum  allowable  message  delay. 

When  p -*>  1,  the  equations  for  average  message  delay  given  by  (6.10)  and  (6.1 1 , Dive  very  large  values  for 
IVj  and  h'2  respectively  (these  equations  arc  not  defined  for  p ^ 1 ).  The  probability  that  a message  will  expire, 
hence  renege,  becomes  very  probable  as  p -*  1.  The  numerical  number  of  messages  that  have  expired  can  be 
approximated  from  the  total  messages  generaitJ  per  unit  time,  times  > MDTmax).  When  the  input  is 

Poisson  the  number  of  messages  reneging  per  unit  time  becomes: 

A’ 

X = ^ ' (6  28) 

1=1 

Then  yn  will  be  the  expiring  messages  per  unit  time  in  a A'  node  network. 

6.3.2  Throughput  factor  $(t,k) 

The  utilization  factor  also  directly  influences  the  throughput  factor  0(r,/v)  as  defined  in  section  (3.5).  The 
defining  equation  lor  the  throughput  factor  given  by  equation  (3.43)  show  the  implicit  dependence  of  ' )on  p. 
Although  the  equations  arc  for  values  . . < I,  if  can  be  readily  shown  that  the  throughput  factor  decreases 
with  incre  ing  p.  1 rom  the  basic  definition  of  the  throughput  factor  it  is  clear  that  an  increasing  p,  increases 
the  percent  of  undelivered  messages  v(t.k),  hence  decreases  <H*,A:).  The  influence  of  p on  is  further 
demonstrated  in  the  next  section. 

6.4  COMPARISON  OF  PRI  DICTHI)  TO  SIMULATION  R!  SUITS 

The  computer  runs  of  Chapters  (IV)  and  (V)  showed  some  tendencies  to  hear  out  analytic  predictions. 
Numerical  calculations  of  average  queueing  delay,  number  of  messages  balking  and  reneging  caunot  be  readily 
assessed  due  to  the  inter-connectivity  of  computer  networks,  and  to  alternate  (adaptive)  routing  procedures. 


76 


There  appears  to  be  an  area  where  the  predicted  values  were  to  within  fractions  of  a percentage  to  the 
simulation  results  for  the  throughput  factor  <p(t,k).  Table  (6.1)  shows  the  computed  $(t,k)  versus  the  simulated. 
The  agreement  was  very  good  from  almost  an  unmeasurable  difference  to  within  1.1%  for  the  worst  case.  The 
comparison  was  made  with  the  cxict  values  of  p and  priority  k for  both  cases.  The  curves  of  Fig.  (6.1)  are  a 
plot  of  the  tabulated  values  of  Table  (6.1). 


TABLE  6.1 

COMPUTED  VS.  SIMULATED  THROUGHPUT  FACTOR  0(U) 


CASE 

NUMBER 

UTILIZATION 
FACTOR  p 

PRIORITY 
k FACTOR 

% UNDEL. 
MESSAGES 
v(t,k) 

THROUGHPUT 

FACTOR 

<P(t.k) 

% DIFFERENCE 
BETWEEN 

COMPUTED  & SIMULA. 

SIMULATED 

1 

.158 

.50 

0.0 

1.0 

No  Signitlca*  ■ 

COMPUTED 

i 

.158 

.50 

0.26  x 10‘6 

.9999 

Difference 

SIMULATED 

.316 

50 

0.0 

1.0 

No  Signn.cant 

COMPUTED 

•> 

.316 

.50 

0.42  x I0‘6 

.9999 

Difference 

SIMULATED 

3 

.474 

.50 

.052 

.948 

.0007 

COMPUTED 

3 

.474 

.50 

.0513  j 

.9487 

SIMULATED 

4 

.632 

i 

.50 

i 

.079 

.921 

.011 

COMPUTED 

I : 

4 

I 

1 .632 

1 

.50 

.068 

i 

.932 

j 

j SIMULATED 

5 

.790 

.50 

1 

j 

.203 

.797 

.008 

1 COMPUTED  ! 

5 

.790 

.50 

.195 

.805  | 

The  simulated  data  comes  *r<mi  the  ARP  A run  using  Backward  Learning  plus  i OC  aging  (BL+5u'/7),  and  these 
are  tabulated  in  Table  v7  of  C hapter  V.  The  computed  data  refers  to  theoretical  calculations  obtained  from  Eqs 
(3.42).  (3.43),  and  (3.44)  ami  is  associated  with  the  3-node  network  of  Chapter  III.  A judicious  choice  of  Aj'  and 
> •>'  m 1 q.  (3.44)  created  the  same  loading  effects  in  the  computed  (theoretical)  3-node  net  as  the  19  node  simulated 
net.  even  though  the  two  networks  had  different  message  volumes  for  the  same  value  of  p. 
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Figure  6.1  - Computed  Vs.  Simulated  Throughput  Factor  <b(t,k) 

It  ’s  also  instructive  that  the  percent  of  undelivered  messages  v{t,k)  is  analogous  to  the  percent  of  messages 
reneging  in  Jhe  system.  Thus  there  is  a way  to  get  a reasonable  prediction  of  the  percentage  of  messages  that 
will  be  undelivered  for  some  selected  values  of  p.  The  actual  predicted  throughput  \J/(r)  for  any  network  loading 
can  be  roughly  approximated  by  the  following  (where  the  input  is  Poisson). 

N 

WO  = \T  , (6.29) 

/=1 

where 


X.  = rate  per  node/unit  time  . 

T = total  time  increment  , 

/V  = total  number  of  nodes  in  network. 

In  an  actual  network  one  could  not  readily  find  the  various  values  of  and  these  would  have  to  be 
approximated  by  appropriate  assumptions  regarding  network  loading  and  average  servicing  rate. 
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VII.  SUMMARY  AND  CONCLUSIONS 


7.1  SUMMARY 

Several  techniques  for  adaptive  routing  algorithms  in  a specific  class  of  networks  have  been  examined.  In 
particular,  it  has  been  demonstrated  by  simulation  how  some  relatively  simple  add-ons  to  already  existing  adaptive 
routines  could  decrease  the  average  message  delay  and  increase  throughput  in  the  selected  networks.  A further 
objecbve  of  this  study  was  to  assess  the  effects  of  priority  assignment  to  messages  that  had  reached  some  speci- 
fied aging  u vshold  and  note  the  effects  of  such  priority  assignments  on  network  performance.  The  perfor- 
mance measures  considered  were  average  message  delay,  throughput,  and  number  of  messages  undelivered. 

A historical  review  of  communication  networks  was  given  in  Chapter  I.  Chapter  II  dealt  with  determinis- 
tic and  stochastic  routing  strategies  and  a brief  overview  of  how  they  operate.  It  was  asserted  in  Chapter  III 
that  there  was  a general  lack  of  knowledge  of  system  specifications  that  renders  extensive  analysis  of  store-and- 
forward  networks  to  be  essentially  intractable.  However,  a measurement  parameter  called  the  throughput  fac- 
tor dU.k)  was  introduced.  This  parameter  could  be  optimized  in  terms  of  the  selection  of  the  priority  thresh- 
old k The  closed  form  solution  for  4{i,k)  was  made  on  a 3-node  network  and  it  was  envisio^d  that  this  analy- 
sis could  be  extended  to  predict  performance  on  some  specific  interconnected  network.  The  problem  of  finite 
buffers,  network  interconnections,  and  adaptive  routing  makes  exact  analysis  intractable. 

Chapters  IV  and  V dealt  with  the  simulation  study  on  an  8 node  and  19-node  network  respec lively.  The 
simulation  was  to  test  the  performance  of  some  selected  algorithms  along  with  the  modifications  suggested  by 
this  study,  and  note  how  the  algorithms  faired  when  there  was  no  network  damage  and  Mien  with  fnk  and  node 
failures.  The  8-node  highly  connected  network  was  a preliminary  run  with  the  roufi.ig  algorithms  and  the  19- 
node  Akf  A net  was  an  attempt  to  verify  the  consistency  of  the  algorithms  on  a larger  net  topology. 

Finally.  Chapter  VI  was  a comparison  of  predicted  performance  with  actual  simulated  results.  This  com- 
parison of  the  analytic  with  the  simulated  can  be  only  viewed  as  cursory  due  to  the  limitations  of  the  mathe- 
matical modeling. 

7.3  CONCLUSIONS 

The  technique  of  using  last  node  visited  (LNV)  and  an  aging  parameter  for  priority  selection  can  be  viewed 
as  an  asynchronous,  stochastic  routing  algorithm.  When  the  network  is  lightly  loaded  and  there  are  no  link  or 
node  outages,  then  the  minimum  path  algorithms  would  work  rather  well.  When  the  network  gets  congested 
at  some  node  or  nodes  and  traffic  intensity  becomes  heavy,  the  adaptive  louting  algorithms  suggested  in  this 
paper  attempt  to  unblock  the  system  and  keep  average  message  delay  at  a minimum  and  throughput  at  a maxi- 
mum. Hence,  there  appeared  to  be  a significant  improvement  over  all  routing  algorithms  when  MOD^  and  50% 
aging  was  applied  to  them.  It  is  noteworthy  that  the  50%  thn:hold  value  was  the  optimum  setting  in  the  simu- 
lation and  corresponded  closely  to  the  predicted  value. 
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There  are  a couple  of  caveats  that  should  be  stressed  at  this  point.  First,  although  Chapter  VI  showed  a 
very  close  comparison  (within  1%)  of  the  predicted  throughput  factor  4(t,k)  with  that  of  the  simulated  data  of 
Chapter  V,  it  cannot  be  considered  conclusive.  The  predicted  </(/,&)  was  made  on  a net  topology  where  the 
average  path  length  did  not  correspond  to  that  of  the  simulated  networks.  However,  by  a judicious  choice  of 
input  messape  rates  \2'  and  ' it  was  possible  to  create  the  same  net  loading  effects  thus  giving  some  very 
close  comparisons  between  predicted  and  simulated  throughput  factor. 

The  second  area  of  caution  deals  with  the  appearance  that  backward  learning  (BL)  can  outperform  the 
ARPA  routine,  as  demonstrated  in  Chapter  V.  Both  ARPA  and  BL  inquire  about  nodes  and  links  via  nearest 
neighbor.  BL  does  this  inquiry  by  incoming  messages,  hence,  under  heavy  loading  it  gets  current  information. 
ARPA  works  better  at  lightly  loaded  nets  using  nearest  neighbor  information,  and  at  higher  traffic  intensities 
its  updating  is  slightly  delayed;  also,  it  adds  to  the  net  loading  with  its  inquiries.  However,  the  BL  routine 
cannot  adjust  effectively  to  nodes  or  'inks  being  restored  and  its  updating  techniques  generally  yield  inadequate 
or  even  misleading  table  delay  values.  The  tendency  to  “ping-pong”  is  very  probable  when  using  BL  techiiiques. 
The  ARPA  routine  is  both  synchronous  and  asynchronous  in  that  it  updates  every  Vi  second,  but  also  if  there 
have  been  link  or  node  changes. 

In  either  case,  it  appears  that  adding  the  asynchronous  technique  of  MODl  and/or  50%  aging,  all  the 
selected  algorithms  performed  better  and  in  some  cases  (Chapter  IV)  the  effects  were  dramatic. 

7.3  FUTURE  WORK 

It  appears  that  additional  simulation  techniques  are  desirable  to  further  evaluate  routing  algorithms  on 
store-and-forward  systems.  The  technique  of  asynchronous  updating  or  adaptation  should  be  further  explored, 
since  such  techniques  minimize  network  loading  by  minimizing  traffic  inquiries  about  the  state  of  the  system. 

The  problem  of  analytic  modeling  still  needs  to  be  addressed  to  obtain  a closer  approximation  of  predicted 
performance  to  simulation  results.  Mot  only  is  there  the  apparent  limitation  of  the  mathematics  but  also  the 
normal  lack  of  complete  knowledge  of  tne  system  specifications  on  which  to  use  such  analysis. 

The  area  of  an  adjustable  threshold  might  also  prove  enlightening.  The  technique  in  this  study  was  to 
select  an  optimal  value  of  k for  all  p < 1.  It  would  be  interesting  to  see  how  selecting  k based  on  the  actual 
message  mtensity  would  effect  average  message  delay  and  throughput. 
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APPENDIX  A 

QUEUEING  THEORY  AND  FORMULATIONS 


A.l  QUEUEING  THEORY 

A queue,  or  a waiting  line,  involves  items  arriving  at  some  service  facility  and  their  subsequent  wait  to  be 
served.  The  items  in  the  queue  may  be  people,  vehicles,  commodities,  or  almost  an>  thing  that  joins  a waiting 
line  for  the  purpose  of  being  served;  there  are  almost  unlimited  specifications  for  the  server. 

The  principal  elements  of  a queue  are  the  inputs,  waiting  line,  and  the  service  facility;  thus  the  statistics 
of  a queue  are  associated  with  the  fluctuations  of  these  three  elements.  The  ratio  of  the  rate  of  input  (X)  and 
the  departure  rate  (/i)  has  been  generally  defined  as  utilization  factoi  p = X/p , and  it’s  this  factor  that  directly 
affects  the  size  (length,  total  number,  etc.)  of  a given  queue.  Optimizing  any  system  of  queues  depends  directly 
on  minimizing  congestion  at  the  node  points  where  queues  develop. 

In  a sture-and-forward  communication  network  the  subject  of  queueing  must  be  addressed  to  handle  con- 
gestion at  the  nodes  within  the  network.  Messages  in  communication  nets  are  sent  from  node  i to  some  node 
/ (7  tj)  via  some  intermediate  nodes  k(k  = 0,  1 2V  . . ,,  Ar)  and  it’s  the  congestion  at  the  intermediate  nodes  k, 
where  queueing  occurs  and  adds  to  the  average  message  delay.  The  source  nodes  /,  destination  nodes  /,  and 
intermediate  nodes  k are  all  random  variables  generally  from  a Poisson  distribution.  The  problem  of  adequate 
waiting  room  within  the  queue  is  important,  especially  when  the  queue  size  is  finite. 

The  queueing  problem  for  nets  is  usually  concerned  with  knowing  the  input  distribution,  the  queue  dis- 
cipline (e.g.,  random,  ordered,  or  priority  selection  for  servicing),  and  the  service-time  distribution  to  determine 
the  desired  measures  o»  effectiveness.  For  problems  dealing  with  finite  queue  length,  one  must  also  detemiine 
the  probability  of  balking  (not  joining  the  queue).  In  addition,  due  to  time  limitations  imposed  on  the  mes- 
sage's stay  in  the  system,  one  must  determine  the  probability  of  abandoning  the  queue  (reneging)  after  joining 
it  due  to  exceeding  the  imposed  aging  limit. 

The  literature  on  queues  is  very  extensive  and  Saaty  |SAo  1 1 lists  over  900  references  and  develops  many 
other  specifications  for  describing  various  forms  of  the  qucuunv.  process;  Svski  [SY 60]  has  over  400  references 
on  queueing  and  congestion  theory-.  Morse  [M058]  and  Cox  (COX61 J are  other  texts  useful  for  the  engineering 
development,  and  Jaiswal  [JA68]  treats  exclusively  the  subject  of  priority  queues.  The  above  mentioned  ref- 
erences are  more  than  adequate  to  find  a very  diversified  treatment  of  the  queueing  process. 

The  queueing  model  generally  used  for  a network  is  the  single  exponential  channel.  Thus  one  can  get 
the  characteristic  behavior  of  a large  class  of  more  complicated  systems  when  the  following  are  assumed. 

1.  The  interarrival  times  are  Poisson-distributed  with  an  average  arrival  rate  of  X messages  per  unit 


time. 
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2.  The  message  lengths  are  exponentially  distributed  with  an  average  length  of  1 /p  bits  per  message. 

3.  Queue  discipline  is  first-come  first-served  with  infir:te  stoiage  capacity. 

4.  There  is  one  channel  available  for  transmission. 

With  the  above  assumptions,  and  with  the  abundance  of  available  literature  on  the  single  exponential 
channel  ([M058]and  [SA61],  the  following  steady  state  statistics  can  be  readily  found: 


a. 

b. 

c. 

d. 

e. 


Pf  [n  messages  in  the  system]  = (1  - p)p”,  p < 1 . 

t\n)  = expected  value  of  number  of  messages  in  the  system  = 

T - average  time  message  spends  in  the  system  = • 


Pr  [more  than  n me€P^c 


in  the  system]  = pAI+  * 


Pr  [total  time  message  spent  in  system  > f]  = . 


P 

1 -p 


(A.l) 


The  above  quantities  are  derived  with  the  aid  of  the  birth-death  process  and  resulting  forward  and  back- 
ward equations,  and  assuming  the  existence  of  a limiting  distribution  for  PfJ(t),  such  that 


lim  Pn  (t)  = Pn  , (A.2) 

[— ► OO 


and  further  that  the  limiting  distribution  is  such  that 


(A.3) 


The  above  mentioned  process  is  Markovian  and  the  techniques  for  solving  them  are  well  covered  in  the 
literature  [Feller  57  and  Kendall  51];  the  use  of  Little’s  [LI61J  equation  allows  for  a quick  calculation  of  the 
expected  time  a message  spends  in  a queueing  system  (in  a conservative  system).  Little  proved  that  in  a con- 
servative system  (i.e.,  messages  are  neither  created  or  destroyed)  the  following  is  always  true: 


n X 7 \ 


(A.4) 


where 

n ~ expected  number  of  messages  in  a queueing  system 
T expected  time  they  spend  in  the  system 

and 


X - average  arrival  rate 
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Thus  from  Eq.  (A.l),  item  b. 


t\n)  = n - fj  td\  = - — , p < 1 
" 1-p 

and  from  the  definition  of  utilization  factor  p 


(A.5) 


p = X/pi  % 


(A.6) 


we  can  find  the  average  time  T from  Little’s  equation 


T='nj\  =-2  -L-=J — 

^ (1-P)  /i  (1-P) 


(A.7) 


which  is  given  by  item  c.  The  foimulations  just  presented  were  for  a Markovian  M/M/1  system  with  infinite 
queueing  room,  where  both  arrival  and  service  are  Markovian. 


A. 2 IMBEDDED  MARKOV  ( MAIN 

The  system  under  consideration  in  this  study  is  a M/G/l  system.  This  type  of  system  employs  a general 
service  time  distribution  and  the  memoryless  property  no  longer  exists;  hence,  the  employment  of  the  birth- 
death  process  is  no  longer  valid.  The  notion  of  an  imbedded  Markov  chain  is  due  to  D.  G.  Kendall  [KE51&53] 
and  was  introduced  by  him  so  that  non  Markovian  processes  could  be  studied  by  extracting  a set  of  regenera- 
tion points  for  winch  the  Markov  property  holds.  In  a M/G/I  system,  the  set  of  departure  instants  from  ser- 
vice is  an  extremely  convenient  set  of  regeneration  points.  Thus,  if  we  specify  the  number  of  me^-^.,0  left 
behind  (in  the  queue)  by  a departing  message  (one  already  served),  we  can  calculate  the  same  quantity  at  some 
point  in  the  future  given  only  additional  inputs  to  the  system.  The  concept  or  technique  of  the  imbedded 
Markov  chain  will  allow  analysis  for  our  priority  queueing  system. 

A. 3 POLLAGZEK-KHINTCMINH  FORMULA 

The  average  number  of  queueing  customers  (those  customers  waiting  in  the  queue)  left  behind  by  a 
departing  customer  is  given  by  the  Pollaczck-Khintchinc  formula 

p 2 + X 2 var  (/) 

i:U\  = P+  (A.8) 

2(1  - p) 

Once  we  know  the  variance  of  the  service  time  / from  its  given  distribution,  the  average  number  of 
messages  in  the  queue  can  be  detcimined.  We  must  note  that  the  above  average  '^as  been  taken  over  instants 
just  following  departures  and  thus  does  not  include  what’s  presently  beuig  serviced. 
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To  obtain  the  average  waiting  time  we  can  follow  the  argument  presented  by  Saaty  [SA61]  who  argued 
as  follows:  If  we  write  /f[ JV]  for  the  average  waiting  time  in  the  queue  (not  including  service),  X [£* [ H7]  + 1 In] 
is  the  expected  number  of  arrivals  during  the  total  waiting  plus  service  of  one  message;  i.e.,  its  stay  in  the  sys- 
tem. But  this  must  be  just  the  number  in  the  system  immediately  after  its  departure;  namely  E[q].  Hence 


WqBE[W}  = 


p2  + X 2 var(t) 
2X  (1  - p)~ 


(A.9) 


We  can  obtain  equation  (A.9)  by  direct  application  of  Little’s  equation  and  substituting  X/p  for  p,  then 


E[q  ] _ I p + Xp  i far {t) 
X p 2p  (1  - p) 


(A. 10) 


where 

n = E[q]  = Lq  , 
and 

Wq  = T - 

Both  equatior  (A.9)  and  (A.  10)  have  the  same  reduced  form  and  state  that  the  average  total  time  spent 
in  the  system  is  the  average  time  spent  in  the  service  plus  the  average  time  spent  in  the  queue.  There  are  num- 
erous other  statistical  parameters  that  can  be  assessed  in  a queueing  problem  such  as  length  of  a busy  period, 
unfinished  work,  probability  of  being  idle,  etc.  The  list  can  get  quite  lengthy,  but  for  the  purpose  of  thic  .tudy 
we  have  limited  our  attention  to  what  has  been  previously  stated.  Our  main  concern  is  minimum  average  mes- 
sage delay  and  maximum  throughput  with  minimum  hardware  and  software  investment.  We  should  further  state 
that  this  minimum  average  message  delay  should  be  realized  from  a real-time  system  and  thus  some  further  con- 
straints must  be  employed  when  analytical  modeling  is  employed.  The  problem  of  priority  assignment  to  mes- 
sages in  the  system  must  be  dealt  with,  along  with  some  type  of  parameterization  to  optimize  the  queueing 
assignment.  The  aggregate  collection  of  nodes,  links,  buffer  size,  queueing  (with  and  without  priorities)  and 
network  topology  must  somehow  be  systematically  developed  and  assessed  to  get  maximum  utilization  out  of  a 
communication  network  with  minimum  overhead.  The  inteiplay  of  the  network  components  and  lack  of  some 
systematic  analytical  analysis  makes  for  a very  formidable  problem. 

A.4  WAITING  TIML  DISTRIBUTION 

The  waiting  time  distribution  for  the  non-priority,  single-server  queue  will  be  used  to  establish  the  priority 
assignment  for  all  messages;  i.e.,  prio*-*ics  arc  established  endogenously.  Messages  arc  given  a priority  if  their 
wait  in  the  system  (queue)  is  longer  than  s'»me  preselected  time  T,  and  (his  priority  threshold  T is  the  same  for 
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all  messages.  Thus,  the  probability  of  jwait  in  Q > T establishes  priority,  and 


Pr  (wait  in  Q>  t\  = 1 - J dW(t)  = priority  assignment  , 


dWU)  = the  Stieltjes  integral  , 


JV(r>  = waiting  time  distribution  for  non-priority  messages  . 


The  waiting-time  distribution  W(f)  for  the  non-priority,  single-server  queue  can  be  obtained  from  the 
Pollaczek-Khintchine  formula 


1 - X/s  (1  - m) 


/3 (s)  - / e ~st  dB(t)  - the  Laplace-Stieltjes  transform  of  the  service-time  distribution, 

O 


B{t)  = service-time  distribution 

liquation  (A.  12)  has  been  derived  by  many  authors  (e.g.,  R.  C».  Miller  [MI60] , L.  Takacs  [TA55],  and 
Riordan  (R162 ) as  being  the  steady  state  solution  for  a single  class  queue  with  Poisson  arrivals  (\)  and  general 
service  distribution  B{t).  Thus  for  a constant  service-time  distribution 


_{!,/>  1 hi 
1/ju  , 


the  derivative  of  B(t)  becomes  a Dirac  delta  function  of  the  form 


dB{t)  = Ht-  I lii)dt. 


OO 

s)=  f e~*r  8(  \ -•  Miiyit  = e~s/» 


I - X/s|  I - <fs/"| 


(A.i5) 
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The  inverse  of  W{s)  as  shown  by  Riordan  [RI]  and  to  have  been  cal  Jated  by  A.  K.  Erland  in  1909  (Ref. 
E.  Brockmeyer  et  al,  1948)  is  given  by 


m - 0 - n) 


[Mr] 

E 


n=0 


(*1  f 


n ! 


(A- 16) 


where 

[pt]  = least  whole  integer  , 


and 

M=  1 

in  this  investigation. 

A. 5 C7  LCULATION  OF  THROUGHPUT  FACTOR 
The  throughput  factor  is  expressed  by 


0(a)  = ! ~ v(t,k)  , 


(A. 17) 


where 

v{t,k)  = l((X2'/2)  r(t,k)  + X/5(a))]/(X272  + X;')  . 

The  solution  of  v(t,k)  is  given  in  terms  of  r(t,k)  and  S{t,k)-  From  Fig.  (3.2)  we  defined  r(t,k)  as  the  prob- 
ability that  a message  with  source  S and  destination  D is  undelivered  and  S(t,k)  is  the  probability  that  a message 
with  source  M and  destination  D is  undelivered.  Then  r{t,k)  is  the  sum  of  the  following  mutually  exclusive 
independent  events: 


r(t.k)  = g(m)  + (1  - X2r)/>(/l)  + p(k)P(B)  + u(k )P{c)  , (A.  18) 


and 

S{trk)  = P(A)  , 

where 

gim)  - IX  {t  > m)  - probability  that  a message  with  source  S and  destination  D,  expires  before  it 
reaches  node  M. 

(’  - X^'jPiA  ) the  joint  probability  that  a message  finds  the  queue  empty  at  S'  and  expires  at  D. 
p(k)PUi)  ^ the  joint  probability  that  a message  belongs  to  a priority  class  and  expires  at  IX 
u(k)P{C)  = the  joint  probability  that  a non-priority  message  wi’h  time  on  it  expires  at  IX 


86 


We  should  note  that  when  dealing  with  a single  class  system,  H^t)  = //2(f)  for  p2  = 0 and  pi  = 0 respec- 
tively. The  distribution  functions  H^t)  and  //2(0  are  given  by  Eqs.  (3.1)  and  (3.3),  denoting  priority  and  non- 
priority  distributions  respectively. 

Unless  otherwise  noted,  it  has  been  assumed  that  the  average  service  rates  are  unity  and  given  by 

p = Pj  = p2  = 1 (A.  19) 


and  hence 


p " X,  pj  = Xj , and  p2  = X2 


The  three  probabilities  P(A ),  P{B ),  and  P(C)  have  been  defined  as  the  conditional  probability  of  some 
event  given  the  event  that  the  observed  value  of  X is  equal  to  x,  denoted  in  symbols  by  P{Y\X  = x).  Where 
the  event  Y and  the  random  variable  X are  both  defined  on  the  same  probability  space  [PARZEN  62]. 

From  a knowledge  of  P(Y\X  = x)  one  may  obtain  P{Y)  by  the  following  Stieljes  integral  formulas: 


OO 

/ P(Y\X  ^ x)dFx(x) 


P(D=  \ f P{Y\X  = xVx  (x)dx 

1 mOO 

T,  P(A \X  = x)px(x). 


(A. 20) 


over  ail  x such  that  px  (jc)  > 0 


in  which  the  last  two  equations  hold  if  X is  respectively  continuous  or  discrete.  The  solution  to  P(A)  and  P(B) 
are  as  given  in  the  text;  however,  P{C)  needs  further  clarification.  P(C)  is  the  probability  that  a message  that 
had  to  w'ait  at  node  S (Fig.  3.2)  is  undelivered  at  D. 


P(C)  = lY{r  > m)  - 


(A.21) 


where 

T - 7’(  + 7\  = total  time  spent  in  system 
Tj  = time  ^pent  at  node  S ) 
jT2  = time  spent  at  node  D \ 


random  variables  assumed  to  be  independent  . 
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P3(t)  = Pr(T  >m\Tx  = /) 

= Pr{Tx  + T2  > m | Tx  = 0 
= />cr.  > m - t | T,  = t) 


and  due  to  the  independence  of  7’1  and  T2 


P3(t)  = Pr(T2  > m - /), 


also 


!0,r  > km 

(1  - X2>'<i'X2,)  ' I [*2'(l-e(,‘V)*m))>  otheiwise 

0,f  < 0. 


Then  by  application  of  the  distribution  function  for  non-priority  messages  given  by  7/,(r)  in  Hq. 
obtain: 

P3  (f)  = P- P|<r«"i-r)  + — f ert  JJcu  - r)(r-  a.Vrir  - a) dr,  for  p22  p.  , 

i P - P,  2ir  *&,  v 2 1 1 

where 

<*  = (P  - Pj)  0 “ P)/P, 

= o Vn)2  • 

02=0  +v/n)2- 

Then 

km 

m=  f p^tVj  m 

o 

a - V) 

(1  - p*u'X2^*m) 

where 

Hr)  ^ cr' - r)  (r  a,)  Mr-  a). 


r™  e-<lV»  ( p2~_ple-a(-n-()t  |-pJ*  /v'p?-ld/ 
P - P,  2ir  ^ 


(A. 22) 


(3.3)  we 


(A.23) 


(A. 24) 


Let 


A,  =(1  - X2')  / (1  - e-(|-V>*m)  , 


(A.25) 


then 


no  - a,  (-— -i  )*■ 


^ ^_1  )e~Ctm  km 

P _ Pj  J*  e*(i-X2  - 


0\.26a) 


Arm  , OU 

A.  (-^P)  f C(I‘X2  *'  f f(r)drdt 
2v  Jo  J 


for  p1  > p. 


and 


Arm 


A(0  = A.(---)  / .•(1'X2,)'  / f(r)drdl  (A.26b) 

o “l 

for  p2  < Pj 


We  define  /j  and  /2  as  follows: 


km 


/,  = A,  ( «■“"  f c < 1 V“> ' <// 

‘ 1 P-P,  J 


(A* 27) 


P-P, 


Arm 


/2=A1  (J-^)  j e-O-Nt'w  f f(r)drd, 


(A.28) 


Arm 


'1 


= *!  J >'  e'rimml\/{a2-~F)(r-tti)l  \r{r-Q))drdt  , 

and  reversing  the  order  of  integration  /2  becomes: 

1-P 


2 = *i<  -- 


, “2  f''ny  (a,  - r\r  - a.)  , 


(A-29) 
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The  solution  to  Eq.  (A.29)  can  best  be  handled  by  numerical  techniques  using  Simpson’s  rule  (see  (Krylov 
62] . pp  94).  Thus  we  can  express  P{Q  as: 


/j  + /2  , for  p2  > pj 
W)s<  0<k<  1 

|/j  , forp1  <p,  , 


(A.30) 


When  k = 0,  /j  = /2  = 0,  hence  P(C)  = 0.  From  Eq.  (3.33)  when  k - 1 ,P(B)  * 0,  since  the  limits  on  the 
integral  become  equal. 


t 
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APPENDIX  B 

ADAPTIVE  ROUTING  PROGRAM  FOR  COMPUTER  SIMULA  HON 

The  computer  simulation  for  adaptive  routing  was  initially  run  on  a DDP-24  small  general  purpose 
computer  built  by  Computer  Control  Company  (3C’s).  A CDC-3800  was  used  for  a more  expanded  look 
at  th'  effectiveness  of  the  selected  algorithms. 

The  framework  for  the  program  was  set  up  from  the  RAND  program  of  Boehm  and  Mobley  [B066] 
and  was  extensively  modified  to  allow  for  all  the  adaptive  routing  techniques  addressed  in  this  study.  The 
RAND  program  originally  had  no  queue  mg  (buffers)  or  PAL  (Preassigned  Linksj,  hence  messages  were 
generated  such  that  the  network  had  a constant  loading.  Without  buffers,  the  messages  generated  had  a 
dependency  on  the  current  status  ot  the  system;  therefore  message  volume  (MV)  would  not  generate  the 
same  total  number  of  messages.  Incorporating  the  buffers  at  each  node  helped  to  establish  a common 
basis  for  comparing  the  selected  routing  algorithms.  Finally,  the  RAND  program  had  to  use  the  “Hot- 
Po'.aio”  technique  since  there  was  no  queueing. 

B.l  FLOW  DIAGRAM  OF  ADAPTIVE  ROUTINE 

The  flow  diagram  for  the  simulation  routine  is  given  in  Figs.  (B.l a)  and  (B.l b).  The  program  assumes 
that  the  network  topology  and  minimum  delay  tables  arc  already  on  magnetic  tape.  There  were  two  pro- 
grams used  in  this  study,  RANSIM  and  ARPSIM.  Their  difference  was  only  in  the  manner  of  updating 
and  type  of  doctrine  used  (NDOCT).  The  Fortran  listing  of  the  ARPSIM  program  is  given  in  (B.2).  After 
some  prescribed  running  time  Tt  the  program  prints  out  the  following: 

1.  Total  messages 

2.  Undelivered  messages,  numbers  that  are  too  old 

3.  Number  of  messages  finding  queue  full 

4.  Delivered  messages 

5.  Average  time  for  message  in  network 

6.  Total  number  of  messages  handled  by  network,  and 

7.  Tabulation  of  routing  tables. 

The  RANSIM  program  had  enough  flexibility  to  try  various  updating  strategies  such  as  backwards 
learning,  negative  reinforcement,  etc.,  with  different  combinations  of  routing  procedures.  Both  programs 
could  incorporate  options  for  priority  assignment  of  messages,  servicing  of  queues,  and  a tag  on  the 
message  to  indicate  last  node  visited,  last  two  nodes  visited,  etc. 


Fig.  Bla  - Flow  Diagram  of  Program  Routine 
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Fig.  Bib  flow  Diagram  of  Program  Routine 
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B.2  PROGRAM  ARPSIM 

COMMON  I WORK ( 19, 1 9 ) ,NL( 19 ) ,LN ( 19 #6 ) 

COMMON  LOT ( 19*6) » L I D ( 19,6) » IROUTl (4) , NRtK ( 19) 

COMMON  I WR I Tc ( 2 1 » 30  ) , I RT 1 ( 1 v ,2 , 1 9 ) 

COMMON  IWARK(  19,19)  » IG(  19,50)  ,NLREM ( 19 ) 

COMMON  MESGI (A),NM(  19)  » I P ATH ( 1 9 , 50 ) ,NQ ( 19,6) 

COMMON/bLKl /AA ( 19, 10), A(  19) 

NMAX=19 
N L M A X = 6 
MAXMS9=3Q 
REWIND  2 
CALL  RANFCitT(b) 

2 CONTINUE 

READ  (60*1902)  (A(I),I=1,  NMAX ) 

CALL  rstart 
CALL  RANESET(S) 

READ (60  * 1901 ) MV ,M T F , MUTMAX *MDT ,KPLF,KPCL ,M I ND , NDEL T * NDOC T , INDPRM, 
1 IMIX*MMmX 

IF  (MUT.EU.O)  00  TO  170 

READ (60, 1901  ) NTR Y ,L0OP ,MTC ,MtbPcR , I DP 

ME  S PER  = ME  SPE  R*  MDT MAX / 1 00 

READ  (60,1901)  NMM 

WRI TL (61 ,1930)  L INC , NDOC  T 

W R I T t ( 6 1 , 1 9 b 2 ) M V , M T F » M D T M A X . M D T , k P L F , K P C L , M I N D , N D E L T , 1 M I X 
WRITE  (61,1902)  (All),  1=1,  NMAX) 
i IF  (UN  IT, 2)  3 *130 
130  Rt  AD  (2)  LlNCtNL  ,LN,LDT  ,L  ID,  I wAKf>. 

4 IF  (UNIT ,2)  4,160 
1 6 C REWIND  2 
NLF  = ^ 

NCL  = w 

DO  430  1=1, NMAX 
IF  ( NL ( I ) • EG • 0 ) GO  TO  430 
NLF=NLF+1 
DO  420  J= 1 » NLMAX 
IF  ( LID(  I , J ) .tQ, 1 NCL  = NCL+  1 
420  CONTINUE 
43  0 CONT  INUE 
NCL=NCL/2 
Rs-ANF ( - 1 ) 

AkN  = R 
AARN=ARN 

MLF  = ( NLF* ( 10G~kPLF I +39 )/ 100 

MCL=  ( NCL*  ( 100-KPCD+30)  /100 

NP0S=0 

NPOSA=0 

MRATE=0 

NDELT I=NDELT 

MREC=0 

Ml XR=0 

MX  = 0 

MY*0 

MZsO 

NZ  = 0 

DO  303  1=1,  NMAX 
NM ( I )=0 

DO  306  J=1 , NLMAX 
306  ND(  1 ,J)=0 

DO  303  J = l*  MAXMSO 
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I 0 ( I • J ) =0 
505  IPATH(I*J)«0 
KLF=NLF“MLr 

IF  (KLF.LE.O)  GO  TO  530 
DO  520  I«l#  KLF 
N*NLF*I“1 
R«kANF(-l) 

KlLL=MlN0(Nt P I ( N ) +1  • ) ) 

N“  0 

DO  510  J” ' * .x 

IF  (NL(J).LE*0)  00  TO  510 

N«N+1 

IF  (N.LT.KILL)  GO  TO  510 
NL( J)=-NL( J) 

GO  TO  520 
510  CONTINUE 
520  CONTINUE 
530  DO  580  1=1*  NMAX 

IF  <NL( I ) .GE.U)  GO  TO  580 
N=-NL ( i ) 

DO  570  J=l,  N 
LCLN  ( I *J)/100(J00 
KbMOD ( LN ( I .J) *100) 

L I D ( I * J ) =-L I U ( I iJ) 

570  LID(LtK)=-LlD(L,K) 

560  CONTINUE 

KCL=NCL-MCL 

IF  ( KC  L • LE • 0 ) GO  TO  711 
DO  710  1=1*  KCL 
KCLl=2#(NCL-i+l  ) 

R=KANF ( - 1 ) 

KlLL=MINOUCLI  * I NT  ( REFLOAT  ( KCL  I ) + l.)  ) 
DO  690  J=l*  NMAX 
NLN= IADS ( NL ( J) ) 

IF  (NLN.EQ.O)  GO  TO  690 
DO  680  < = 1*  NLN 
IF  (LID( J*K)  .nl.1 ) GO  TO  680 
K ILL=K ILL-1 

IF  (KlLL.GT.O)  GO  TO  660 
L ID  U*K  ) *-L  ID(  JtK  ) 

L-LNC J#K )/100000 
N = MOD(LN{J  * K ) *100  ) 

LIDIL*N)=-LID(L»N) 

GO  TO  695 
680  CONTINUE 
690  CONTINUE 
695  CONTINUE 

710  CONTINUE 

711  DO  720  1=1,  NMAX 
NLRLM ( I ) =0 

IF  (NL< I ) .LE.U)  GO  TO  720 
DO  718  j=i,nlma;: 

IF  ( LID(  I , J)  .to.  1 ) NLRIHI  I )=.*LRLi.(  I )+l 
718  CONTINUE 

IF  (NLKLM(I) .EO.O)  NL(I)=-NL(I) 

720  CONTINUE 
MCL  = o 
MLF  =u 

DO  725  1=1, NMAX 
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IF  (NL( I ) .LE*Q)  GO  TO  725 

MLF  =MLF+ 1 

DO  724  J=  1 * NLMAX 

IF  ( LlDII *J) #E0.1)  MCL=MCL+1 

724  CONTINUE 

725  CONTINUE 

DO  729  I«1#NMAX 
IF  ( NL ( I ) .LE.O)  GO  TO  729 
DO  727  J=1#NMAX 
I RT I ( J * 1 * I ) *0 
I RT I ( J * 2 * I ) =0 

IF  ( J • EC  I.OR.NL(J).LE.O)  GO  TO  727 

DO  726  K = 1 * NLMAX 

I R = t K + 2 ) /3 

L*LN  ( I *K  ) / 100000 

IF  (L.EO.OI  GO  TO  9000 

I RT I ( J • I k • I ) =M 1 NQ ( I WARM J *L)+LDT ( I *M  »999 ) +10CO* IK  I I l J*  IR* I ) 
GO  TO  726 

90u0  IRTItJ»lk»I)  =999  + 1000*  I R T I ( J*  IK*  I ) 

726  CONTINUE 

727  CONTINUE 
729  CONTINUE 

CAUL  IDEAL 


WR I TE ( 6 1 * 1 800 ) 

DO 

7115  I=1*NMAX 

IF 

(NL(I).LE.O)  GO 

TO 

7115 

Ks0 

DO 

7110  J=1*NLMAX 

IF 

i L I D < I *J) • N L • 1 ) 

GO 

TO  7110 

K*K+1 

IWRI  ItlKil )=LN(  I *J)/1QQ000 
7110  CONT INUE 

WRI TE(61*1801)  I * ( URITE(  J*1  ) * J=l*  K, 
7115  CONTINUE 

18uQ  FORMAT  (12H  NODE  LINKS) 

18 J 1 FORMAT  (1H  *I3*4X*3013) 

MT--MDTMAX 

7295  DO  8 1 u 1 = 1*  NMA  X 
NLN  = NL ( I ) 

DO  310  J=l*  NLN 

If-  ( I 0 ( I » J ) / luG  • tUo  0 ) GO  TO  810 
IQ(  I * J ) = I Q I I * J ) -MD  T * 1 00 
b i 0 CONTINUE 

LRN  = MDT:-.AX  + MT 
DO  88u  11=1*  NMAX 
NLN  = NL  UI) 

IF  (NLN.LE.O)  GO  TO  tt 80 

DO  8 70  J R N = 1 * NLN 

J = MOD ( LKN  + JKN  * NLN ) + 1 

If-  ( IO(  I I * J ) •lCJ.OJ  GU  10  8 70 

IF  ( lu(  I I * J ) /100.NL .0  ) GO  TU  870 

1 0 ( I I * J ) =0 

K.  K = L N ( I I *J  ) / 1O0000 

JJ=MAXMSG+f-‘OD(  LN(  I I * J ) • lOO)-NL(KK) 
L*M0D(  IO  ♦ JJ ) /100 jOOOOOO * 10'  )+l 
MESGI ( 3)-  jui Io(KK*JJ)/100^o*1000) 

I = KK 

IF  ( L.LU. I ) GO  TO  828 

IF  (Ml SGI ( 3) •LT.MdTMmX)  GO  Tu  850 
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826  IF  (MT.LE.O)  GO  TO  827 
MY»MY+1 
GC  TO  829 

828  IF  (MT.LE.O)  CO  TO  827 
MREC=MKt01 

M I XR=M I XR+ 1 
MRATE=MKATE+1 
MZ«MZ+MtSGI (3) 

829  MX*MX+1 

IF  ( iQ(KK»JJ)/iOOOOCOOOOOO.LQ.2 ) NPOS«NPC$+l 
e27  CONTINUt 

NQ(  i I *J)*NQ(  II  * J ) - 1 
IU(Kft»JJ)=0 

IF  (MT.LE.O)  GO  10  84U 
IF  (MOu(MX.Ii'.IX)  .NL.O)  GO  TO  84<J 

830  NPOSA=NPOSA+NPOS 

WRITE(  61 .1961 ) MXiMRF.C.NPOSA.  IMlX.MIXR.NPObtMl 
NPO$=C 
M I X R = 0 

IF  (MX.GL.MMAX)  GO  TO  930 
8 A 0 CONTINUt 
GO  TO  870 
860  N L N - N L ( I ) 

DO  rbt  .1  11  = 1*  NTRY 
NN* 32768 
KGU  T r . 

DC  763  ft=l.  NLN 
I R=  ( K ••  2 ' / 3 
I KX  - * ' >loD(K“l*3  ) 

A A ft  N = -•  '■  0 \ A A K h + « 3 • 1 • ) 

I ROUT  Li i )=MODl  I KT I ( L » I ft  * 1 ) / lUUO** I KX • 100U ) +NO l I • K ) 

IF  ( NN • L ; s’  T KOU T t ( 1 ) ) GU  10  763 
IF  ( L 1 U ( I »•  N t • 1 ) GO  TO  763 

IF  ( 1 I 1 .LCJ.l  i GO  TO  766 

IF  ( I I 1 • Gl  • 2 • A,, * LN  ( I • k ) / 1uuuuu.Lu.nl  ) Gu  i u 7 8 3 

IF  ( I I 1 . GL  - 3 . And  • li\  i I #k  ) / 1 uuuuu  . Lu  • m.l  ) Gu  iv  763 

766  CONTINUl 

IF  ( NN  • L(J.  1 KOJ  T L ( 1 ) • a.  iL.AARm#Li..6)  GO  IO  763 
NN= 1 ROUTE ( 1 ) 

K.0UT  = K 

763  CONTINUL 

IF  (KOUT.EO.Q)  GO  TO  826 
IF  ( LOOP • EQ. 0 ) GO  TO  766 
DO  766  IK=1»  LOOP 
KLL  = LN ( 1 tKOUT ) / 10000U 
IF  (nL.lu.nLL)  Gu  Tu  7t 6 

IF  (KLL. EU. MUO(  IPATHt  I»JJ)/lUl""l  In-1) »lul) ) Gu  lu  766 
766  CONTINUt 
GO  TO  766 
766  CONTINUt 

IF  (II 1 .EO*  ' ) K L = K L L 
IF  ( IlI.LO.NTftY)  GO  TO  826 
76A  CONTINUL 
766  DO  873  K=  1 ♦ NMh 
KL*K+NLN 

IF  ( lU(KK.KL).NL.O)  GO  10  873 
I0(KK*KL)=IQ(KK»JJ) /100*100+KOUT 
IO(KK*JJ)=0 

NO ( KK »KOU  T ) *NU(KK.KOUT ) + l 
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NQ ( I I >J )=NQ( I I t J ) -1 
IPATM(Mc.,kL)  = iHATH(KK*JJ) 

GO  TO  870 
873  CONTiNUL 

I Q ( I 1 t J * -iuO*MDT  + 2 

MESGI  < 3)=MIN0U*.0D(  iQ  ( KK*  JJ  ) / 10000  * lUOO  ) +MD  I *999) 

I Q ( KK*JJ)  = IOUK*  JJJ  + 1C000*  C.tSGl ( 3) -MODI  IQ(KK*JJ)  / 10000*1000 ) ) 
870  CONTINUL 
880  CONTINUL 

DU  78o  JJ=1*  NiiAx 

IF  ( NL  ( J. 1 ) * LL  • 0 ) uu  lu  / L»  U 

K*KANF (“1 ) 

CALL  POlS(RtJJtNMSG) 

AkN  = K 
AAR,\  = k 

IF  (NMGG.EU.O)  GO  TO  780 
KtSGl ( 2 ) = JJ 
DO  741  1=1*  NMSG 

R=RANF (-1 ) 

n*nlf- 1 

IF  (MIND.NE.O)  NLF  = .\LF-1 
N s M I N 0 ( N * I NT (REFLOAT (N5  + 1 . ) ) 

DO  748  J=l*  NMAX 
IF  ( J.t-Q.MESGI  ( 2 ) ) GO  TO  748 
IF  <NL(J)*GT*u)  n=n-1 
IF  ( NL  < J ) • L T • 0 « AND  • K I NL  • Nl  • (J  ) h = N— 1 
IF  (N.LL.O)  GO  TO  747 
7a8  CONTINUE 
747  MESGI ( 1 ) *J 
N“MESGi (2 ) 

7 A 8 NLN  = NL(m) 

do  781  ii  = i *n 
M * N L N+ 1 I 

IF  ( !u(N*M) *Nl.O)  GO  10  781 
I PATm  ( N *M  ) =!  L.SGI  (2) 

IUN  *M)  = 100w0v000U*J+lUUUUlol*;sfci>Gl  i 2 ) - 1 u 1 wuuu  1 u 1 
IF  ( I /.  0 1<  k ( N * J ) • N t • 3 2 7 6 8 ) i o * N » i « i u i r » f",  / t 2 u u u u u u u u u u u 
NLN-NL ( N ) 

NN* 32768 
DO  782  k=l*  NLN 
Ik* ( K+2 l / 3 
I nXs  2*viuD  I K“  1 * 3 ) 

AARN  = a,v,0D(  AARN+.3*1  • ) 

I ROUTE ( 1 )=KOD<  I KT 1 ( J * I K • N ) / 1000** I KX * 10  00 ) +NQ l N * K ) 

IF  ( NN.LT. I ROUTt ( 1 ) ) GO  TO  ' j2 
IF  ( L 1 D ( N * k ) • i*  L • 1 ) g^j  lu  f \ i. 

IF  ( NN.EO. i ROUTE! 1 ) .AND.AARN.L I • *5 ) GO  IU  782 
NN-IRQUTc ( 1 ) 

KOUT  =K 

782  CONTINUL 

1 0 ( N ».v. ) s I U ( N *M  ) / 1 00*  1 00  + i^uuT 
NO < N *KOUT ) -NO( N*kUU I ) +1 
GO  TO  741 
781  CONTINJL 

IF  (MT-GT.O#  NZ=NZ+1+NMSG“I 
GO  TO  780 
741  CONTINUL 
780  CONTINUL 

AARN=AM0D ( AARN  + R » 1 • ) 
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DO  881  1=1 i \MA  X 

IF  ( N L R t M ( I ) *LL.O ) GO  10  881 

NLN=NL ( I ) 

DO  882  K= 1 i NLN 
IF  ( I U ( I i K ) • i\  l • 0 ) 00  Tu  8 82 
Du  8 6 3 J = li  N i'M-i 
J J = J + hLh 

IF  (MODI  10  ( I i JJ  ) ♦ IOC  ) .tG.k ) 00  10  88A 
88  3 COMJMUl 
GO  TO  882 

88A  IF  ( M T 0 • 1 0 • 0 ) Oo  TO  892 
Ni\  = NL  N + uMM 
DO  891  NL-JJiNN 

I F ( MuD  ( I U ( I iKL  ) *100)  Ou  I vy  oV  1 

IF  (MODI  I0(  I iKL  )/ 10000  » 1000)  .LT.MLoPtR)  00  lo  891 
JJ  = KL 
GO  TO  892 
391  CONTINUE 
892  KOU  T = K. 

IQ  (I  ) =100*LdT  ( I ik  ) + l 
Mi-.  = LN  I I ) / 1 CD  QUO 
LL=/,UD(LN(  I *k)  * 100) 

f'.'N  = MAXt‘.  bG-NL  I * * » i • i ) A’LL 

MtSGI I 3 ) =MI NO  I MOD  I I0(  I * JJ )/ 10000 » 1000 ) +LD I (MMiLL)  »999) 

I P A Th  I MM  » NN ) = 1 0 1 *M0D  ( IPATmI  I iJJ)  »104Q6QA01  ) + I 

I 0 ( MM*NN  ) = I 0 1 I ♦ JJ  ) / lQOoGCQo-"-lUOUUOU  J+lUO*  I +MM- 1 0 1 + 1 U00U*ML  uG  I l 3 ) 

I C < I tJJ  ) =0 
663  CUNT  I Nut 
8 82  CON T I NO t 
886  KsimLN 

NN  aNMM+i\u\ 

JJ=NLN+1 

DO  887  J = J J *NN 

IF ( I0( I • J ) • EU. 0 ) 00  TO  887 

MtSGI  ( 3 ) =M  I NO  I MOD  I l0(  I * J ) / 1 UOUO  * 1 000 ) +MDT  *999) 

888  K s\+  1 

I Pm  T h i 1 y k ] = i HmTm  ( * * j / 

lot  I ♦ l\  ) = I 0 I I < J I / lUOuliUbU#  l^«ww'jijuu  + 1*1  tool  ( 3 ) #luwbu+i'iub  I I o \ I *J)  tluuuj 
1 ) 

IF IJ.LU.K)  GO  TO  88  7 
I Q ( 1 ♦ J i = C 
8 8 7 C 0 f T I N u L 
881  CONTlNUt 

IF  (MT .LT  *0  ) 00.  ro  92  3 
IF  ( NDtLT I • GT«  C ) OO  TO  922 
NDEL  f I =uDELT 

921  WRITE  I 61  1 1938)  MT*MKmTl 
MKATL=0 

922  CONTINUE 

NDEL  T I -*•  N D E L T i~ M u T 

923  CONTlNUt 

IF  (MT.Gt.MTF)  Ou  lu  930 
MT=MT+KdT 

IF  ( MOD (MTiMDT, mAX/A ) .Nt.O)  Ou  lu  729  3 

IF  (NDOCT.E0.Q)  oO  TO  7299 

DO  113  1= 1,  NMAX 

DO  112  J=li  NMAX 

IWARKI I .J)“32768 

DO  112  K = 1 * NLMAX 


OQ 


IF  ( LID{ I *K) cNt. 1 ) GO  TO  112 

IK=(K+2)/3 

IRX=2-M0D(K~1*3) 

I ROUTE ( 1 ) = MOD( IRTI l J t IR » I) / 1000** I RX# 1000 ) +NQ ( I sK) 

IF  ( I ROUTE  ( 1 ) *L7 • ! WARK ( I • J ) ) I W ARK  (I  » J ) = I ROUTE  ( 1 ) 

112  CONTINUE. 

113  IWARK(I#I)«0 

DO  124  I-l.  NMAX 
DO  124  J=l.  NMAX 
I RT I ( Jtltl )«0 
IRTI ( J.2.I ) =0 
DO  124  Kr  1 • NLMAX 
I R= ( K+2 ) / 3 

IF  (LlD(I#K)«Nt*l!  GO  TO  12b 
L = LN ( I *K ) / 1 00000 

IRTI  U*  IP,  • I ) = 1000*  IRT  I ( J » IR  1 1 ‘+M I NO (LOT  i I #K  ) + 1 WARK  ( L * J / ♦ I DP  * 999  ) 
GO  TO  124 

12  3 IRTI ( J • I R . I ) =1000* lKTI(J»lR*i)+999 
124  CONTINUt 
GO  TO  7295 
930  CONTINUt 

AMZ  = FLOAT  ( MZ ) /FLoAT * MREC ) 

9 50  M = MQD ( MX  * I M I X ) 

NPOSA=NPOSA+NFO$ 

WR  I T t ( 6 1 . 1 96 1 ) MX  *hHLC .NpUSA  *K»n I Xk . NPOG • M T 
WKl  rt  (61,1952)  MX.MY.KRt-C.AHZ 
WRITL (61 .1952)  NZ 
IF  ( iNDPRM.tU.O ) GO  TO  995 
NROW=NMA\+2 
DO  955  1=1.  NPOW 
DO  555  J= 1 .30 
955  I.  . I T l.  ( ! . J ) * 0 
N = 0 

DO  990  1=1. NMAX 

IF  (NLklM(  INLt.O)  GO  TO  972 

N = N«2 

I’fRITt  ( 1 *N-1  )=  I 
M = 2 

DO  957  J= 1 , NMAX 
F ( uLRcM ( J ) • L l • 0 ) GO  TO  95  7 
M = i'K  1 

I W R I T t ( M * N ) = J 
957  CONTINUt 

DO  9 7u  J= 1 • NLMAX 

IF  ( L I D ( I . J ) * N t • 1 ) GO  TO  970 

N*N+i 

XK  = 2 

ItoRl TL ( 2.N ) =LN( I .J)/100000 
I K=  ( J + 2 ) /3 
IRX=2-MOD( J-1.3) 

DO  c.  60  K = 1 • MM«X 

IF  (NLKLM(K ) *LL.0)  Gu  TU  96 C 

KK=KK+1 

IwRITt (KK.nJ=M0D( I K r I (k.Ik.I )/100Q*MRXtlQOO) 

960  CONTINUt 
970  CONTINUt 

972  IP  ( I • L T #NMAX*  ANU.N  + i^LKLM  ( I + 1)+2«LL«3Q)  GO  TO  990 
W R I T k I o 1 . 1 9 6 0 ) ( (Ii’.RITt(J.K)  • X = 1 * 3G)#J=1.  XK) 

N*C 
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DO  980  J=lt  NROW 
DO  980  K*1 *30 
980  I WR I TE I J*K ) *0 
990  CONTINUE 
995  CONTINUE 
60  TO  2 
170  STOP 

1901  FORMAT  (1216) 

1902  FORMAT ( 12F6*  3 ) 

1802  FORMAT  <10112) 

1950  FORMAT  (3H1*10X*18HLINK ^CONFIGURATION*! 8 *5X  *8HDOCTR INE  #18) 

1952  FORMAT  '17HOTOTAL  MESSAGES  *18/ 17H0UNDEU I VERED  - I 8/ 1 7H0DEL I VERE 

ID  *i3/17H0 AVER AGE  TIME  *F8.2) 

1953  FORMAT  (1X*I6*30H  MESSAGES  FOUND  THE  QUEUE  FULL  ) 

1955  FORMAT  < 3H0T* * 1 6 • 5X > 1 2HKESSA6ES/DT* * 1 6/ ) 

1960  FORMAT  ( 1H0 *30X * 14HR0UT I NG  TABLES/ / ( 30 1 4 ) ) 

1961  FORMAT  ( 2 1HOTOT AL  MSGS  EXPIRED  = 1 5 • 5X * 1 6HT0T AL  MSGS  DLV  =I5*5X*10H 

1P0SS IdLt  = I 5/2 1H  OF  LAST  I5*5X*16H  MSGS  DLV  *15*5 

2X  * lOHPOSS I PLE  =I5/4H  T *15) 

1962  FORMAT  <7H0MV  I4*5X*6HMTF  1 4 . 5X * 6HMDTMAX 14  * 5X *6HMDT  I4»5X»6 

1HXPLF  1 4 • 5X  * 6HKPCL  l4/7ri  MIND  1 4 * 5X  * 6HNDELT  I4*5X*6HIM1X  14) 

END 


SUBROUTINE  RST ART 

COMMON /BLK 1 / AA ( 19*10) *A  < 19) 

NMAX* 1 9 

DO  2 uO  1 = 1 > NMAX 
AA ( I .1 )=EXP(-A( I ) ) 

T* A A ( 1*1) 

DO  200  J = 2 • 10 
T J = J-l 
T = T* A ( I ) /TJ 

200  AA ( I * J ) =AA ( I *J-1 )+T 
RETURN 
END 


SUBROUTINE  PulS  (R*I*NMSG) 
COMMON /bLKl /AA { 19*10) *A< 19) 
DO  2oO  J=l*10 

IF  IR.LT.AAI  1 * J ) ) GO  TO  100 
200  CONTINUE 
NMSG=1Q 
RETURN 
100  NMSG= J~ 1 
RETURN 
END 


SUBROUTINE  I Dt AL 

COMMON  I WORK,  l 1 9 » 1 9 ) *NL'i9)  *LN<lv*6) 

COMMON  LDT  < 19*6)  »LIL(\9*6)  » I ROUTE  I 4 ) *NREM ( 1 9 ) 
DIMENSION  NC ( 1 9 ) 

NMA:  -19 

DO  5oO  I = 1 *NMAX 


DO  20  Js 1 tNMAX 
I WORK ( ItJ)*32768 
' CU)=0 
2t-  CONTINUE 

IF  (NL( I ) .LE.G ) GO  TO  500 
NODE' = I 
IDT -0 

50  CONTINUE 
N»NL (NODE ) 

DO  100  J«ltN 

IF  ( LlD  ( NODE  t J ) .NE • 1 ) GO  TO  100 
L*LN (NODE  tj)/100000 

I WORK ( 1 »L)*MIN0( I WORK ( I iL) *LDT ( NODE # J ) + \ DT ) 

100  CONTINUE 
NC ( NODt ) - 1 
IDT=32766 
DO  200  J = l>  NMAX 

IF  ( (N-(  J)  ./;l.O)  .UK.  ( 1WOKM  i *J)  .EU.32768)  ) 00  TO  200 
IF  (IWORKCiJ).GT.lOT)  GO  TO  200 
NODE  = J 

IDT=  1 WORK  I I tJ) 

200  CONTINUE 

IF  ( IDT. NE. 32768)  GO  TO  50 
I WORK ( I 1 1 ) =0 
500  CONTINUE 
RETURN 
END 


B.3  EXPLANATION  OF  VARIABLES 

Suppose  we  are  given  a communication  network  17  with  nodes  J / j , links  j iik  | where  iik  is  the 
Ath  link  out  of  node  / and  t.k  is  the  time  to  transmit  along  the  klh  link  out  of  node  /.  Each  node  is 
assumed  to  have  a storage  area  for  messages  and  the  assumption  is  made  that  if  there  is  a link  of  length 
from  nod  ; / to  node  /,  then  there  is  also  a similar  link  from  node  j to  node  i. 

The  program  ARPSIM  was  designed  to  run  a simulation  of  the  above  communication  network,  with 
and  without  damage  under  several  methods  of  routing  procedures.  To  aid  in  the  understanding  of  this 
program,  the  following  list  of  variables  is  given. 


NMAX 
N LAI  AX 
IWORK(l,J) 
NL{D 
LN(I,  K) 


LDTJ.K) 

NRh'M(I) 

NLREM(I) 

IWRfTE(I,J) 

NMM 

MAXMSG 

MV 


MZ 

NZ 

AMZ 

MT 

IRTliJ,  IR,  f) 


- maximum  number  of  nodes  in  the  undamaged  network 

- maximum  number  of  links  which  can  come  out  of  any  node 

- minimum  path  length  from  node  / to  node  J 

- number  of  links  out  of  node  / in  undamaged  network 

= 100000  V + L where  J is  the  end  node  of  the  A'th  link  out  of  node  / and  the 
backward  link  is  the  Ith  link  out  of  node  J 

- the  length  of  the  A'th  link  out  of  node  / 

- number  of  links  out  of  node  / in  the  undamaged  network 

- number  of  links  out  of  node  I in  the  damaged  network,  in  particular  if 
NLREM{I)  < 0,  node  / is  isolated 

- is  used  for  outputting  data  only 

- number  of  storage  locations  in  each  queue  available  for  messages  waiting  to  be 
transmitted 

- NMM  + 2*  MM  AX 

- percent  of  network  loading.  The  number  of  messages  generated  in  each  time 
interval  is  equal  to  M V*  (number  of  links  in  network)/ 1 00  (used  only  for 
constant  input  ) 

number  of  delivered  messages 

number  of  generated  messages  finding  the  storage  queue  full 
number  of  messages  it  war,  possible  to  deliver 

- accumulated  simulation  time 

- contains  the  routing  tables  from  node  / to  node  J,  in  particular  the  time  to  go 
from  node  I to  node  ./  using  the  A'th  link  out  of  node  / is  given  by  ( IRTJ(J,1R , 
/ >/ 1 00  • • IRX),  modulo  100  where  IRX  - 2 - Kmodl  and  IR  = (A  + 2)/3 
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IQ(l,  50)  - the  storage  queue  for  node  /.  IO(I.  50)  contains  three  parts 

(Part  H for  A'  = 1 SLMAX,  IQ{I,  fC)  = 100V  + 1 indicates  the  status  of  the 

Ath  link  out  of  node  /.  If  IQ(I,  K)  F 0 the  link  is  busy.  J is  the  remaining 
time  required  to  complete  transmission  of  the  message. 

(Part  IP,  for  A = SLMAX  + 1 , NLMAX  + NMM,  IQU,  A)  is  the  storage  area  for 
messages  waiting  to  be  transmitted;  if  IQU,K)  = 0 this  space  is  empty. 

(Part  III)  for  K = 50  - NLMAX  t 1 50,  IQ(L  A)  contains  the  message  being 

transmitted  from  the  node  L = LN(I,  Ay/ 100000  along  link  number  LNfJ,  K) 
mod lOO  to  node  /.  If  1Q(I,  K)  = 0,  there  is  no  message  being  transmitted 
along  this  link. 

JPATIHL  A)  contains  the  past  history  of  the  message  presently  at  IQ(I,  A').  IQU,  A')/l  \ **IRX 

modulo  10  is  the  (IRX  + 1)  previous  node  visited. 

A{1)  te  average  number  of  messages  arriving  at  node  / during  a time  interval 


The  actual  message  M being  transmitted  is  the  following: 


Mmod  JO2  + 1 

Mmod\Qr  + 1 

3//102modulo  102  + 1 
A/./103  modulo  103 
3// 107  modulo  102  + 1 
Ml  I09  modulo  102  + i 

Mj  10 11  modulo  10 


is  the  present  node  of  the  message  M if  the  message  is  in  Part  III  of  the 
queue  or  has  not  been  placed  in  the  queue 

contains  the  number  link  over  which  the  message  will  be  transmitted  if 

the  message  is  on  Part  11  of  the  queue 

is  the  last  node  visited  by  the  message 

is  the  time  the  message  lias  been  in  the  network 

is  the  source  node  of  the  message 

is  the  destination  node  of  the  message 

if  this  is  even,  message  is  deliverable;  if  odd,  undeliverable 

if  this  is  greater  than  1 the  next  node,  the  message  is  going  to  is  the 

san^  as  the  !act  n • le  visited 


MTI- 


MDTMAX 


MDT 


KPLF 

KPCL 


maximum  simulation  time 

warm-up  time  for  the  system;  maximum  time  a message  can  remain  in  the  network 
before  being  considered  unuelivcrable 

length  of  simulation  time  unit;  must  divide  LDT[I,K)  evenly  for  all  nodes  I and 
links  K 

percent  of  nodes  to  be  disabled 
percent  of  links  to  be  disabled 
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MIND 


NDELT 

NDOCT 


ISDPRM 

MM  AX 

LOOP 

XTRY 


MTC 

MESPER 

MX 

MY 


RSTART 

RANDOM 

IDEAL 


- a>-  indicator  to  determine  whether  or  not  generate  - messages  can  have  disabled  nodes 
as  if  stinations.  If  MIND  = 0,  we  can  use  inopera  )le  nodes  as  destinations,  if  MIND 
--  1 we  cannot. 

- time  increment  between  summary  prints 

- indicates  the  method  of  routing  table  updating  to  be  used.  If  NDOCT 


RANSIM 

AliPSIM 

= 1 

no  updating 

= 0 ii  > updating 

= 1 

backwards  learning 

^ C update  every  A/DFmax 

= 3 

negative  reinforcement 

= 4 

superposition  of  (2)  and  (3) 

= 5 

bi-adaptive 

= 6 

dynamic  programming 

- indicator  for  printing  routing  matrices.  If  INDPR V - 1 we  print  matrices;  if  INDPRM 
= 0 we  do  not. 

maximum  number  of  messages  to  be  delivered  before  simulation  is  stopped 

- number  of  previous  nodes  that  a message  will  have  remembered  visiting 

number  of  links  out  of  a node  that  a message  will  try  in  order  to  find  a node  differ- 
ent from  a previous  node;  it  the  search  is  unsuccessful  the  message  is  considered 
undeliverable 

--  indicator  as  to  whether  or  not  old  messages  should  be  given  a priority;  if  MTC  - 0 
no  priority  is  given,  if  MTC  = 1 a priority  is  given 

- on  input  percent  of  A IDT^^  which  determines  when  a message  is  old;  later  it  is 
changed  t > the  age  when  a message  is  considered  old 

- total  number  of  messages  handled  by  network 

- number  of  undelivered  messages 

SUBROUTINES 

- initializes  laiidom  number  generator 

- generates  a random  number  in  [0.,  .999] 

- calculates  for  a given  communication  network  the  minimum  path  lengths  between 


nodes 
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B4.1  SIMULATION  RUNS  FOR  TH£  DDP-24  (No  Destruction  of  Networks) 

The  following  list  comprises  the  various  routing  algorithms  used  on  the  DDP-24  computer.  All  runs  were 
made  on  an  8-node  network  with  40  links  fully  duplexed  and  no  network  damage.  Runs  were  made  for  identi- 
cal net  loadiigs  under  the  following  constraint": 

1.  Messages  were  generated  from  a uniform  distribution 

2.  Single  destination  for  each  message;  all  destinations  equally  likely 

3.  Finite  storage  opacity  at  each  node 

4.  Defections  from  the  system  were  allowed  as  network  loading  increased 

5.  All  messages  were  of  equal  length 

6.  Servicing  was  on  a first-come-first-served  basis  (FCFS) 

7.  Service  distribution  was  a constant 

8.  Perfect  transmission  and  reception  assumed 

The  program  algorithms  are  defined  along  with  their  abbreviations  as  follows: 

a.  Dynamic  programming  (DY) 

b.  Basic  program,  no  updating  (NU) 

c.  Backward  learning  (BL) 

d.  Negative  reinforcement  (NR) 

e.  Superposition  of  BL  and  NR,  where  BL  can  increase  or  decrease  table  entries  (SP) 

f.  Last  M nodes  visited  (LMNV) 

g.  Best  three  (output)  links  (BTL) 

h.  The  combination  of  LMNV  and  BTL  will  be  further  defined  as  ( MOD M)  where  M = number  of 
past  nodes  visited 

i.  Aging  parameter  (X%  aging  on  queue) 

Table  B1  lists  the  computer  runs  for  the  8-ncde  network  configuration  (see  Fig.  4.1). 


TABLE  B1 

COMPUTER  RUNS  ON  AN  8-NODE  40  LINK  NETWORK  WITH  NO  DESTRUCTION 
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1. 

NU 

12. 

NU,  LNV,  BTL,  10%  aging 

2. 

BL 

13. 

NU,  LNV,  BTL,  50%  aging 

3. 

SR 

14. 

i 

NU,  LNV,  BTL,  157c  aging 

4. 

SP 

15. 

BL,  LNV,  BTL,  10%  aging 

5. 

BA 

16. 

BL,  LNV,  BTL,  50%  aging 

6. 

SU,  LNV,  BTL 

17. 

BL,  LNV,  BTL,  75%  aging 

7. 

NU,  L3SV,  BTL 

18. 

BL,  A 3 AT,  BTL,  10%  aging 

8. 

NU.  t.SNV,  BTL 

19. 

BL,  A 3 AT,  BTL,  50%  aging 

9. 

Fl,  LNV,  BTL 

20. 

BL,  A 3 AT,  BTL,  75%  aging 

10. 

BL,  A 3 AT,  BTL 

21. 

ARP  A,  UP  - 0 

11. 

BL,  L5NV,  BTL 

22. 

ARPA,  UP  = 0,  A AT,  BTL,  50% 

B4.2  SIMULATION  RUNS  FOR  THE  DDP-24  (With  50%  Destruction  on  Links,  No  Node  Damage) 

The  algorithms  listed  in  Table  B2  were  run  on  the  8 network  of  section  B1  with  50%  link  deduction  (see 
Fig.  4.2).  The  definitions  of  the  algorithms  and  constraints  were  identical  with  section  B4.1. 

TABLE  B2 

COMPUTER  RUNS  ON  AN  8-NODE  NETWORK  WITH  50%  LINK  DESCRUCTION 


1. 

DY 

16. 

NU,  LNV  BTL,  10%  (aging 

2. 

NU 

on  que) 

3. 

BL 

17. 

NU,  LNV,  BTL,  50% 

4. 

NR 

18. 

NU,  LNV,  BTL,  75% 

5. 

SP 

19. 

BL,  LNV,  BTL,  10% 

6. 

BA 

20. 

BL,  LNV,  BTL,  50% 

7. 

NU,  LNV,  BTL 

21. 

BL,  LNV . BTL,  75% 

8. 

NU,  A 3 AT,  BTL 

22. 

BA,  BTL,  10% 

9. 

NU,  L5NV,  BTL 

23. 

BA,  BTL,  50% 

10. 

BL,  LNV,  BTL 

24. 

BA,  BTL,  75% 

n. 

BL,  A3AT,  BTL 

25. 

ARPA,  NU 

12. 

BL,  lSNV,  BTL 

26. 

ARPA  (standard  with  updating) 

13. 

BA,  LNV,  BTL 

27. 

ARP\  i^NV,  BTL,  50% 

14. 

BA,  A 3 AT,  BTL 

28. 

ARPA,  DP*  = 1 

15. 

BA,  LSNV,  BTL 

29. 

ARPA,  DP  = 1,  A AT,  BTL,  50% 

* DP  = a bias  term  hat  is  added  to  the  routing  delays  to  prevent  looping 
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APPENDIX  C 

SIMULATION  ON  THE  CDC  3800 

C.l  SIMULATION  RUNS  FOR  THE  CDC  3800  (Nr  Destruction  of  Network) 

All  runs  were  made  for  identical  network  message  intensities  on  a 19  Node  ARPA  Network.  The  algorithm 
abbreviations  and  program  constraints  are  the  same  a.  given  in  Appendix  B. 

Runs  C2.1  and  C2.2  were  made  with  messages  generated  from  a uniform  distribution.  Runs  C2.3  and 
C2.4  were  made  with  Poisson  input  distributions,  such  that  the  average  message  rate  per  node  was  A^.  The 
total  message  rate  into  the  network  per  time  was  19  \.  (when  all  19  nodes  were  operable). 

The  periodic  updating  (PUD)  was  employed  on  some  of  the  runs  to  make  the  ARPA  routine  conform  more 
closely  to  its  present  configuration  and  a bias  term  of  DP=2  was  used. 

After  some  preliminary  runs  were  made  to  see  the  extent  of  link  or  node  destruction  that  could  be  used 
without  disconnecting  the  network  these  percentages  were  used  for  setting  up  Tables  C2.3  and  C2.4.  The 
message  intensities  were  also  adjusted  to  prevent  net  saturation  for  any  given  set  of  runs. 

TABLE  C2.1 

COMPUTER  RUNS  ON  THE  19  NODE  ARPA  NET  WITH  UNIFORM 
MESSAGE  INPUT  DISTRIBUTION 

1.  NU 

2.  BL 

3.  NU,  5 0%  aging  on  queue 

4.  NU,  LNV,  BTL,  50%  aging 

5.  ARPA,  NU 

6.  ARPA,  DP=3 

7.  ARPA,  DP=0,  50%  aging 

8.  ARPA,  DP=3,  50%  aging 

9.  A RPR,  LNV,  DP=3,  50%  aging 


TABLE  C2.2 

COMPUTER  RUNS  ON  THE  19  NODE  ARPA  NET  WITH  POISSON 
MESSAGE  INPUT  DISTRIBUTION 

1.  NU 

2.  BL 

3.  NU,  50%  aging 

4.  NU,  LNV,  BTL,  50%  aging 

5.  ARPA,  NU,  DP=0 

6.  ARPA,  DP=3 

7.  ARPA,  DP=0,  50%  aging 

8.  ARPA,  DP=3,  50%  aging 

9.  ARPA,  LNV,  DP=3,  50%>  aging 

10  ARPA,  LNV,  BTL,  DP=3,  40%  aging 

1 1 . ARPA,  LNV,  BTL,  DP=3,  60%  aging 

1 2.  ARPA,  PUD,  LNV,  BTL,  DP=3,  50%  aging 

TABLE  C7.3 

COMPUTER  RUNS  ON  19  NODE  ARPA  NET  WITH  10%  LINK  DESTRUCTION 
(POISSON  INPUT  DISTRIBUTION) 

1.  DY 

2.  BA 

3.  BA,  LNV.  BTL,  50%  aging 

4.  DY.  LNV,  BTL,  50%  aging 

5.  ARPA 

6.  ARPA,  50%  aging 

7.  ARPA,  LNV,  BTL,  50%  aging 

TABLE  C2.4 

COMPUTER  RUNS  ON  19  NODE  ARPA  NET  WITH  1 NODE  DISABLED 
(POISSON  INPUT  DISTRIBUTION) 

1.  ARPA,  DP=2,  PUD 

2.  ARPA,  DP=2,  PUD,  509?  aging 

3.  ARPA,  DP-2,  PUD,  LNV,  BTL.  50%  aging 
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